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SUMMARY 


Consider  the  constrained  part  (denoted  as  the  Plant)  of  a 
control  system,  consisting  of  n cascaded  sections,  each  of  whose 
outputs  can  be  sensed  for  feedbek  purposes.  Feedback  from  these 
points  is  to  be  used  to  achieve  apriori  specified  tolerances  on  the 
system  response,  despite  great  uncertainty  in  plant  parameter 
values.  In  this  first  quantitative  work  of  its  kind,  the  feedback 
is  permitted  to  proceed  directly  to  internal  plant  variables, 
constituting  rplant  modification*.  The  reason  is  that  the  Internal 
plant  signal  levels  needed  to  achieve  a specific  output  are  now 
affected  by  the  feedback  loops.  This  is  in  contrast  to  previous 
quantitative  research  in  which  the  feedback  was  confined  to  the 
plant  input.  The  plant  signal  levels  needed  to  achieve  a specific 
output  were  then  not  influenced  by  the  feedback. 

Plant  modification  feedback  permits  greater  reductTbu  in  the 
"cost  of  feedback",  in  terms  of  feedback  loop  bandwidths  and 
effect  of  sensor  noise,  at  the  cost  of  increase  in  plant  internal 
signal  levels.  In  this  work,  the  maximum  permitted  increase  Q 
in  signal  level,  is  part  of  the  design  specifications.  A step  by 
step  design  procedure  is  presented  for  satisfying  this  Q require- 
ment and  the  system  response  tolerances  over  the  given  range  of 
plant  uncertainty,  and  doing  so  at  sensibly  minimum  "cost  of 
feedback".  This  permits  the  designer  to  achieve  desired  trade-off 
between  increased  plant  signal  level  and  cost  of  feedback. 
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STATEMENT  OF  PROBLEM  AND  PRELIMINARY  BACKGROUND 

1.1  Introduction 

This  work  deals  with  the  problem  of  making  a system  perform 
satisfactorily  despite  uncertainty,  in  the  following  context: 

Equipment  has  been  assembled  by  specialists  in  the  area  of 
concern,  in  order  to  achieve  certain  objectives,  for  example,  an 
airframe  and  engine  for  accomplishment  of  certain  aeronautical 
objectives,  a chemical  plant  for  production  of  certain  chemicals, 
etc.  This  assembly  of  equipment  is  denoted  as  the  plant.  The  plant 
has  the  ability,  the  muscle  so  to  speak,  to  achieve  the  objectives. 
However,  it  does  not  have  the  accuracy  needed.  This  is  manifested 
by  uncertainty  in  the  parameters  of  the  mathematical  relations  des- 
cribing the  plant.  For  example  in  Figure  1.1,  suppose  the  relation 


Plant 


Figure  1-1.  Plant 

between  plant  output  y and  input  x is  given  by  a linear  time- 


invariant  (tti)  ordinary  differential  equation 
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(1.1-1) 


i.e.,  the  A^  , B^  are  constants.  There  are  q physical  para- 
meters k^,...,kg  and  K , A^  , B^  are  functions  of  these  para- 
meters. The  values  of  these  parameters  are  not  known  precisely,  but 

1 2 

it  is  known  that  they  lie  within  certain  bounds  k.  € [k.,k.]  . 

ill 

Hence,  each  possible  parameter  vector  k , each  combination  of  k^ 
values  i = l,...,q  , gives  a different  plant  transfer  function, 
generating  a set  V = (P(s) } of  possible  plant  transfer  functions. 

Such  a formulation  of  the  uncertainty  problem  may  appear  naive 
because  one  might  argue  that  often  the  parameter  values  change  with 
time-giving  uncertain  linear  time-varying  relations,  because  the 
rate  of  variation  is  uncertain.  Also,  the  S-ti  description  is 
usually  an  approximation  of  a nonlinear  relation.  We  are  really 
assuming  iti  relations  with  the  above  uncertainty  form,  in  order 
to  be  able  to  rigorously  use  Laplace  transforms  and  frequency  res- 
ponse methods.  However,  it  has  been  rigorously  proven  [16]  that 
uncertain  linear  time-varying  plant  problems  are  reducible  to  the 
above  iti  uncertainty  form,  and  even  uncertain  nonlinear  time- 
varying  plant  to  a certain  extent  [15],  can  be  so  reduced.  Hence 
the  above  modelling  is  highly  justifiable.  Finally,  one  must  begin 
somewhere  with  the  development  of  a scientific  synthesis  theory  for 
uncertain  systems,  and  the  ftti  case  is  obviously  where  to  start. 

The  objective  is  to  achieve  certain  apriori  specified  per- 
formance objectives  V P € *P  . Since  the  overall  system  is  to  be 


Iti  , it  can  be  characterized  by  its  response  to  any  input,  and  the 
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. 


step  response  is  very  popular  because  it  combines  within  it  both  the 
fastest  kind  of  input  (an  abrupt  change)  and  the  slowest  (no  change) . 
Time  domain  specifications  are  reasonable  in  many  cases,  as  in 
Figure  1.2a,  the  step  response  is  to  be  inside  the  bounds 
blfb2ype1P  , with  additional  bounds  of  similar  nature  on  the 


C (t) 


lower  bound 
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Figure  1-2 (a).  Time  domain  step  response  specification. 

first  and  perhaps  higher  derivatives.  Our  design  technique  is  in 
the  frequency-domain,  so  we  must  translate  such  t-domain  bounds  into 
"equivalent"  w-domain  bounds  on  the  system  frequency  response 
T(ju>)  . If  the  system  is  minimum-phase  [2],  | T ( jto)  | suffices  and 

we  restrict  ourselves  here  to  such  systems.  This  translation  is,  as 
of  this  date,  an  engineering  art  rather  than  a science.  Advice  on 
how  to  translate  is  scattered  in  the  literature  [2,6,14].  Very  good 
results  have  been  obtained  with  only  moderate  effort.  We  shall 
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Figure  1-2 (b).  Frequency  domain  specification. 


assume  in  this  work  this  translation  has  already  been  done.  It  is 
worth  noting  that  it  has  been  shown  [15]  that  for  minimum-phase 
systems,  time-domain  specifications  on  the  step  response  and  on  its 
derivative  of  the  following  nature 

(t)  $ c(i)(t)  $ b*i}(t)  , i - 0,1,. ..,n  (1.1-2) 

can  always  be  satisfied  by  means  of  u>-domain  bounds  of  the  following 
nature 

B2(w)  { |c(jw)/R(jw) | $ B (w)  . (1.1-3) 


•L  • 
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In  our  work,  the  bounds  on  system  performance  will  have  this  form. 


s 

1.1.1.  Previous  work 

The  quantitative  aspect  of  our  work  cannot  be  over-emphasized. 
The  sensitivity  reduction  capability  of  feedback  is  very  well-known. 
Hundreds  of  books  and  thousands  of  papers  have  been  written  on  the 
subject,  but  the  number  of  these  which  are  quantitative  in  nature  is 
extremely  small,  i.e.,  with  uncertainty  bounds  and  performance 
bounds  explicitly  included  in  the  problem  statement.  It  is  as  if 
the  mere  use  of  a feedback  configuration  around  the  uncertain  plant, 
suffices  to  scare  it  into  docile  behavior.  In  the  vast  majority  of 
the  techniques  the  uncertainty  is  completely  ignored,  and  there  are 
no  or  extremely  crude  performance  specifications.  One  presumably 
emerges  with  the  same  design  whether  the  parameter  uncertainty  is 
x%  or  1000x%,  and  irrespective  of  whether  the  bounds  B1  , B2  in 
Figure  1-2 (b)  are  narrow  or  wide  apart.  There  is  no  concern  with 
the  ’cost  of  feedback'  - which,  aside  from  the  sensors,  lies  in  the 
bandwidth  of  the  loop  transfer  function,  and  little  concern  with  the 
extremely  important  matter  of  sensor  noise  [see  Sec.  1.21  . These 
points  have  been  emphasized  in  [5] . 

Our  work  follows  closely  in  the  tradition  of  'quantitative 
synthesis'  recently  established  [2,5,8].  To  appreciate  the  present 
work,  it  is  important  to  be  aware  of  the  highlights  of  this 
previous  work.  Quantitative  synthesis  was  first  developed  for  a 
plant  with  only  one  variable,  the  plant  output,  c(t)  in 
Figure  1-3,  available  for  feedback  [2] . The  system  command  input 
r(t)  was  also  assumed  accessible,  so  the  processing  of  these  two 
signals  provides  two  independent  compensation  functions  to  the 
designer.  An  infinitude  of  canonical  two-degree-of- freedom 
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Figure  1-3.  Structures  of  2-n.O.F.  system. 


structure  may  be  used  [4] . The  design  procedure  developed  in  [2] 
used  Figure  1-3 (a),  but  suppose  the  sensor  transfer  function  is 
H (s)  , then  one  can  use  Figure  1-3  (b),  letting  G-^H  (of 
Figure  1-3 (b) ) ■ G (of  Figure  1-3 (a)),  in  order  to  have  the  same 
loop  transmission  function  L(s)  = GP  - G^H  , and  F1G1  = FG  in 
order  to  have  the  same  system  transfer  function 


T(s) 


FGP 

1+GP 


F1G1P 

1+G.jPH 


(1.1-4) 


1.1.2.  A 2-degree-of- freedom  structure  with  2-loop 
implementation 

Suppose  large  loop  feedback  bandwidth  is  needed  and  it  is 
found  that  an  independent  sensor  measuring  c(t)  (e.g.,  a tacho- 
meter in  a position  servo)  gives  less  noise  than  the  differentiation 
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of  a position  sensor,  so  both  sensors  are  used,  as  in  Figure  1-4, 
with  the  two  sensor  transfer  functions  , and  say  the 

structure  in  Figure  1-4  is  used.  This  is  a two-loop  structure 


/ 

Figure  1-4.  2-loop,  2-D.O.F.  structure. 


physically,  but  in  terms  of  fundamental  feedback  design  it  is  a two- 
degree-of- freedom  system,  so  the  quantitative  design  theory  of 
Figure  l-3(a)  is  used,  giving  G and  F . It  is  required  that  the 
loop  transmission  around  P , be  the  same  in  both  cases,  i.e., 

L = PG  (Fig.l-3(a)  ) =*  P GJHj+G^)  (1.1-5) 

and 

T " F l+GP  (Fl9-1_3(a))  “ Fb  l+p[Ga(H1+GbH2) J (1.1-6) 

SO 

G - Ga(H1+GbH2)  , FG  - FbGbGa  . (1.1-7) 

Hl»  H2  are  *tnown»  so  one  must  decide  how  to  split  G = Ga^Hl+GbH2^ 
between  Gfl  and  G^  . This  is  done  by  considering  the  effect  of 


given  that 


• | ■ y. 
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is 


G = |(jw)  = Ga<H1+H2Gb>  (1-1-9) 

is  fixed  by  the  quantitative  design  technique  of  [2] . 

00  2 

The  objective  is  to  minimize  / 1x^1  dw  , subject  to  the 

o 

above  constraint.  This  is  a straightforward  optimization  problem 
which  can  be  solved  outside  the  realm  of  quantitative  feedback 
synthesis.  The  latter  only  provides  the  design  with  the  feedback 
loop  transmission  (L)  needed  around  the  plant,  and  the  prefilter 
(F)  needed  to  process  the  command  input  r(t)  . The  state-of-the- 
art  in  sensors  and  in  filter  synthesis  determines  how  L and  F 
are  to  be  realized.  In  fact,  in  the  above  context  one  might  con- 
sider use  of  an  accelerometer  in  a 3-loop  feedback  structure.  But 
from  our  point  of  view  the  structure  remains  that  of  a two-degree- 
of- freedom  system  and  we  shall  continue  to  associate  the  latter 
with  a single-loop  system. 


: 


1.2  Review  of  Two-Degree-of-Freedom  Quantitative  Design 
Theory 

GP 

Figure  1-3 (a)  is  used  with  T = F . It  is  assumed  that 

the  compensation  network,  whose  power  level  can  be  very  low  (as 
the  plant  contains  the  power  elements) , can  be  constructed  with 
negligible  uncertainty  in  their  transfer  functions.  Hence,  due  to 
the  uncertainty  in  P , 

CP  T 

4 in  T » A tn  A £n  , L - GP  . (1.2-1) 

and 

a in  | T ( jc*»)  I - A in 


L(jo>) 


1+L( ju) 


(1-2-2) 
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Given  that  A in  | T ( ju»)  | ( db  for  example  at  in 

Figure  1-2 (b) , what  are  the  resulting  constraints  on  L(jw)  ? It 
is  convenient  to  pick  a "nominal"  plant  Po(s)  • and  derive  the 
bounds  on  the  resulting  "nominal"  loop  function  Lq  **  P^G  . These 
bounds  can  be  found  by  means  of  a digital  computer,  but  it  is  very 
useful  for  insight  to  see  it  done  on  the  Nichols*  chart 
(logarithmic  complex  plane  with  abscissa  in  degrees,  ordinate  in 
decibels  • 20  log^)  . The  procedure  is  illustrated  for  the  case 

V a 

P(s)  = ~s"('s+a)  ' 1 5 k f 10  ; 1 J a $ 10  . (1.2-3) 

At  a)  = 2 rps  , P(j2)  lies  within  the  boundaries  given  by  A BCD  in 

Figure  1-5.  Since  (n  L ■ In  G + In  P , the  pattern  outlined  by 


Figure  1-5.  Derivation  of  bounds  on  L(jw)  on  Nichols  chart. 
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ABCD  may  be  translated,  but  not  rotated,  on  the  Nichols'  chart, 
the  amount  of  translation  being  given  by  the  value  of  £n  G(j2)  . 
For  example,  if  a trial  design  of  L(j2)  corresponds  to  the  tem- 
plate P(j2)  at  A'B'C'D'  in  Figure  1-5,  then 

lG«J>ldb  ■ lL(i2)lab  - [p<i2> ldb 


= (-2.0)  - (-13.0)  = 11.0  db 


(1.2-4) 


Arg  G(j2)  = Arg  L(j2)  - Arg  P(j2) 

= (-60°)  - (-153.4°)  = 93.4° 


(1.2-5) 


1.2.1.  Bounds  on  L(j<*>)  in  the  Nichols'  chart. 

The  templates  of  P(jw)  are  manipulated  to  find  the  position 
of  L( ju)  which  results  in  the  specifications  of  Figure  1-2 (b)  on 
fcn  | T ( ju>)  | being  satisfied.  Taking  the  w = 2 template,  one 
tries,  for  example,  positioning  it,  as  shown  in  Figure  1-5,  at 
A'B'C'D'  . Contours  of  constant  tn|L/(l+L) | are  available  on  the 
Nichols'  chart.  Using  these  contours,  it  is  seen  that  the  maximum 
change  in  Jln|L/(l+L)  | which  from  (2),  is  the  maximum  change  in 
tn  | T | is,  in  this  case,  very  closely  (-0.49)  - (-5.7)  = 5.2  db  , 
the  maximum  being  at  point  C'  , the  minimum  at  point  A'  . 

Suppose  that  the  specifications  tolerate  a change  of  6.5  db  at 
w * 2 , so  the  above  trial  position  of  | L ( j 2 ) | is  in  this  case 
more  than  satisfactory.  The  template  is  lowered  on  the  Nichols' 
chart  to  A"B"C"D",  where  the  extreme  value  of  tn|L/(l+L)|  are 
at  C"  (-0.7  db)  , A" (-7.2  db)  . Thus,  if  Arg  Lft(j2)  = -60°  , 
then  -4.2  db  is  the  smallest  magnitude  of  L^(j2)  which  satis- 


L. 


i 


fies  the  6.5  db  specification  for  A An  |t|  . Any  larger  mag- 
nitude is  satisfactory  but  represents  over-design  at  that  frequency. 
The  manipulation  of  the  oj  = 2 template  is  repeated  along  a new 


vertical  line,  and  a corresponding  new  minimum  of  | ( j 2 ) | found. 

Sufficient  points  are  obtained  in  this  manner  to  permit  drawing  a 

continuous  curve  of  the  bound  on  Lft(j2)  ' as  shown  in  Figure  1-5. 

The  above  is  repeated  at  other  frequencies,  resulting  in  a family 

of  boundaries  on  L. ( j 2 ) . 

A 


1.2.2.  Nature  of  the  bounds  on  L(ju>)  . 

A typical  set  of  bounds  is  shown  in  Figure  1-6.  The  bounds 


-360°  -240®  -120°  0° 

DEGREE 

Figure  1-6.  Typical  bounds  of  L in  Nichols'  chart. 
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tend  to  move  down  in  the  Nichols'  chart  (become  less  onerous), 
obviously  because  as  u>  increases,  greater  change  in  | T ( jo>)  | 
is  permitted,  as  in  Figure  1-2 (b) . It  is  in  fact  essential  that 
at  large  enough  u>  , the  uncertainty  in  | T ( ju>)  | (i.e.,  the  bounds 

on  |T( ju)  | ) be  greater  than  the  uncertainty  in  P(ja>)  , because 
the  net  sensitivity  reduction  is  always  zero  in  any  practical 

system  as  was  long  ago  (1]  shown  by  Bode, 

00  00 

/ *n  |Sp(jw)|du)  = - / fen  |l+L(jw)|  dw  = 0 (1.2-6) 

o o 


where  is  the  sensitivity  function. 

P O P/P 


In  the  above  example  as  u -*•  00  , P -*•  — j , so 


ka 

s' 


A fen  |p|  -*•  A fen(ka)  = 40  db  . Note  in  Figure  1-2  (b)  that  the 
permitted  A fen  | T( ju) | » 40  db  for  u > 50  . Such  large  toler- 
ances on  | T ( jcu)  | at  large  w are  tolerable  because  |T(jto)|  is 
negligible  at  large  w , e.g.,  if  |p(ju>)|  can  change  at  most 
by  40  db  at  large  u>  but  | T ( joj)  | changes  by  52  db  , who  cares 


if  this  52  db  change  is  from  | T| min  = 10 


to  T to 

' 1 max 


400  x 10  .In  return,  one  can  concentrate  the  sensitivity 
reduction  over  the  bandwidth  of  T(jw)  . Thus,  although  |p(jw)| 
in  this  region  varies  by  say  40  db  , | T ( ju>)  | may  be  controlled 
to  vary  by  only  4 db  , or  0.04  db  if  desired. 


1.2.3.  Universal  high-frequency  boundary. 

As  noted,  in  the  high-frequency  range  A fen  | T ( jw) | must 
realistically  be  allowed  to  be  » A fen  |P(jw) | , and  this  is 
reflected  in  the  bounds  on  Lo(ju>)  tending  to  a very  narrow  pencil. 


■ip  m mh 


In  Figure  1-7,  is  drawn  for  the  case  A fcn  L * A An  k = 20  db  , 

A in  | T ( ju>)  | = A in  |l/(1+L) | 5 23  db  at  wv  • However,  the 

resulting  peak  value  of  1l/(1+L)|  is  23  db  >=  14.1  arithmetic  at 

k = k x , indicating  a highly  under-damped  pole  pair  at  the 

corresponding  frequency  with  damping  ratio  £ = 0.034  , when 

k = kmax  • This  tremendous  peaking  does  not  appear  in  the  system 

response  to  the  command  inputs  R , because  it  is  filtered  out  by 

the  pre-filter  F in  Figure  l-3(a) . But  the  system  response  to  a 

C -1 

disturbance  D is  given  by  **  — » (1+L)  . Disturbance 

attenuation  generates  its  own  requirements  on  L , which  may  lead 

to  more  stringent  bounds  on  L than  those  due  to  T(ju>)  . The 


L L (DEGREE) 

Figure  1-7  Bounds  on  L(jw)  on  Nichols'  chart 
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final  contours  used  in  the  design  [2]  must  be  the  most  stringent 
composite  of  the  two.  However,  even  if  D is  very  small,  it  is 
usually  certain  that  a peak  |Tj|  of  14.1  is  intolerable.  It  is 
reasonable  to  add  a requirement  |t^|  $ y some  constant,  for  all 
u)  and  over  the  whole  range  of  P parameter  values.  The  resulting 
constraining  contours  denoted  by  are  shown  in  Figure  1-7 

for  the  case  A Jin  k = 20  db  , and  for  y = 2.3,  3.5,  5 db  (all 
these  contours  are  symmetrical  with  respect  to  the  vertical  line 
Arg  L = -180°  on  the  Nichols'  chart) . If  y = 5 db  is  used,  then 
B(wv)  indicates  the  composite  contour  shown  in  Figure  1-7.  For 

w > u>v  , | A T(jui)  | increase  while  y remains  the  same,  so  that 

sooner  or  later  there  is  reached  a frequency  w 9 B (gj)  is  a fixed 


boundary  , effective  V uj  > 


u 


1.2.4.  The  optimum  L(jh>)  . 

It  has  been  shown  [3]  a realistic  definition  of  optimum  in 

the  fti  system  is  the  minimization  of  k , defined  by 

lim  L(s)  = k s , where  e is  the  excess  of  poles  over  zeros 
S-x» 

assigned  to  L(s)  . 

It  has  been  proven  [3]  that  the  optimum  L lies  on  its 
boundary  B^  at  each  ah  and  that  such  an  optimum  exists  and  is 
unique.  Most  important  for  the  present  purpose,  is  that  in  signi- 
ficant plant  ignorance  problems  the  ideal  optimal  L has  the 
properties  shown  in  Figure  1-8,  i.e.,  over  a significant  range  it 
follows  B^  along  UV  up  to  the  point  J at  which  it  abruptly 


jumps  to  infinity  along  WW'W"  and  returns  on  the  vertical  line 


i 

N 
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YZ  , whose  phase  is  (-90°) *e  . Such  an  ideal  L(jw)  is,  of  course, 
impractical.  A practical  suboptimum  L is  shown  in  Figure  1-8. 


L L (degree) 


Figure  1-8.  Bounds  on  L and  Optimum  L on  Nichols'  Chart. 


Some  results  of  a numerical  design  example  are  shown  in 
Figure  1-9.  They  were  derived  for  the  following  problem. 

1.2.5.  Numerical  example.  (Figure  l-3(a)). 

Plant:  P = P^ 

Plant  ignorance:  P.  - k,/s  , /2  $ k,  $ 10  /2 

11  1 

P2  - k2/s  , SI  $ k2  * 10  /2 

Performance  Specification:  Shown  in  Figure  1-2 (b)  were 

originally  derived  from  time 
domain  bounds  of  Figure  1-2 (a) . 


Disturbance  response: 


Y S 2.0  db  . 


■J 


1 


16 


Figure  1-9.  Single  loop  Lg  and  bounds  of  a numerical  example. 


This  example  is  used  later  in  Chapter  2 as  a vehicle  for  presenting 
the  plant-modification  design  technique. 


1.3  Cost  of  Feedback  and  Effect  of  Sensor  Noise 
In  significant  plant  ignorance  problems,  there  is  a strong 


tendency  for  the  design  to  be  such  that  N , in  Figure  1-3 (a),  is 
so  highly  amplified  as  to  saturate  the  plant  input  at  X . The 
noise  response  function  is  (see  Figure  l-3(a)) 


« -L/P  in  h.f.  range.  (1.3-1) 

The  noise  response  of  the  numerical  design  example  of  the  last 
section  is  shown  in  Figure  1-10.  Notice  that  the  noise  component 
at  x , in  Figure  1-3 (a) , is  most  important  in  the  high-frequency 
range  where  the  useful  command  and  disturbance  components  due  to 
D , are  relatively  small,  rather  than  in  the  low  frequency  range 
where  the  latter  are  relatively  large.  This  is  further  enhanced 
by  the  fact  that  arithmetic  scales,  shown  in  Figure  1-11,  must  be 

• used  to  find 
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Hence,  it  is  desirable  to  decrease  |L|  vs  u>  , as  fast  as 
possible  in  the  high  frequency  range.  Even  a saving  which  is  small 
in  the  logarithmic  scale  near  A in  Figure  1-10,  can  be  signi- 
ficant in  rms  sensor  noise  effect. 


U) 


\ 

Figure  1-11.  Arithmetic  plot  of 

1.4  Reduction  in  Cost  of  Feedback 

1.4.1.  Linear  time-varying  compensation  and  nonlinear 
compensation . 

To  reduce  the  hf  sensor  noise  effect,  one  way  is  by  linear 
time-varying  compensation  if  the  problem  has  time-varying  features 
[7] . Another  is  by  nonlinear  compensation.  Actually  the  so-called 
"adaptive’  system  is  in  the  category  of  nonlinear  compensation. 

They  may  or  may  not  be  better  than  tti  compensation  in  reducing 
the  'cost  of  feedback'.  It  is  noteworthy  and  scandalous  that  in 

L ...  - A 
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the  vast  literature  on  adaptive  systems,  there  is  hardly  ever  any 
quantitative  comparison  between  the  adaptive  design  promoted  and  a 
proper  £ ti  design  accomplishing  the  same  design  objectives.  One 
could  excuse  this  not  being  done  in  a general  manner,  because 
there  is  hardly  any  'adaptive'  method  permitting  quantitative 
design  in  the  sense  here  defined.  However,  it  could  at  least  be 
done  experimentally.  Occasionally  one  sees  a comparison,  with  an 

I 

'ordinary'  or  so-called  'classical'  design.  But  the  comparison 
is  usually  greatly  biased,  because  generally  some  very  naive  £ti 
design  is  used,  and  there  is  no  statement  of  specifications  - 

I 

even  made  up  after  the  fact.  There  is  not  recalled  a single  com- 
parison, on  the  part  of  the  proponents  of  adaptive  systems,  with 

I 

the  £ti  quantiative  design  technique  [2]  discussed  here.  Some 
nonlinear  compensation  techniques  for  which  a quantitative  design 
theory  exists  to  a greater  or  lesser  extent  have  appeared  in  the 
literature  [8,9,10,11]  for  which  such  comparisons  are  possible. 

It  is  notevrarthy  that  these  were  expressly  motivated  by  the  desire 
to  reduce  the  'cost  of  feedback' , so  that  such  comparisons  were  a 
natural  by-product. 

1.4.2.  Multiple- loop  feedback. 

Another  method  of  'cost  of  feedback'  reduction,  in  the  con- 
text of  ilti  design,  is  by  means  of  multiple-loop  feedback, 
restricted  to  those  cases  where  additional  plant  variables  (besides 
the  plant  output)  are  available  for  feedback  purposes.  Such  a 
multiple  loop  design  technique  was  first  developed  [5]  for  the  cas- 


caded structure  of  Figure  1-12. 
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for  the  above  non  P.M.  multiple  loop  designs.  Hence,  it  is  very 
important  to  first  thoroughly  understand  the  essentials  of  the 
above  design  philosophy,  which  is  therefore  next  presented,  since 
this  present  work  is  confined  to  the  cascaded  plant,  only  the 
latter  is  presented.  It  suffices  also  to  consider  the  cascaded 
two-section  case. 

1.4.3.  Cascaded  2-loop  design  — no  P.M. 

The  basic  idea  is  to  use  the  inner  loop  L2  = G2P2  in 
Figure  1-15 (a),  to  minimize  the  effect  of  sensor  noise  at  the 


Figure  1-15 (a).  Cascaded  2-loop  system  — no  P.M. 


plant  input  • This  effect  is 
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1 + P.  G P 

1 1 2e 


where  P_ 
2e 


G2P2 

1+G2P2 


V(P1P2) 

1 + L, 


PlGlP2e 


(1.4-3) 


(1.4-5) 


w - - ■ - in  the  hf  range  where  | L | « 1 . (1.4-6) 

P1P2  1 

2C 

Hence  to  reduce  — jj-  in  the  hf  range,  one  must  try  to  reduce  L . 

1 

But  L must  cope  with  the  uncertainty  in  P p = P. (P?C  /1+P  G ) ) . 
Assuming  the  worst  case  of  uncorrelated  uncertainties  in  P^  and 
P2  , the  best  that  can  be  done  by  the  inner  loop  P2G2 
wipe  out  the  uncertainty  of  P2e  , so  that  need  only  cope 

with  the  uncertainty  on  P^  . Physically,  this  makes  sense  - for 
obviously  the  inner  loop  cannot  take  care  of  the  uncertainty  in 
P^  . There  is  then  left  the  single-loop  system  of  Figure  1-15 (b) 
and  Lj  can  be  designed  to  handle  the  ignorance  of  P^  only.  The 
resulting  is  therefore  more  economical  in  bandwidth  than  its 


Figure  1-15 (b).  Equivalent  single-loop  structure 
of  a cascaded  2-loop  system. 
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counterpart  in  the  single  loop  system,  for  the  realization  of  the 
same  C/R  specifications.  For  example,  compare  the  appropriate 
Tjj  in  Figure  1-16 (a)  (logarithmic  scale)  and  in  Figure  1-16 (b) 
(arithmetic  scale)  of  the  numerical  example  in  Section  1.2-5. 


UN,!  single  loop 


20  h |TN,|  cascaded 


Figure  1-16 (a).  Comparison  of  noise  response  |tn^|  on  Bode  diagram. 


This  looks  very  good,  but  the  obvious  question  is:  What  of 
the  effect  of  sensor  noise  N2  not  present  in  the  single  loop 
design?  It  would  first  appear  that  the  inner  loop  G2P2  would 
have  to  be  enormous  in  magnitude  and  bandwidth,  leading  to 
tremendous  effect  of  N2  . But  this  is  not  so.  And  the  basic 
reason  is  available  if  one  studies  the  mechanics  of  sensitivity 
reduction  by  frequency  response  methods.  The  reason  is  that 
design  of  the  outer  loop  to  handle  a certain  definite  amount  of 
uncertainty,  even  though  designed  optimally,  is  nevertheless  able 


Figure  1-16 (b).  Comparison  of  noise  response 

| Tn  | on  arithmetical  scale  plot. 

to  handle  a "much  larger"  amount  of  uncertainty.  This  is  nicely 
seen  in  the  Nichols'  chart  in  Figure  1-17. 

Thus,  in  Figure  1-17,  suppose  the  uncertainty  in  P^jw)  is 
given  by  the  template  shown  of  which  is  not  a point  (it  would 

be  a point  if  there  was  no  uncertainty)  but  a region.  Suppose  the 
specifications  require  the  closed  loop  response  uncertainty  to  be 


various  frequency  ranges. 

Thus,  after  has  been  designed,  one  finds  what  ignorance 

of  P2e  can  actually  be  tolerated  in  the  abrve  design  of  1,^  . H 


is  found  that  in  only  one  frequency  range  it  is  important  to  com- 


promise  (overdesign)  the  outer  loop  , in  order  to  ease  the  L2 
design  problem.  This  is  in  the  range  IVK  in  Figure  1-8,  where 

trade-off  between  the  two  should  be  made. 

l2 

Since  P2e  *«  ~L~  in  Figure  1-15 (a)  , with  L2  * G2P2  ’ the 

final  step  is  to  determine  the  L2  needed  so  that  the  resulting 

P2e  does  indeed  stay  within  the  bounds  found  by  the  previous  step. 

This  is  precisely  a single  loop  design  problem  with  P„  and  its 

2e 

tolerances  replacing  the  T function.  It  was  demonstrated  (51 
that  the  resulting  two- loop  design  could  be  highly  superior  to  a 
single  loop  design,  in  the  sense  of  achieving  the  seme  quantitative 
sensitivity  specifications,  but  with  considerably  less  effect  of 
sensor  noise. 

The  numerical  single  loop  design  example  of  Section  1.2  was 
also  done  by  a 2-loop  cascaded  design.  The  outer  loop  with 

its  bounds  is  shown  in  Figure  1-18  and  the  inner  loop  with  its 


bounds  in  Figure  1-19. 
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Step  2 . Plot  |Pn|  = |P^nP2n^  ' T^e  difference  between 

I L I and  I P | in  the  region  where  I L I < 1 , is  the  N, 

1 sn1  1 n1  3 1 sn 1 1 

sensor  noise  amplification  (for  its  effect  at  the  plant  input). 


Step  3.  Note  the  hf  uncertainty  of  P2  . This  is  the  maxi- 
mum amount  that  can  be  saved  in  the  hf  range,  giving  | L | in 
Figure  1-20.  It  is  always  worth  having  a few  db  overdesign  of 

, giving  I^n  in  Figure  1-20.  Note  that  Lan  is  obtained  by 
simply  introducing  the  high  frequency  characteristic  XYZ  of  Lgn  , 
at  point  X which  is  20  db  higher  than  point  X . For  example, 

cl 


in  Figure  1-20,  P_  = k_/(s  ) at  hf  with  k. 


../k- 


= 20  db  , 


then  lLan|  is  obtained  by  inserting  the  XYZ  pattern  at  x&  . 

For  5 db  gain  margin,  point  X^  is  5 db  lower  than  Xa  and 
| | is  the  curve  obtained  by  inserting  XYZ  at  Xfa  . 

Step  4.  L2n  can  already  be  approximately  drawn,  as  follows. 
|l^nl  is  used  for  lL^nl  • Then  max  |l2iJ  is  near  u>b  in 

Figure  1-20,  i.e.,  near  the  middle  of  the  fairly  flat  region 
preceding  the  large  high-frequency  slope.  There  is  hardly  any 
obligation  on  L2  in  to  < tob  . The  shape  of  lL2rJ  for  w > wb 


is  fairly  standard.  Its  slope  is  ~30  db/decade  from 


% 


to  to. 


bl 


in  Figure  1-20  until  |L,  ( ju>)  I = -[k_ 


/k. 


+ Gain  Margin] 


For  to  > to.  _ , the  shape  of  L_  depends  on  the  excess  of  poles 
b2  2n' 


over  zeros  of  L2  which  is  due  to  G2P2 


Because  the  excess  of 


poles  over  zeros  of  is  due  to  = GiG2P2Pl  , the  magni- 

tude of  the  negative  slope  of  |l2|  is  less  than  that  of  | | 


Step  5.  Sketch  | P 


2n' 


The  difference  between 


J2n 1 


and 


( 

i 


i 
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?2n  is  N2  sens°r  noise  amplification  because  in  Figure  1-15 


(a) 


N2'n 


A (JL, 

= VKT  • 


N2  n P+G2P2n+GlG2PlnP 


2n 


(1.5-1) 


• y<i«y2„> 


P+GlG2Plnp2n^ ^ p+G2P2n^ 


(1.5-2) 


- L2n/P2n 


(1+Lln) (1+L2n)  ' L2n  = G2P2n  ' 


L,_  i?  G, P,  P 


In  = Tin  2en 


(1.5-3) 


• 2n 


'2en  “ 1+L 


2n 


2n  . 


w “ p in  the  hf  range  where  |l,  I *■  1 
*»-  In' 


2n 


Il2„!  * 1 • 


(1.5-4) 


6-  0ne  can  now  decide  whether  to  use  , or  to 

compromise  further.  Suppose  that  more  compromise  is  wanted  - 

e.g.  x more  db  overdesign  giving  L,  for  outer  loop.  Then  the 

inner  loop  |l2J  maximum  level  can  be  reduced  by  x db  relative  to 

I I*2b  I max  * as  shown  in  Figure  1-21. 

It  is  important  to  note  that  the  above  is  based  on 

P1,P2  ^ ki/s  in  the  “-range  at  which  L.  is  at  B of 

1 h 

Figure  1-8.  This  is  reasonable  in  large  hf  uncertainty  problems, 
because  these  require  slow  decrease  of  |l|  over  large  w range. 

It  is  in  these  problems  that  the  complexity  of  multiple  loop  design 
is  warranted. 


J 


1.6  Plant  Modification  System  Synthesis 
1.6.1.  Plant  modification  (P.M.)  structure. 

It  is  again  emphasized  that  all  the  multiple-loop  systems 
discussed  above  [5,12,13]  are  restricted  to  non-P.M.  structures. 

In  this  case,  the  degree-of-freedom  available  for  the  design  is 
limited.  For  a n cascaded  plant  system,  ( n plant  variables 
available  for  measurement) , there  are  n independent  feedback  ioops 
which  may  be  used  — if  no-P.M.  is  allowed,  i.e.,  if  the  output  of 
all  the  feedback  processors  are  allowed  only  to  the  plant  input. 
Together  with  the  prefilter,  this  gives  a (n+1)  degree  of  freedom 


system,  as  shown  in  Figure  1-11.  If  feedback  to  internal  plant 

variables  is  allowed,  then  the  number  of  independent  feedback  loops 

can  be  greatly  increased.  In  a cascaded  plant  with  n-1  internal 

variables  and  one  output  variable  the  number  is 
n 

£ (n-i)  » (n+1) *n/2  . Now  the  total  numbers  of  freedom  available 
o 

is  (n+l)n/2  + 1 , including  the  prefilter  as  shown  in  Figure  1-22. 


Figure  1-22.  The  P.M.  multiple-loop  with 

n(n+l)/2+l  D.O.F.  structure. 

1.6.2.  The  plant  modification  problem. 

Consider  a plant  consisting  of  two  cascaded  sections,  as  in 
Figure  1-23 (a)  and  suppose  a certain  maximum  output  signal  level  is 
required,  with  Laplace  transform  Cjls)  . The  signal  level  at  the 
input  of  PL  is  then  X ^ - Cj/Pj^  * c2  (out  of  p2  ) * This 
remains  true  in  the  cascaded  feedback  structure  (no-P.M.).  Thus, 
the  feedback  designer  does  not  affect  the  signal  level  in  the  plant. 


r:TA. 


needed  to  obtain  a specific  output  level.  This  is  true  for  any 
feedback  structure  in  which  all  the  feedback  paths  return  to  the 
plant  input.  If  a feedback  loop  is  put  around  plant  , by 

means  of  in  Figure  1-23 (b) , then  the  signal  at  is  still 


X2  p2  c2  I X,  P Cl 

A.  • * /N 


c,(t) 


p2 S?  i x‘  Pi  5l 


Figure  1-23  (a)  A simple  two  plant  system. 

(b)  Adding  a P.M.  loop  = PiHi 


Cj/P  * But 


C2  = Xi+Hici  “ Cipi+H1ci 


(1+H^) = (l+H^J-Xj. 


(1.6-1) 


(1.6-2) 


For  the  same  signal  level  of  and  therefore  of  X.^  , the 

signal  level  at  C2  is  now  multiplied  by  (1+P^)  , which  con- 
stitutes a modification  of  the  plant.  It  is  concievable  that  this 
new  level  of  C 2 may  be  so  much  larger  than  the  old  one,  that  P2 
may  have  to  be  rebuilt  to  be  able  to  handle  this  larger  signal 


level. 
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If  the  'feedback  expert’  is  working  together  with  the  'plant 
expert'  in  the  design  of  the  plant  itself,  then  the  trade-offs  in 
such  significant  plant  modification  may  be  seriously  considered. 

But  in  many  cases  this  is  not  so.  It  is  advisable  to  find  a design 
method  which  uses  the  P.M.  loop  to  improve  the  sensor  noise  problem, 
but  keeps  the  increase  in  the  signal  level  in  the  plant  within  a 
tolerable  range.  Recall  that  the  noise  problem  is  significant  in 
the  hf  range  (high  relative  the  bandwidth  of  the  useful  control 
signals) . This  gives  us  a very  good  opportunity  to  achieve  signi- 
ficant improvement  in  sensor  noise  effects  with  only  small  or 
moderate  signal  level  increase.  In  addition,  the  insight  obtained 
from  such  a synthesis  procedure  is  very  useful  for  those  cases 
where  the  feedback  expert  is  called  in  to  help  the  plant  expert  in 
the  actual  plant  design  stage. 


1.6.3.  Reduction  in  cost  of  feedback  by  P.M.  structure. 

The  cost  of  feedback  can  be  tremendously  reduced  by  a P.M. 
system.  This  is  illustrated  by  the  simple  2-plant,  3-loop  P.M. 
system  shown  in  Figure  1-24.  The  noise  effect  at  X1  and  x2  due 
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Figure  1-24.  A P.M.  2-plant,  3-loop  system  structure  with 

unit  feedforward  element  between  plant  section. 
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(1.6-5) 


(1.6-6) 


(1.6-7) 


« -L.  /P_  in  hf  range  where 
X2 


lL0l  « 1 * I4  I « 1 


(1.6-8) 


1 A X1  "P2H2 
Tf,2  " "T  = ^i^  (l+I^J+GP^ 


(1.6-9) 


- L / [ (1+L . ) (1+L.  )) 
2 1 12 


^ - L^  in  hf  range  where  |Lq|  « 1 , 

|L  | « 1 , |l . I « 1 
1 2 


(1.6-10) 


Following  the  same  design  philosophy  as  in  a non-P.M.  cas- 
caded system  design,  let  L.  cope  with  the  uncertainty  in  p . 

*■2  1 
But,  for  the  sake  of  the  significant  signal  level  variation  due  to 


L.  , it  is  impractical  to  let  L.  cope  conpletely  with  the 
1 1 

uncertainty  in  P.  . Let  L.  cope  with  P uncertainty  in  the 

1 1 

high  frequency  range  where  sensor  noise  is  significant,  while  the 

uncertainty  in  the  low  frequency  range,  where  the  control  signals 

dominate,  is  taken  care  of  by  the  outer  loop  L . Then,  L can 

o o 

be  designed  as  a highly  economical  loop,  in  terms  of  bandwidth. 

So  the  sensor  noise  effect  at  x2  due  to  in  (4)  becomes  very 

small  compared  to  the  single  loop  design. 

The  results  of  a numerical  example,  taken  from  Chapter  4, 
are  shown  here  to  illustrate  the  huge  improvement  in  sensor  noise 
effect,  i.e.,  excellent  saving  in  the  cost  of  feedback. 

Fioure  1-25 (a)  snows  the  Bode  plots  of  noise  response  and  the 
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Figure  1-25 (a) . Sensor  noise  effect  at  and  in 

single-loop  system  and  P.M.  3-loop  system. 

arithmetical  plots  are  in  Figure  1-25 (b) . The  effects  of  intro- 

2 1 

ducing  a new  sensor  noise  source  N_  , i.e.,  T and  T in 

N2  N2 

(8) , (10)  and  the  effect  at  , input  to  P^  , by  introducing 

a new  feedback  loop  H , i.e.,  T^_  in  (6)  are  all  shown  in 

1 Nx 

Figures  l-25(a),  (b) . It  is  known  that,  by  a moderate  trade-off 

between  inner  loops  L.  , L.  and  outer  loop  L , the  effects 

X1  12  ° 

2 1 1 

of  T , T and  T are  reasonable  small,  although  they  are 
2 ™2  1 

bigger  than  their  counter  parts  in  a single-loop  system.  Note, 
is  bigger  here  than  the  single  loop,  even  though  Lq  band- 
width « Lg  bandwidth  in  single  loop  system.  The  reason  is  that, 
in  the  3- loop  P.M.  system,  there  is  an  extra  path  — from 
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Figure  1-2 5 (b) . Sensor  noise  effect  at  X^  and  x2  in  single- 
loop system  and  P.M.  3-loop  system  — 
arithmetical  scale. 


to  . This  is  also  seen  in  (6),  as  the  |l^  | bandwidth  is  not 
small  compared  with  that  of  |Lg|  , even  though  the  |lq|  band- 
width is  relatively  small. 

It  appears  that  the  P.M.  structure  is  very  good  in  decreasing 
sensor  noise  effects.  However,  the  signal  levels  in  the  plant  to 
achieve  a specific  output  are  affected  by  the  introduction  of  the 
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P.M.  feedback  loop,  so  that  no  longer  does  one  really  have  a fixed, 
constrained  plant.  This  is  the  plant  modification  problem. 


1.6.4.  P.M.  structure  with  unit  feedforward  elements 
between  plant  sections. 

Note,  in  Figure  1-24,  we  put  a unit  feedforward  element 
between  P^  and  P2  . This  makes  it  somewhat  simpler  to  develop 
a systematic  synthesis  procedure.  Now  we  prove  this  system  is 
equivalent  to  the  system  without  separation  between  plant  sections, 
as  in  Figure  1-26.  The  transfer  function  in  Figure  1-24  is 


fgp1p2 


1+P1H1+P2H2+P1P2H1H2+GP1P2 


The  transfer  function  in  Figure  1-26  is 


f*g*p1p2 

1+P1H1+P2H5+P1P2G* 


(1.6-11) 


(1.6-12) 


Figure  1-26.  A P.M.  2-plant,  3- loop  system  structure 
— without  separation  between  plants. 


Compare  (11)  and  (12).  Let 


H1  = H1 


“5  - h2 


G = G+H1H2 


F*  = FG/tG+HjHj) 


(1.6-13) 


(1.6-14) 


(1.6-15) 


(1.6-16) 


and  then  the  two  systems  are  equivalent  with  respect  to  T and  the 

effects  of  , P2  uncertainty  on  T . Note  also  that  the  actual 

C, 

output  of  P2  is  not  xx  = C1/P1  but  is  X^l+PjH*)  * ^(1+PjH*) 

cl  1 

which  is  the  same  as  that  of  X^l+p^i^)  - — (1+p.^)  in 
Figure  1-24.  And  the  noise  responses  in  Figure  1-26  are 


2 A 2 

T..  £ — - 


-G*+H*H* 


Nx  Nx  1+g*p1p2+p1hJ+p2h* 


i Xi 

n^*  A 


- H*  - G*P2 


Nx  NJ_  i+g*p1p2+p1h’+p2h2 


t2  a x2 

N2  = N2  * 1+G*P1P2+P1H^+P2H2 


1 

T1  £ — . 


- H*P„ 


N2  N2  1+G*P1P2+P1H^+P2H^ 


Substitute  (13) , . . . , (16)  into  (17) , . . . , (20) 

2 -(G+HjH^+H^ 

\ " 1+ (G+H1H2 ) P1P2+P1H1+P2H2 


1+P1H1+P2H2+GP1P2+P1P2H1H2 


_x  -Hj^-  (G+H1H2)P2 

rNi  " i+(g+h1h2)p1p2+p1h1+p2h2 


(1.6-17) 


(1.6-18) 


(1.6-19) 


(1.6-20) 


(1.6-21) 


r 

. 

u 

p 


11 

* 


- [Hj^d+PjHjJ+GPj] 

(1+1*  ) (1+L.  )+P.P~G 

il  2 
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(1.6-22) 


2 -^(l+P^) 

rN2  " 1+  (G+H^)  P1P2+P1H1+P2H2 


Hjd+P^) 


(1+L.  ) (1+L.  )+P.P,G 
1 2 


(1.6-23) 


-h2p2 


N2  1+(G+H1H2)P1P2+P1H1+P2H2 


-P2H2 


(1+L  ) (1+L.  ) +GP. P_ 

X1  2 


(1.6-24) 


(21) , . . . , (23)  are  exactly  the  same  as  in  (3),  (5),  (7),  (9)  of 
Figure  1-24. 


For  simplicity,  henceforth,  we  always  choose  the  structure 
with  unit  element  between  plants,  for  development  of  the  P.M. 
synthesis  theory. 


1.6.5.  The  RMS  signal  level  problem. 

The  P.M.  system  synthesis  theory  in  the  present  work  is 
restricted  to  the  RMS  signal  level  problem  — not  to  peak  values. 
That  is,  this  synthesis  theory  is  based  on  the  amount  of  signal 
level  variation  Q defined  by 
00 

/ |c4  (oj)  1 2 du> 

* i,m  'max 

Qi  * • 

/ |C.  (o))|2  do) 

JQ  1 i , 8 ' max 

where  |C.  (u) I lathe  maximum  signal  level  of  P.M.  multiple-loop 

i , m max 

system  at  stage  n-  i such  that  there  is  at  least  one  P.M.  inner 


i - 2,...,m  (1.6-25) 

for  m-plant  system 
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loop  in  stage  n-i+1  dr  in  later  stages;  |c  (w)  | Is  the  maximum 

1 i s max 

S.L.  of  the  single-loop  system  at  the  corresponding  n-i  stage. 

1.6.6.  The  P.M.  synthesis  philosophy. 

P.M.  design  based  on  controlling  the  RMS  value  of  the  effect 
of  plant  modification,  has  the  nice  property  that  the  individual 
contributions  at  each  value  of  u in  [0,®)  add.  Now  the 
sensitivity  reduction  properties  in  the  multiple-loop  feedback 
system  — the  trade-offs  between  the  loops,  etc.  — vary  considerably 
in  different  parts  of  the  frequency  spectrum.  This  is  evident  from 
Sections  1.4,  1.5.  This  permits  the  designer  to  divide  up  the 
spectrum  into  significant  portions  and  in  each  portion,  concentrate 
on  the  crucial  properties.  Consider  how  complex  the  problem  would 
be  otherwise  — quantitative  design  for  significant  parameter 
uncertainty,  multiple  loops,  sensor  noise  effects  and  plant  modi- 
fication. Using  frequency  response,  one  can  see  the  forest  from 
the  trees  and  obtain  good  engineering  solutions.  Modern  control 
theory  with  its  contempt  for  frequency  response  and  non-analytic 
solutions,  has  not  been  able  to  cope  quantitatively  with  the 
uncertainty  problem,  even  in  single-loop  design,  despite  its  vast 
number  of  competent  researchers. 

It  is  assumed  that  the  design  specifications  dictate  the 
permissible  Qi+1  in  RMS  signal  level  at  the  outputs  of  P2,...,P 

in  Figure  1-22,  which  is  allowed  in  the  P.M.  design.  The  design 
technique  exploits  this  permitted  increase  so  as  to  simultaneously 
decrease  the  net  effects  of  the  various  sensor  noise  sources. 
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in the  very  low  frequency  range,  the  control  signals  are 
generally  large,  so  even  a small  increase  percentage-wise  can  have 
a significant  effect  on  the  signal  level  change.  Fortunately,  the 
outer  loop  (as  is  evident  from  Sections  1.4,  1.5)  can  do  the 
sensitivity  reduction  job  in  this  range,  with  little  help.  There- 
fore, the  P.M.  loop  is  designed  in  this  range  to  be  as  small 

as  possible,  in  order  to  guarantee  that  the  relative  increase  in 
signal  level  (of  P2  , at  C2  in  Figure  1-24)  is  very  small.  In 
the  middle  frequency  range,  the  signal  level  is  allowed  to  increase 
by  amounts  related  to  . And  the  P.M.  loop  is  assigned  to 

participate  to  some  extent  in  the  uncertainty  problem.  In  the 
high  frequency  range,  the  signal  level  problem  is  relatively 
unimportant,  so  the  P.M.  loop  is  designed  to  minimize  the  cost  of 
feedback,  i.e.,  to  minimize  the  effect  of  sensor  noise  without 
worrying  over  the  S.L.  effect. 
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CHAPTER  2 


SIGNAL  LEVEL 


2.1  Introduction 

The  signal  level  (S.L.)  problem  is  the  basic  one  in  plant  modi- 
fication (P.M.)  system  design.  If  there  is  no  limitation  on  the 
output  of  each  plant  section,  then  P.M.  is  easy.  In  a practical 
system,  the  designer  must  understand  the  effect  of  the  local  feed- 
back loop  P_^H^  on  the  signal  level  of  the  preceding  stage. 

This  chapter  presents  a detailed  analysis  of  this  effect.  The 
signal  level  variation  ratio  (SLVR)  p is  defined.  Then  the 
problem  is  simplified  by  dividing  the  frequency  spectrum  into 
distinct  frequency  ranges.  A single  simple  factor  |l  + Linl » Is 
shown  to  dominate  p in  each  frequency  region.  A power  increase 
tolerance  level  is  assinged  to  any  P.M.  design  and  it  is  finally 

shown  how  to  relate  p and  |l  + L,  I to  this  tolerance  level. 

in 


2.2  Definitions  of  Signal  Level  Variation  Ratio  (SLVR)  and  p 
SLVR  is  defined  with  repsect  to  a single  loop  design  for  the 
same  plant.  Let  fi  **  {P}  be  the  plant  set  due  to  uncertainty. 
Single  loop  system  (see  Figure  2-1  for  notation) 

The  maximum  output  signal  level  is 


IV^Lax 

r 


i 


■ 


(2.2-1) 
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and  the  maximum  signal  level  at  the  output  of  P.,  is 


lC2s<J“>Lax 


|RFs(ju>) 


La(ja>) 


P,  (jo))*  (1  + L ( jio) ) 'max 
i s Jp 


(2.2-2) 


The  ratio 


Ys(Jw) 


|C2s(ju>)  | max 
[ C (ja)>  | max 

s r 


I Ls  ( ju)  / (P-L  ( ju>)  • [ 1 + Lg  ( ju)  ] ) I max 

[L_(»/[l  + L (jw)]  | max 
s s Jp 


(2.2-3) 


5*  is  dropped  henceforth,  i.e.,  maximization  and  minimization 
is  always  over  P unless  otherwise  specified. 

P.M.  2- loop  system  (see  Figure  2-2  for  notation) 

Similarly, 


|C.  (ju>)|  - 1 RF(joi)  1 • [ V I 

1 1 J 'max  1 J ' 1 1 + Lptju))  'max 


(2.2-4) 


C22«">L*  ' lRr<3“> 


L0(ju)) 


'Ple(jw)  [1  + Lq( jco) ] 'max 


(2.2-5) 


5 |c22(jM)Im<  |L0(jM)/<ple(wu^u(jM)i)|„, 

Y0(J<*>)  I r 1 4,.^  I I t rTTwTm  TTmTT!  ' 


~ C^jio) 


max 


lL0(»/[l  + L0(ja»)l,max 


such  a case  we  shall  use  at  each  u)  , that  parameter  combination 


Figure  2-3  Division  of  the  frequency  spectrum 

high-frequency  (denoted  by  hf)  uncertainty  of  be  db  (20  in 

our  example)  and  that  of  P2  be  Mj  db  (also  20  db  in  this 
example) . 

(1)  The  very  low  frequency  range  * [0,  is  that  in 

which  | Lg ( juj)  | 25  db  over  the  entire  plant  parameter 

it 

1 


space. 


In  Figure  2-3,  u)^  » 0.08,  because  Lg  there  is 

I Ls I zn±n  ' *n  ran8e  in  the  multiple  loop  design 

(cascaded  no  modification,  or  P.M.)  the  inner  loops  are 

not  used  to  help  the  outer  loop  Lq  , i.e.  in  R^, 

L ~ » where  L and  L designate  the  nominal 

sn  on  sn  on 

loop  transmissions. 

(2)  Middle  frequency  range  Rj  = (iu^ , u^]  defined  by 
25  db  > | Lgn ( jco)  I _>  -M2  db  where  M2  is  the  hf 
uncertainty  of  P2  , which  is  20  db  in  this  example. 

In  Figure  2-3,  (*>2*8. 

(3)  High  frequency  range  1,  R-j  *=  (u^,  defined  by 

-M2  db  > |Lsn(jio>|  > -(M1+M2)db  where  (M1+M2)  is  the 
hf  uncertainty  of  P = P^P2*  In  Figure  2-3,  w-j  = 44. 

(4)  High  frequency  range  2,  R^  = (co^,  where  00^  = 10  • 

In  Figure  2-3,  u)^  ■ 440. 

(5)  The  very  high  frequency  range  R^  * (u>^»  00 ) . 

It  is  helpful,  for  interpretation  of  SLVR,  to  point  out  here 
orders  of  magnitude  of  the  P.M.  inner  loop  (nominal  value  at 

|P|min)*  before  going  into  detail.  Two  numerical  examples  (from 
Chapter  3)  of  typical  P.M.  2- loop  designs  are  given  in  Figure  2-4. 

In  R^,  lLinl  very  small  in  general,  though  the  actural  values 

depend  on  the  specific  problem.  Thus,  in  Figure  2-4,  Design  A has 
jL^n(j(u)|  <_  -30  db  and  in  Design  B,  |L^n(ju))j  < -22  db  for 
oj  € Rx  . In  R2»  the  SLVR  is  allowed  to  increase  and  | I becomes 
bigger.  In  Figure  2-4,  -30  db  < |Lin(juj)|  < -7  db  in  Design  A 


Figure  2-4  Division  of  frequency  spectrum  and  corresponding 


L , L and  L.  of  design  examples, 
sn  on  in 


and  -20  db  < | (Jw)  | < -5  db  in  Design  B,  for  go  C R^  . In 
Ry  L is  used  to  cope  with  the  parameter  uncertainty  of  Pj  , so 


|Lln(ju))|  cannot  be  very  small  and  actually  tends  0 db  . In 


Figure  2-4,  -7  db  < |L^(ju>)|  £ -2.5  db  in  Design  A and 


-5  db  < |L^(ju))|  < -2  db  in  Design  B,  for  u>  « R^.  Note  that 


here  |L^n(ju>)|  is  much  bigger  than  the  corresponding  | L (Ju>) 


and  |Lgn(ja))[  - see  Figure  2-4.  In  R^,  in  which  LQn  design  was 
based  on  no  uncertainty  in  P^p*  is  nevertheless  satisfactory  for 
large  uncertainty  (cf  section  1.4.3),  so  Il^J  can  be  small. 
In  Figure  2-4,  -6  db  > |L^n(jw)|  > -50  db  in  Design  B and 

-4  db  > |L^(Ju>)|  > -44  db  in  Design  A,  for  0)  £ R^.  And  the 
last  range  R,.  is  the  very  high  frequency  range  where  both  Lgn  , 


Lon  * Lm  are  vory  small*  In  Figure  2-4,  |Lgn(ja.i)|  < -84  db  ; 


- 


obtained  by  changing  the  sign  of  constant  magnitude  curves  (see 
Figure  2-5).  These  loci  will  be  very  useful  in  P.M.  design. 

The  simplification  of  (2.2-7)  for  SLVRp,  follows. 


r-  -J» 
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-180*  *120*  -60°  0*  60*  I2cf  i80° 

£jLi  (DECREE) 

Figure  2-5  Loci  of  constant  |l  + Lj  on  the  Nichol's  chart 


2.4.1  Rx  - [0,  lx).] 

In  this  range,  the  single  loop,  unmodified  cascade  and  P.M. 

designs  are  almost  identical,  so  |l  I,  |l  | >25  db  (see 

Section  2.3  and  Figures  2-4,6).  Thus,  the  inner  loop  1,^  is  not 

used  in  to  help  the  outer  loop,  and  it  is  therefore  small  with 

typical  maximum  (over  ix)  ) values  of  the  nominal  L_,  (chosen  as 

in 

the  smallest  over  ?P|)  ^ ~30  db  (see  Figure  2-7  which  is  a 

typical  numerical  example  from  Chapter  3).  So  |Lq/(1  + Lo)|  = 1 
over  the  range  of  plant  uncertain,  and  in  (2.2-7) 


1 + L 1 

o near  max 


(2.4-4) 


I iLi  ( DEGREE! 


Figure  2-7 

where  'JP^ 
Similarly 


With  this, 


A typical  inner  loop  L - P^  (from  Chapter  3, 
Design  A) 


{Pj^}  , due  to  uncertainty  in  P^  parameters. 


/Ip. 


max 


1 + L • 1 * l 'min 

s near  max  ^ 


(2.2-7)  becomes 


(2 . i 
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1+L 


P - 


s max 


1+L 


1+L 


o near  max 


s near  max) 


'1+L 


o max 


1 max 


s max 


(2.4-6) 


£ 6, 


/ Ip  1 \ 

/ ' l^min  ' 

/ i 

V, 

/ c,  1 \ 

\ 1 

1 max  ] 

\ 'Ple^min  i 

1 | C | J 

\ 1 s'  max  J 

' rj 

(2.4-7) 


lPlLln  ^ 


le 1 min 

*l/ 


'1 'max 


(2.4-8) 


s 1 max , 


Consider  the  error  Introduced  In  letting  6^1  in  (7,8). 


For  example,  let  L have  its  minimum  value  at  25  db  ^■-90“ 
sn 


with  gain  uncertainty  of  P.  , 1 to  10.  |L  /(1  + L )|  - 

a s s 

| jl7.8/(j  + J17.8)|  - 0.99842  at  |P  | , and  | J178/(l  + J178)  | - 


0.99998  at  | | • So  the  percentage  error  in  cancelling  the 

x max 

*i 

first  factor  of  (6)  is  (0.99998-  0.99842)/0.99842  - 0.157%. 


Similarly,  the  second  factor  in  (6)  has  the  same  error  due  to 
cancellation.  But  note  that  in  the  cancellation,  the  second 
factor  in  (6)  has  an  error  opposite  to  that  of  the  first  factor. 
So  due  to  both  cancellations,  the  total  maximum  error  <<  0.157%. 


I 


So  in  (7),  6 *■  1 . Note  again  that  we  choose  the  nominal  plant 

to  be  at  the  minimum  value  over  V. 

Consider  |P.  I . / 1 P,  | . - |P.|  . /|P./(l  + Ld]  . in  (8). 

1 l'min  1 le  min  1 l'min  1 1 i 'min 


want  to  justify  (with  1^  = P^H^), 


1 jmin 

V, 


i 'min 
'Pi 


i min 


°1  i " 1 

.)  . I ~ H,  ll+  (L.)  . 
i min'  1 ' i mi 


(2.4-9) 


from  Figure  2-5  as  follows.  Recall  that  in  , in  the  P.M.  design 

philosophy,  the  inner  loop  is  not  used  to  help  the  outer  loop,  so 

IL  | ■ |l  I and  L.  is  rather  small.  For  example,  in  the 
1 sn1  1 on1  i 

design  in  Figure  2-7,  the  biggest  |Linl  over  u>  C R^  is 


-30  db  ^.61°  at  co  = 0.08.  This  values  is  |L 


in 1 min 


= {P^} . A typical  template  of  due  to  P^  uncertainty  is 
ABCD  in  Figure  2-8  in  which  the  loci  are  of  constant  1^/(1 + 1^)1 
It  is  obvious  that  both  ] L^/ (1  + L^) |min  and  J | min  are  at  A, 


•°  IVH+li>l^n  ‘ ‘ -20'73  db-  Eve”  1£ 


the  template  of  L^  is  AB’C’D  instead  of  ABCD  in  Fig.  2-8, 

IV(1+L1>U  ; lLiU/l1+(Li>„J  * lV(1+Li>»i„  ■ 

9°i  '•’i  n 

-20.8124  db  at  A,  IL. I . can  be  any  point  along  AB’ . Let  us 

l min 


take  B’ , the  furthest  from  A,  with  | | g, / | 1 + (L^g, I “ 

-20.S636  db.  So  62-|{|L1/(l+L1)loln}-{|Lilmii>/|1+tt1)min|}| 

f>l  Pj  V,  *»’ 

| 20. 8124  - 20. 5635 | » 0.2488  db  . Note  that  62  depends  on  the 


L Li  (DEGREE) 


Figure  2-8  Templates  of  in  Nichol's  chart  in 

position  of  template  . Thus,  if  A in  Fig.  2-8  is  at 
-30  db  /—  30°  instead  of  -20  ^.30°,  then  6^  =*  0.058  db  and  when 
A is  at  -10  db^_  30°,  62  “ 0.382  db.  Our  design  philosophy  is 
to  make  as  small  as  possible  in  , in  order  that  the 

minimum  signal  level  of  C22  be  very  closely  equal  to  that  of  C2 

This  is  because  the  control  signal  level  is  highest  in  this  low 

frequency  range,  so  that  even  small  percent  changes  could  lead  to 
large  absolute  differences.  So  62  is  small  in  R^  . Here  we 
took  /L-  ^ 0°  , because  R^  is  the  low  frequency  range.  Even 

if  the  system  is  type  1,  it  is  the  outer  loop  which  is  preferably 
made  type  1,  as  the  primary  function  of  the  inner  loop,  , is 

to  help  relieve  the  burden  on  the  outer  loop  in  the  hf  range.  But 

even  if  Z.  L.  in  R,  'V-  -90°  , the  error  60  , 


in  place  of 


TJ  . r 


0.2488  db  for  A at  -20  db  /L  30°  is  still  not  great,  e.g.  in 
Figure  2-8,  |L^/(1  + L^)j  ^ * -20.19  db  at  E and  the  furthest 

value  from  E according  to  |L1lm±n/|l+  (Li>mlnI  is  “19.892  db 


at  F'  and  m 0.298  db.  So  (8)  becomes 


p * a 


1 1 1 min 
yi 

PlLin^1+  *Li\niJ 


'1  'max 


(2.4-10) 


a • 1 + L, 


J1 'max 


(2.4-11) 


where  a - 1 in  R, 


Note  again  that  L^n  is  the  nominal  value  of  defined  as  the 

minimum  of  [ ( over  TP ^ . 

At  worst,  in  (11),  a very  slight  adjustment  can  be  made  in 
^l^max  respect  to  l^slmax  to  “‘ske  a - 1 exactly  or  even 

< 1 . Curves  of  ot  for  two  numerical  examples  from  Chapter  3, 
are  shown  in  Figure  2-9. 

The  conclusion  is  that  in  R^  , the  SLVR  p « |l+Lt  | • 

( I CjJ  max^ I Cs I max^  with  very  good  accuracy> 


2.4.2  R2  3 (u^,  u2]  (Figures  2-3, 4, 6, 7) 

Recall  R2  is  the  middle  u>  range  in  which  the  signal  level 
is  allowed  to  increase,  p > 1 , and  with  larger  Il^I  than  in 
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Figure  2-9  Nature  of  a - from  examples  in  Chapter  3 

From  section  2.3,  R,  ■ {to 1 25  db  > |l  (loo)  1 > -M„  db}  where 

2 sn  1 — 2 

M2  is  the  h.f.  uncertainty  of  P2  , i.e.  P2  - k2Ph2  , 

e2 

lim  P « 1/s  , e.,  » the  excess  of  poles  over  zeros  of  P-  , 

s-~*>  h2  2 2 

M2  - 2o  log  t't2max^2min^  * Note  c^at  the  M2  uncertainty  is 

handled  by  the  outer  loop  Lq  in  the  P.M.  2- loop  design.  If  P2 

has  significant  uncertainty,  R2  covers  a rather  big  range,  e.g. 
(0.08,  8.J  in  Figure  2-10(a),  which  is  replotted  from  Figure  2-6. 
It  is  important  to  note  that  even  at  a>2  (-8)  , at  the  end  of  R2  , 

the  template  of  L » G PjP2  at  lP]J  “ |PjJ  ^ fixed,  extends 

Pi 

above  the  horizontal  0 db  line,  in  the  Nichol's  chart  of 
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Figure  2-10(a).  Similarly,  in  P.M.  2- loop  system,  Lq  = GP^eP2 
at  |Plel  * lPie^min  fixed»  extends  above  the  horizontal  0 db 

Pi 

line,  in  Figure  2- 10(b),  which  is  replotted  from  Figure  2-6. 

Consider  the  first  two  factors  in  (2.2-7).  We  want  to  justify 


|L  /(1  + L J 
1 s s'  max 


rvtpi(1+VJ 


maxi 


|L  /[P  (1  + L )]| 

0 le o 'max 

| L / (1  + L)  | 

1 O o' ‘max 


1 1 min 


- p 


le  'min 


or: 


(2.4-12) 


'1  + L 


s max 


■ifr, 


min 


1 + L 


o max 


le 'min 


(2.4-13) 


'P.U  + L )» 

1 s max 


'P.  (1  + L )' 
le  o max 


This  is  done  by  illustrating  typical  templates  of  single-loop  and 

2-loop  P.M.  designs,  in  each  part  of  . 

It  is  concluded  from  Appendix  I that  there  are  three  kinds  of 

P^  , P2  patterns,  shown  in  Figure  2-11.  In  the  first  kind,  the 

template  extends  upward  and  to  the  right  from  the  nominal  point 

Pin  (see  section  2.2.1),  as  in  Figure  2-11 (a).  The  second  kind 

extends  upward  and  to  the  left  from  P , Figure  2-11 (b).  The 

in 

third  extends  both  to  the  right  and  to  the  left  from  Pin  , as 
in  Figure  2-11 (c). 


m.r 


! P.  I 
I " I 

Pin 


Figure  2-11  Three  patterns  of  template  of 


Consider  the  first  kind  of  pattern.  Typical  templates  of 

P.  , P_  at  to.  < o)D  < o)„  < U)_  are  shown  in  (a)  and  (b)  of 

Figures  2-12,...  15.  A typical  locus  of  Lin  = P^H^  in  the 

frequency  range  is  shown  in  Figure  2-16  with 

Lin(wA)  « -30  db /_  60°  , Lin(wB)  = -22  db  Z.  50°  , 

Lin(to^,)  * -10  dbZ_40°  and  L^n(uip)  = -4  db  Z_30°  . The  resulting 

templates  of  P^g  * P^/(l  + L^),  for  P^  in  (a)  of  Figures  2-12,...  15, 

are  shown  in  (c)  of  the  same  figures,  and  this  shape  of  is 

next  explained.  Recalling  section  2.4.1,  we  can  suppose 

P,  (defined  as  min  |P,  I = P.  /(1  + L.  ) with  error  6„  depending 
len  -p  • le'  In  in  2 

“l 

on  the  position  of  , but  very  small.  For  example,  if 

| I = -30  db,  6^  - 0.058  db  for  a 20  degree  wide  P^  template, 

6_  = 0.2488  db  at  |l,  | * -20  db  and  6.  = 0.382  db  if 
2 'in'  2 

| L^n | = -10  db  , for  the  same  P^  template.  In  (a)  and  (c)  of 
Figures  2-12,.  ..15,  the  nominal  P^n  » Pien  are  at  A'  ancl  A" 
respectively.  Note  that  the  P^  template  extends  upward  and  to 
the  right,  but  P^e  extends  both  to  the  left  and  right.  This  can 
be  explained  by  Figure  2-17  where  A,B,C,D,  correspond  to  A' , B' , 
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(l+U> 

X 


le'at  A 


Hence  D"  is  to  the  left  of  A" 


C' , D'  in  Figure  2-13(a).  Note  XO  = 1 , so  the  vector  from  X 

to  any  point  in  ABCD,  represent  1 + and  argL^  >0  in  R?  . 

Arg(l  + L^)  ^ at  A = -ZOXA  and  Arg(l  + L^)  ^ at  D = -ZoXD. 

Obviously,  Z.  OXD  > ZOXA  , soZPlelat  D - plelat  A > 

Z ”,  | - Zp,  | . . Hence  D"  is  to  the  left  of  A" 

1 at  D 1 1 at  A 

in  Figure  2-13(c),  which  is  the  template  if  . Summarizing, 

if  the  template  of  has  the  shape  of  Figure  2-11  (a)  in  R.,  , 

then  P^e  has  the  shape  of  Figure  2-11 (c). 

In  this  case,  typical  templates  of  Lg  = GgPjP2  (denoted  by 
®7)  in  single  loop  design  and  Lq  «*  GPieP2  (denoted  by  ^7')  in 
2- loop  P.M.  design  are  shown  in  (d),  (e)  of  Figures  2-12,... 15. 
These  are  constituted  from  the  corresponding  P^  , P2  and  P 
in  each  figure.  The  shaded  regions  in  7 and  represent  TJ  at 
lPlLin  due  to  P?  (denoted  by  72>  and at  |p  | . due  to 


1 1 min 


le 1 min 


y2  (denoted  by  respectively.  These  templates  “J,  °J' , and 
®j2  in  Figures  2-12, ...15  are  used  to  justify  (13). 


There  are  five  possibilities: 


'p.U  + L )•  " |p,  [ , 'i  + l 

1 s max  1 1 1 min  s 


(2.4-14) 


(because  [ ^ L~-|  occurs  at  I Pj  I mln  >» 
s max 


and  similarly 


1 P,  (1  + L )'  “ Ip,  [ . 'l  + L 1 

le  o max  1 le'min  o 


(2.4-15) 


Then  both  sides  of  (13)  equal  1 and  (13)  is  satisfied  with  ’ - ’ 
sign.  The  situation  is  shown  in  Figures  2-18,19. 
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-180*  -120*  -60*  -180*  -|20* 


Figure  2-19  °j'  (<*>A)  at  |AT(jo»)|  - 
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. . ..  . 


(2a) ; (14)  occurs,  so  £hs  of  (13)  - 1 as  before,  but  not 

so  for  Lq/  [P le (1  + Lq)  1 , from  which  it  follows  that 


(1*L„) 


max 


le 1 min 


L 

o 


1 + L 


o max 


(2.4-16) 


So  (13)  is  satisfied  with  ' < ' sign.  This  situation  is  shown 
in  Figs.  2-20,. . .22. 

(3a) : Neither  (14)  nor  (15)  occurs  and  such  that  for  each 

sides  of  (13),  the  numerator  7X.  is  > the  denominator  c 6.  However 
(^/^))^^s  < (^$)r^g,  because  of  the  shift  of  the  template  of  P 
to  the  left,  previously  emphasized  (this  is  to  be  proven,  of  course). 
Examples  of  this  situation  are  shown  in  Figures  2-23,... 25. 

We  also  have  to  prove  that  cases  (2b, 3b)  below  do  not  occur. 


(2b): 


‘Pjd  + L 


s max 


min 


1 + L 


s max 


(2.4-17) 


but  (15)  is  valid  for  the  right  of  (13). 

(3b) : Neither  (14)  nor  (15)  occurs,  so  for  each  side  of 

(13),  7l>«£).  But  (#/-0)ehs  > (>?/.0)fli8  . 

Case  1 will  occur  in  the  lower  u>  range  when  both  | | and 
|Lg|  are  so  large  as  to  give  the  situation  shown  in  Figures  2-18,19. 
Here  I L„/ d + Ls)  lmax  ls  Pol«t  X , lLs/d  + Lg)  lmax  at 

‘ *>2 

lPlLin  18  at  Y"X*  IV^  + VLax  18  at  X’  and 

lLo/(1  + Vlmax  at  lPleLin  ls  at  Y’  " X’  • So  both  <14> 

*2 

(15)  occur  and  1 - ’ sign  in  (13)  is  valid. 


-120' 


I at  (jw)  | 


. . 


Case  2a  tends  to  occur  in  the  middle  u)  range  giving  the 
situation  shown  in  Figures  2-20,... 22.  Here  points  X and  Y 
coincide,  but  not  X’,  Y’  . So  (14),  (16)  are  valid  and  ’ < ' 
sign  in  (13)  is  valid. 

Case  3a  (Example  Figure  2-23) 


s max 


is  at  point  X , which  does  not  correspond  to 


1 1 min 


.At  |P 


l^min  max^mum  *-s  ^ue  to  point  Y . If  we 


try  a larger  | P. | , 


is  increased  but 


1 + L 


1 Ls  P2Gs 

decreases.  Since  — • -+  ^ ^ 4.  p'  p 'q  ' » the  maximum  over  of 

1 s 12s 


(max  I1--,--;8- — r] } can  be  found  by  considering  the  loci  of  constant 

v2  v1+v 

|l  + L|  in  Figure  2-5.  In  the  entire  shaded  region,  |l  + L| 
increases  as  |l|  is  increased.  Hence  |y-yy|  decreases  as  |l| 
is  increased.  Therefore  max  is  at  point  Y.  So  the  ratio 


|P.  | . >1  + L I 

1 1 1 min  s max 


I § 

'Pl(l  + L8) 


• {Mag.  of  ly---L—  1 at  pt.  X} 


Tp'T  * {MaR*  of  ly+y  ~1  at  pt.  Y} 
1 ' min  s 


2 db  - 1.4  db  - 0.6  db 


(2.4-18) 


Now  consider  rhs  of  (13).  The  equivalent  points  X,  Y are 

L 

point  X’,  Y'  and  again  |p — . | is  at  point  Y*  , giving 

le'1  o'  max 


<•'  ■ .. 


TpH  < ’ |ttH  Tph- . * {Mag- of  Ii+l"I  at  ptX'} 

1 J.6  min  mav  1 1 o 'min  a 


1 P,  (1  + L ) 1 
le  o max 


le 'min 


lPleLin 


u 

{Mag.  of  [y +°L  | at  ptY'} 


2 db  - (-0.7  db)  = 2.7  db  (2.4-19) 


Note  this  value  of  2.7  db  > 0.6  db  in  (18),  and  the  reason  is 
due  to  the  leftward  shift  of  previously  emphasized.  This 

creates  a bigger  "effective"  separation  between  X’ , Y'  than 
between  X,  Y . This  case  also  occurs  in  Figures  2-24,25  and  it 
will  be  seen  later  that  this  is  so  in  general  in  R2  , when  the 
template  of  has  the  shape  of  Figure  2-11  (b). 

Cases  2b, 3b  do  not  occur  for  exactly  the  same  reason  as  in 
case  3a,  due  to  the  leftward  shift  of  ' creating  a bigger 
"effective"  separation  between  X',  Y'  than  between  X,  Y. 

Consider  the  second  kind  of  pattern  of  template  P 

(Figure  2-llb).  Typical  templates  of  P^  P2  at  coA  < u>B  < coc  < 

are  shown  in  (a)  and  (b)  of  Figures  2-26,..., 29.  With  the  same 

typical  locus  of  = P^H^  in  the  frequency  range  R2  in 

Figure  2-16,  the  resulting  templates  of  P^g  = P^/d  + L^)  , for 

P^  in  (a)  of  Figures  2-26,... 29  are  shown  in  (c)  of  the  same 

figures  and  this  shape  of  P^g  is  next  explained.  Again,  we  can 

suppose  P^en  ” P^n/(1  + L^n)  w*-ch  error  62  very  small.  In  (a) 

and  (c)  of  Figures  2-26,... 29  the  nominal  P,  , P,  are  at 

In  len 

A'  and  A"  respectively.  Note  that  both  templates  P^,  P^e  extend 
upward  and  to  the  left,  but  P^e  bends  more  to  the  left  than  that 
of  P^  in  Figures  2-26,... 28.  But  P^g  bends  less  to  the  left 
than  that  of  P^  in  Fig.  2-29.  This  can  be  explained  by 


5 


Figures  2- 30(a),  (b).  Very  similar  to  Figure  2-17,  ABCD  in 
Figure  2- 30(a)  corresponds  to  A'B'C'D'  in  Figure  2-27(a).  Note 
XD  = 1,  so  the  vector  from  X to  any  point  in  ABCD,  represents 
1 + L^  and  Arg  > 0 in  R2.  Arg(l  + Li)  ^ at  A = - ZoXA  and 

Arg(l  + Li)“1  at  C = -Z0XC.  Obviously , Z OXC  > Z OXA,  so 
Z Ple|at  C - Z ple|at  A >Zp1|at  C -Z-P-Jat  A.  Hence,  C"  is 
more  to  the  left  of  A"  than  is  C'  with  respect  to  A' , as 
shown  in  Figures  2-26, ...28.  On  the  other  hand,  if  |L  | is  big 
enough,  e.g.  |L^|  “ I (^p)  | > then  ABCD  in  Figure  2-30(B) 
corresponds  to  A'B'C'D'  in  Figure  2-29(a).  Note, Z OXC  < Z OXA 
here.  So  Z P^g | at  C - Z P^e | at  A < Z | at  C - Z P^ | at  A.  Hence 
C"  is  less  to  the  left  of  A"  than  C'  with  respect  to  A' , as 
in  Figures  2-29(a),(c).  Fortunately,  this  is  the  case  with  large 
|Lin|,  in  which  u)  is  large.  So  if  it  does  occur,  it  is  only  at 
the  very  end  of  f^.  In  most  cases,  it  will  not  occur  in  R2* 

Concluding,  if  the  template  of  P^  has  the  shape  of  Figure  2- 11(b) 
in  R2»  then  P^e  also  has  the  same  shape  but  wider  than  that  of 

pi-  ; 

Similarly,  typical  templates  of  Lg  = GgPjP2  (denoted  by  ) 
in  single  loop  design  and  Lq  = GP^eP2  (denoted  by  *J')  in  2- loop 
P.M.  design  are  shown  in  (d) , (e)  of  Figures  2-26,... 29.  The 
shaded  regions  of  and  represent  "J  at  |P]Jmin  due  to 

^2  (denoted  2)  and  at  !Pielmin  due  t0  v2  (denoted  by 

t 

2).  The  five  possibilities  in  relation  (13),  as  pointed  out 
before,  are  considered. 


Cases  1,  2a,  2b  do  not  occur,  because  the  shapes  of  P^,  P2> 

P,  are  as  in  (a),  (b) , (c)  of  Figures  2-26,... 29  and  in  (d) , (e) 
le 


of  Figures  2- 26,... 29;  the  rhs  of  TJ  (or  °J  ' ) coincide  with  the 
» 

rhs  of  *7.  (or  *J~).  But  |L/(1  + L)|  occurs  at  S,hs  or 
2 2 ' 'max 

bottom  of  *7,  rj'  (see  Figures  2-31, . . . 37),  so  neither  (14)  nor 
2 

(15)  can  occur. 

Case  3a  (Figures  2-31,... 36) 


I is  at  point  X , which  does  not  correspond  to 

s max 

lPll  * lPllmin  • At  lPlLin  the  maximum  of  ItmH  1s  at  Y 

1 Ls  P2G2 

Similarly,,  — ^ + L ■ y yy - and  by  considering  the  loci  of 
1 s s L 

constant  |l  + L|  in  Figure  2-5,  | — ■.  £ - — r-|  is  at  point  Y. 

1 s max 

So  the  Jlhs  of  (13)  is 


fpTT 4 * 'ttH  jh  < * {Ma8‘  of 

I l'min  s max  1 l'min 

II  - — — - - ' - - — ■ as  -■■■■■  — - 


'p.U  + L )' 

1 s max 


TjTT  * (Mag.  of 
1 l'min 


(2.4-20) 


And  similarly,  the  yhs  of  (13)  ls 


. i ^ 

|P,  | . 1 + L 

' le'min  r 


p.  (1  + L ) 1 
le  o max 


• (Mag-  °f 
• {Mas-  °f 


(2.4-21) 


For  example,  in  Figure  2-31,  (20)  gives  1 db  - 0.32  db  - 0.68  db  and 
(21)  gives  1 db  - 0.32  db  = 0.68  db,  so  '='  is  valid  in  (13). 

And  it  is  also  so  in  Figure  2-32.  In  Figure  2-33,  (20)  -*•  2 db  - 
0.9  db  “ 1.1  db  and  (21)  -*■  2 db  - 0. 7 db  ■ 1.3  db,  so  '<'  is  true 
in  (13),  and  also  in  Figures  2-33, ...36.  The  reason  is  that  °J  ' 


in  the  2-loop  P.M.  case  is  wider  than  the  corresponding  in  the 
single  loop  case. 


82 


Case  3b  (Figure  3-37) 

Here  (20),  (21)  are  also  true,  but  as  ' is  narrower  than  , 
(20)-*  2 db  - 0 db  - 2 db  and  (21)  -*•  2 db  - 0. 7 db  = 1. 3 db.  Note 
this  occurs  at  higher  oj  and,  fortunately,  it  is  at  the  beginning 
of  R3  or  at  worst  at  the  very  end  of  R2  , where  the  signal  level 
is  not  important. 

Actually  the  third  pattern  of  Pj^  in  Figure  2-ll(c)  is  the 
combination  of  the  first  two  of  Figures  2-ll(a),  (b),  so  the 
reasoning  and  conclusions  for  the  first  two  also  applies.  Summariz- 
ing all  the  above  discussions,  (13)  is  true  in  R2  for  any  kind  of 
template  of  P^  . So  (2.2-7)  becomes 


p(u>)  < 1 


le 'min 


"1 1 max 


(2 . 4-22) 


(2.4-23) 


where  a £ 1 in  R2 

Recall  that  I \nLJ  I ple  lmln  " «2 1 1 + Linl  • where  «2  «*  1 in  R2 

Vl  V, 

Since  a < 1 in  Rj  , the  resulting  design  is  conservative 
if  it  is  based  on  |l  + L^n|  . The  factor  a is  shown  in 
Figure  2-9  for  two  designs  and  it  is  seen  that  a < 1 in  R„  . 
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-180*  -120*  -60*  /Ls  -180*  '120*  '60*  £Lo 


Figure  2-37  GJ ( Wp),  (u>D)  in  Nichols'  chart 

with  |AT(ju))|  * 24  db 

A very  Interesting  point  is  that  based  on  such  a conservative 
approach,  one  can  design  a 2-loop  P.M.  system  with  the  signal  level 
variation  Q defined  in  (1.6-25)  even  less  than  1.  This  will  be 
seen  in  Chapter  3. 

2.4.3  R3  M (b>2 » u)^] 

Recall  R_  covers  {to I - M„  db  > |L  (jto)  I > - (M  +M  )db) 
where  M^,  Mj  are  the  hf  uncertainties  of  and  P2  respectively. 

In  problems  with  large  parameter  uncertainty,  R^  is  well  above 


i 
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the  control  system  bandwidth,  so  the  control  strength  is  small  in 
Rj  . We  can  therefore  tolerate  SLVR  >1  in  R^  . 

Since  R^  is  relatively  high  in  frequency,  it  is  reasonable 

ei 

to  assume  that  over  most  of  R^,  P^  -*■  k^s  , where  e^^  is  the 
excess  of  poles  over  zeros  of  , ki  is  the  hf  gain  of  Pj  , e.g. 
in  Figure  2-6,  R^  ■ (8,44],  while  the  maximum  control  system 
bandwidth  allowed  (see  Figure  1-2 (b))  is  2.1.  The  single  loop 
transmission  Lgn  definitely  has  reached  the  universal  hf  cylinder 
boundary  B^a  (corresponding  to  y - 2 db  in  this  example  - 
Figure  2-38). 

In  Figure  2-38,  the  typical  template  AF  is  used  to  discuss 
the  first  factor  in  (2.2-7).  Let  |Ls/(l  + Lg)|  = Mg.  Point  B in 
Figure  2-38,  corresponding  to  Mgmax  can  achieved  at  an 
infinitude  of  combinations  of  P^  and  P£  , giving  say  a set  S. 

i i * 

Let  the  smallest  |P^|  in  S be  denoted  by  P^  . Obviously,  from 

Figures  2-38,  for  any  fixed  P^  , Mg/P^  is  maximized  at  P^  = 

|P  | , and  with  respect  to  P at  | P | <_  | P | . Thus, 

z max  l l l 


GJP,I 


‘p.a+L  )' 

1 s max 


s'  2 1 max 

1 + G P_  P.  | . 
s 2max  1 1 min 


(2.4-24) 


From  the  constant  |l  + L|  curves  in  Figure  2-39,  it  is  seen  that 

to  minimize  the  denominator  of  (24),  |P^|  should  be  chosen  as 

small  as  possible  when  the  template  of  L at  |P-|  is  in  the 

s z max 

shaded  area  and  as  large  as  possible  when  L at  |P  | is 

s z max 

below  the  shaded  area.  Hence  the  minimum  point  is  at  D in 
Figure  2-39,  which  corresponds  to  point  D in  Figure  2-38.  In 


the  early  part  of  R0,  L at  |P„]  is  on  or  above  D V P.  , 

J r 3s  '2 'max  1 

so  | P^ | ” |P  I is  used.  But  eventually  in  R^  , the 
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template  at  does  not  extend  to  or  above  D,  so  I Pj | is 

chosen  big  enough  to  make  L reach  D at  |P„|  . The 

s 1 2' max 

required  |P^|  Increases  with  u>  , until  at  » the  end  of  R^  , 

it  is  I P I - AD.  In  this  later  part  of  R..  , L /(1  + L ) 

1 'max  v 3 s s 

corresponds  to  point  D and  |P^|  I Pi I min'  var*es  from  lp^lm£n 

to  l^j_lmax  ~ AD*  earlY  Part*  |Lg/(l  + Ls)|  lies  somewhere 

in  BD  while  |pj  - IpJ^  . so 


'1  + L 


s max 


1+L 


s B 


1 P. (1  + L )' 

1 s max 


T^T. 


min 


1 + L 


s in  BD 


‘ V |PlLin  : 6iel°  db*  |BD|dbl 

(2 . 4-25) 

where  | BD | ■ 7 db.  In  the  later  part,  | L^/ (1  + Ls ) | is  always 
at  D , so 


'1  + L 


s max 


P. (1  + L )‘ 

1 s max 


'1  + L 


s B 


|P1I*  varies  from  |P1lmln 

to  I P . I - AD  at  the  end 
' 1 ' max 

of  R3 


|P1ldb  + [2  db  - (-0.2  db)] 
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(2.4-26) 


I 


where  6 varies  from  2.2  db  to  2.2  db  + I P I - Ip  I _ am 
1 1 l'rnax  1 l'min  AU’ 

since  the  smallest  AD,  U it  u ■ Uj  , - 3 db.  So  6 e [2.2  db  , 

|AP1|  - 0.8  db].  Combine  (25),  (26) 


— I 

1 + L 1 


P . (1  + L )' 
i s max 


, 6 - [o  db,  Jap  |-o.8  db J 


(2.4-27) 


Consider  the  second  factor  in  (2.2-7).  A typical  outer  loop 
Lq  and  template  of  Lq  in  R^is  shown  in  Figure  2-40.  Just  as  in 
the  first  factor  (24). 


^ Lon 

TferftpUti 

of  Lo 


-40  1^- 
-3oo° 


l 1 1 

-180*  -60* 

£ Lo  ( decree; 


Figure  2-40  L__  and  typical  template  of  L in  R. 
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'P,  (1+L  )' 
le  o max 


|g|  • |p,| 

1 1 2 ' max 

1 1 + GP»  P,  I , 
2max  le'min 


(2.4-28) 


And  by  inspection  of  the  constant  1 1 -f-  L.  | curves  of  Figure  2-39, 
the  corresponding  Lq  at  ! P£ | m^x  , |Plelmln  in  R3  , is  in 
A'B'C'.  In  the  early  part  of  R3  , GP2inaxPl  ls  above  B»  so 


'P,  (1  + L )' 
le  o max 


le 1 min 


o in  AB 


(2.4-29) 


Thus,  the  second  factor  of  (2.2-7)  is 


|L  /[P.  (1  + L )]| 

1 o L le  o 1 ‘max 

|L  /(1  + L )| 

1 o o ' max 


1 

TTJ, 


le'min 


|Lo/(1  + Lo)'lnAB 

IV<1+lo)Iy 


3 1 le'min 


, 63e[Odb,  3 db] 


(2.4-30) 


In  the  latter  part  of  R_,  in  Figure  2-39,  L at  lpj  , 
|Plelmin  al°n8  PC,  so  |P^e|  is  chosen  big  enough  to  make 

Lq  reach  the  y locus  (see  section  1.2.2),  in  order  to  minimize 
the  denominator  of  (28).  The  required  | P^  | depends  on  the 
amount  hf  overdesign.  In  the  numerical  example  of  Figure  2-40 
(Design  A chapter  3),  with  5 db  hf  overdesign  on  Lq  , 

- lpJ.i„  + 5db-  s° 


WBxr  1 


1 iLo/<1+Lo>i 


L / (1  + L ) 

' o o' 'max 


Tf^T*  |Lo/<i  + i.o)[y 


TpTT*  ‘ «,Tp/|  , • V0db' Sdb' 

le1  4'  le'min 


(2.4-31) 


where  6 db  Is  the  hf  overdesign  of  Lq  . Combining  (30),  (31) 


iyiti.'in."i. 


max  j 

" &'  • Ip 1 • 6'e[0db,  max  {3db,  3 db } ] 

* 1 le'min 


(2.4-32) 


Substituting  (27),  (32)  into  (2.2-7), 


p‘  dT* 


?1  min  . C1  max 


le  min  ' s max 


a • 1 + L . 


'“I  max 


; a > 1 in  R. 


(2.4-33) 


(2.4-34) 


From  (33,34),  the  relation  |P  | /|p  | = 1 1 + L.  \ still 

l‘min  1 le'min  ' in' 

holds.  Although  is  not  small  in  (see  Figure  2-7),  the 

difference  between  I px lnin/  I ple !mln  and  1 1 + L±n | (=0.4db)  is 
fairly  small  compared  with  the  error  in  5/6'.  In  the  previous 
example  (Design  A,  chapter  3),  a e[0  db,  |APj  - 5.3  db]  = 

[0  db,  14.5  db].  And  in  Design  B (chapter  3),  a c[0  db,  14.5  db]. 
The  actual  results  for  a are  shown  in  Figure  2-9,  and  they 
confirm  the  above  argument. 

Note  that  a > 1 in  R^  , but  fortunately,  the  frequency 


I 


here  is  relatively  high  where  the  design  philosophy  is  to  emphasize 


I 

1 1 
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the  saving  in  sensor  noise  effect,  instead  of  considering  the  signaj 
level  variation.  In  addition,  if  we  use  an  additional  inner  loop 
around  P ^ (Figure  2-2),  the  a > 1 region  disappears.  This  is 
shown  in  Chapter  4. 


2.4.4  R4  = (to3,  u>4] 


Recall  R4  is  the  range  from  u>3  to  w4  = 10  u>3  , where  o>3 


is  the  frequency  at  which  L turns  the  bottom  corner  of  B.  at 

sn  hs 


X in  Figure  2-6.  In  R4  , Lg  is  under  B^g  , decreasing 
moderately  in  magnitude  and  considerably  in  phase. 

Consider  the  first  factor  of  (2.2-7).  In  the  early  part  of  7 

R,  , 


L /(1  + L ) reaches  the  forbidden  y disturbance  locus 

s s 1 max 


(y  = 2 db  in  Figure  2-6),  so  | L / (1 + L ) | * y db,  and 

s s max 


- G * Pj  / 1 + GP.  P,  | . . Consider 

' max  's'  1 2 1 max  1 s 2max  1 1 min 


iLs/[pi<i+Ls>r 

the  constant  1 1 -+•  L.  | curves  of  Figure  2-39.  In  the  early  part  of 
R, 


4 ’ 


1 1 + L | is  always  under  the  VWXYZ  line,  so 


'P.  (1  + L ) 1 
1 s max 


i L 

1 I S_| 

P,  '1  + L 1 

1 max  s ’ 


(2.4-35) 


so 


'1  + L 


s max 


'1  + L 


a Y 


1 P.  (1  + L ) 1 
1 s max 


1 1 max 


'1  + L 


s y 


1 1 max 


(early  part  of  R4)  (2.4-36) 
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r 


SI. ™ M-r3fflr  ■ '7^ 


In  the  latter  part  of  R,  , |L/(1  + L) 


s 'max 


cannot  reach  y in 


a practical  design,  but  occurs  at  |P, I , I P_ I , i.e.  at 

|LslmaX  which  for  example  Is  at  point  Q In  Figure  38  and  above 

WXYZ  in  Figure  2-39.  Consider  the  factor  I L / IP, (1 + L ) ] I 

1 8 1 s J 'max 

I®.!’  lP2L„/ll  + 0.P2™«Pll.ln-  “d  l1  + C.P2m.,Pll  1"”°t  * 

minimum  at  |P, I . However,  since  |L  I is  small  enough 
x max  e 

(Figure  2-6),  the  variation  of  |l  + L | due  to  P.  Is  fairly 
small.  In  order  to  simplify  the  explanation  later  on,  we  take 


P.  (1  + L )' 

1 s max 


T^T, 


max 


'1  +L  ' P.  , PJ 

s 1 l'max  1 2 'max 


in 


the  latter  part  of 


So 


ItTTTl 


s max 


1 1 + L P.  , P, 

s ' 1 ' max  ' 2 ' max 


'P.U  + L )' 

1 s max 


[77]  1 1 + L 1 I P | , | P | 

' 1 1 max  s 1 1 1 max  1 2 ' 


max 


- P 


Combining  (36),  (38) 


1 + L 


s max 


- | P.  | , in  R. 

' l'max  4 


'P.d  + L )' 

1 s max 


(2.4-37) 


-|  (latter  part  of  R . ) (2.4-38) 

I max  *♦ 


(2.4-39) 


Consider  the  second  factor  of  (2.2-7).  In  R. , L is  after 

4 on 

point  X'  In  Figure  2-6.  Just  as  in  the  single  loop  case. 


j . "»w  • ^z/9*^?****:. ■ '■;,w”'  *7'"  " ' ~7A 
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"P,  (1  + L )* 
le  s max 


le 


max 


1—2-1 
1 1 + L 1 


1+L  ' P,  , PJ 

o 1 le 1 max  1 2 1 max 


o max 


'1  + L 'P.  , 

o le  1 max  1 2 1 max 


• |P,I. 


1/ I P,  I , in  R, 

1 le 1 max  4 


(2.4-40) 


Substituting  (39),  (40)  into  (2.2-7) 


1 ‘max 

> 

le  max 


lCll 


max 


s 'max 


(2.4-41) 


Consider  the  factor  I P , I /|P,  I at  several  go  e R.  . 

In  this  range,  the  P^  template  is  certainly  a vertical  line  and 
L^n  has  negative  phase  angle  V P^  (see  Figure  2-7),  so  the 
template  of  P^e  - P^/(l  + L^)  definitely  looks  like  Figure  2-ll(a), 
and  this  is  so  for  our  large  plant  class.  A numerical  example  from 
Chapter  3 (Design  A)  is  shown  here  to  explain  the  relation  between 


lpil 


max 


and 


| Pie  I max  in  R4'  In  tflis  example  P^,  P ,,  each 


have  20  db  hf  uncertainty.  in  R^  is  shown  in 

Figure  2-41  with  the  template  P.  at  the  frequencies  co  , . . .co„. 

1 A r 


At  u>A  (-50)  , 


AP 


le 1 


0.2  - (-5.8)  - 6 db  (Figure  2-41),  is 


considerably  smaller  than  the  range  of  | AP 


20  db.  So  the 


difference  between 


1 1 max 


and  |P,  I is  14  db  - 
1 le 1 max 


Figure  2-42 (a).  At  coD  (-80)  , 

D 


in1 


larger  - 5.5  - (-11.5)  - 17  db  and 


is  smaller,  so  A|P^e|  is 

|P. I / 1 P-  I - 3 db  in 
1 l1 max  1 le'max 


Figure  2-42(b).  At  larger  go,  L.  is  close  to  -180° 

in 


(Figures  2-7,41).  The  template  of  P 


le 


10  - (-16)  = 26  db  at 


goc  - 100  and  A|P  | 


- 15.  5 - (-20)  - 35.5  db  at  to. 


125. 
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So  in  Figure  2-42  (c),  (pJ^/  |plel 


-6  db  . And  in 


Figure  2-42 (d),  I px I max/  I I*le I max  “ "15.5  db.  At  higher  u>  , 
Z.  Lin  < -180°  and  moves  downward  in  Figure  2-41.  So  lLinl 
very  small  and  P^e  - P^/d  + L^)  = P^  . At  uiE  = 250  , in 
Figure  2-41,  A|PX  | - -8.5  - (-30.5)  - 22  db,  and 


|P.|  / | P 

1 1 1 max  1 


le'max 


- 20  - 22  = -2  db  in  Figure  2-42 (e).  At 


oip  ■ 500,  in  Figure  2-41,  A|P^e| 


-24  - (-43.5)  = 19.5  db,  and 


|PlLax/lPleLax  " 20  " 19'5  = °-5  db  in  Figure  2"42<f>*  The  Plot 

of  P/[|C.|  / | C | ],  (which  is  |P.|  /|P,  | in  R.  -see 

1 l'max  1 s' max  1 l'max  1 le'max  4 

(41))  , for  Design  A,  is  shown  in  Figure  2-43.  Note  that  the 
results  in  of  Figure  2-43  agree  with  the  explanation  and 

Figure  2-42 (a) , . . . (f) . 

In  order  to  have  an  expression  similar  to  that  in  the  other 
frequency  ranges,  (41)  is  rewritten 


'1 'max 


(2.4-42) 


where  a = |P, I /[|P.  | 
1 1 1 max  1 1 le 1 


|P.|  / I p,  I 

1 1 1 max  ' le 


!1+Li„l) 


in  R, 


because  1 1 + L . I is  fairly  small  compared  with  I P , I /|P, 
in  R^  . This  can  be  seen  from  the  corresponding  in  the 

constant  |l  + L|  curves  of  Figure  2-44,  in  which  1 1 + L^l  varies 
from  1 db,  0 db,  -1.7  db,  -0.5  db,  0 db,  0.05  db  to  ^ 0 db. 
So  a in  is  close  to  I Pi  I lPie  lmax  * Tbe  nature  of  a 

in  R^  for  two  numerical  examples,  is  shown  in  Figure  2-45,  which 
is  replotted  from  Figure  2-9.  It  is  seen  that  in  R^  , a varies 


w1 » . 1 1 


100 


from  a large  positive  db  value  to  0 db  , and  then  to  a large 
negative  db  value  and  finally  tends  to  0 db. 

2.4.5  R,.  - (oj^  , *>) 

Recall  that  R,.  is  the  very  high  frequency  range,  where 

|L  I , |L  ] , I L . I « 1.  So  1 1 + L | “ 1 , |l  + L I - 1 

1 sn'  'on'  1 in'  ' s'  1 o' 

and  |Ple|  - IP^d  + L^I  - |PX|  . In  (2.2-7) 


P - 


s 'max 


L /P. 

1 s 1 1 max 


|g  p.p 

1 9 1 


2 ' max 


G P0 

1 s 2‘max 


|l  /p  | 

o le  max 
L 

o max 


| GP  | 

2 max 

gp,  p,; 

le  2 'max 


'1 1 max 


s ' max 


lCll 


max 


s 1 max 


lpi 


max  . 1 1 1 max 


[ p,  | |c  | 

1 le ' max  ' s ‘ max 


(2.4-43) 


Again  (43)  can  be  rewritten  to  have  the  same  form  as  in  the  other 
frequency  ranges 

lc.‘ 


p - a • 1 1 + L 


in' 


'1 1 max 


s max 


(2.4-44) 


where  a « 1 in  Rc 


The  curves  of  p/[|C, I / I C I ] and  a are  shown  in 
11  1 1 max  1 s 1 max J 

Figures  2-43,9  for  numerical  examples.  Note,  both  are  "1  in 


R,.  in  agreement  with  (43),  (44). 
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2.4.6.  Summary 

Concluding  the  above  discussions,  a £ 1 in  and  , 

where  the  control  signal  is  Important.  So  the  design  philosophy  is 
conservative,  a > 1 in  R^  and  in  the  beginning  of  R^  , a < 1 
in  the  latter  part  of  R4  and  a - 1 in  R$  . The  latter  are  the 
high  frequency  ranges,  where  the  signal  level  is  not  important  and 
we  concentrate  on  the  reduction  of  the  sensor  noise  effect. 


CHAPTER  3 


2-PLANT.  2-LOOP.  P.M.  SYSTEM  SYNTHESIS 


3.1  Introduction 


This  chapter  is  devoted  to  the  simplest  (two-loop)  plant- 

modification  (P.M.)  structure — Fig.  3-1.  The  plant  has  two  cascaded 

sections  P^,  P^  and  two  loops  are  used.  (A  more  complex  structure 

with  a third  loop  from  C ^ to  is  considered  in  chapter  4) . The 

loop  from  C to  X is  called  the  outer  loop  L ■ G P„P. /(1+P.H  ) 

1 li-  o 2 1 1 1 

A 

■ G P2Ple  an<*  ot^er  is  t*'e  Plflnt  modification  inner  loop 
Lj«P^H^.  At  best,  the  inner  loop  L^  can  be  designed  to  take  care 
of  all  parameter  uncertainties  in  P^,  as  in  the  cascaded  no  P.M. 
design  philosophy  (section  1.4.3).  But,  as  discussed  in  chapter  2, 
it  is  necessary  to  compromise,  because  of  the  resulting  increase  in 


signal  level  of  C^*  |L^|  is  therefore  made  small  in  the  low  fre- 

quency control  bandwidth  regions  R^,  R^.  So  the  inner  loop  L^ 


G Xn  P,  V*l  P,  \D* 


T A fl  _ V A p A A 

T R 1 + L * Lo  GP2Ple  ’ Ple  1 + L ’ Li  P1H1 
o l 

Figure  3-1.  2-plant,  2-loop,  P.M.  system 


It 
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cannot  handle  all  the  uncertainties  of  in  R^ , R^,  and  the  outer 
loop  must  in  this  sense  be  "overdesigned"  to  make  up  the  differ- 
ence. The  important  point  is  that  in  problems  with  significant 
uncertainty,  the  sensor  noise  effect  is  strongest  in  the  high  fre- 
quency range,  which  is  well  beyond  the  control  frequency  bandwidth. 

Based  on  this,  the  design  philosophy  is  to  concentrate  on  the  signal 
level  increase  in  the  i.f.  ranges  (R^,  R^)  and  let  L^.  do  its  best 
in  the  h.f.  range.  The  means  of  doing  the  former  was  prepared  in 

chapter  2,  with  (1  + L.  ) the  principal  design  tool. 

in 

In  this  chapter,  a systematic  design  procedure  based  on  the 

above  is  presented  with  detailed  numerical  examples.  Several  values 

of  signal  level  increase  ratio  are  chosen,  in  order  to  get  a feeling 

for  the  trade-off  between  the  bounds  on  I 1+L  I, increase  in  C„„ 

in  22 

signal  level  and  sensor  noise  effects.  For  comparison  purposes, 
single  loop  and  cascaded  2-loop,  no  P.M.,  designs  are  also  given, 
which  satisfy  the  same  specifications.  Also,  a hypothetical  P.M. 
design  which  ignores  the  signal-level-variation  problem  (ISLV)  is 
given,  in  order  to  see  the  maximum  savings  in  sensor  noise  effect, 
in  case  the  permitted  S.L.  increase  is  unlimited.  Finally,  a de- 
tailed comparison  is  made  between  the  P.M.  2-loop  system  and  the  no 
P.M.  cascaded  2-loop  system,  with  respect  to  sensor  noise  effects. 

It  is  assumed  throughout  that  each  plant  section  is  minimum- 
phase,  although  it  need  not  be  open-loop  stable  [2,  5]. 

_ _ - L J 


r 


Ej 
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3. 2 Principle  Steps  In  Design  Procedure 


(1)  Translation  of  time-domain  bounds  on  c(t)  into  bounds 
A 


on  | T ( j o>)  | = |c(jw)/R(ju>)  | , if  the  specifications  are 
not  originallyin  the  frequency  domain. 


(2)  Choice  of  bounds  on  1 1 + L^(ju))|  from  a given  value 


of  permitted  C ^ power  increase 


00  „ . 2 GO  O 

Q - / {C00(oj)}  dw/  / fC.  (w>rdw. 

O ZZ  o zs 


(3)  Determination  of  the  resulting  bounds  B^(w)  on  the 


inner  loop  L = t^ie  l°w  frequency  range, 


V V 


(4)  Design  of  inner  loop  in  R^,  R^  which  satisfies 


Bu<“>- 


(5)  Determination  of  the  template  of  P,  = P /(I  + Iv  ) 

le  1 i 


from  L^,  in  R^,  R^. 


(6)  Finding  the  bounds  B (a>)  on  the  outer  loop,  in  order 

Lo 


to  satisfy  the  design  specification  on  T(ju>).  In  R^ , 


R^,  these  bounds  are  chosen  so  that  is  satisfactory 


and  need  not  be  altered.  In  R„,  R, , Rc,  the  bounds 

3 4 5 


are  chosen  as  in  the  cascaded  design  philosophy. 


(7)  Design  of  the  outer  loop  L from  the  bounds  B,  (ui) , 

o Lo 


and  then  finding  the  corresponding  bounds  B (ta)  on  the 

LI 


inner  loop,  in  R^,  R^,  R,_,  such  that  chosen  is 


satisfactory. 


1 
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(8)  Completion  of  the  design  of  the  inner  loop  which 

satisfies  both  B ^ (w)  and  B (w). 

LI  LI 


3. 3 Design  Procedure 

The  design  procedure  is  explained  and  illustrated  by  means 


of  a design  example. 


3.3.1  Specifications  on  Numerical  Example 


Plant: 


Plant  uncertainty: 


Bounds  on  | T(  ju>)  | : 


P - PlP2,  - k1/s1 

/2f  k { 10  Jl,  /2  S k s 10  J2 

independently 

shown  in  Figure  3-2 (b)  were 
originally  derived  from  time 
domain  bounds  of  Figure  3-2(a). 


Disturbance  response:  Y £ 2.0  db 

(All  above  are  the  same  as  in  Section  1.2.5) 
Restriction  on  signal  level:  Q £ 1.05  for  unit  step  input. 


3.3.2  Translation  of  time-domain  specifications  (Figure 

3-2 (a))  into  frequency  domain  bounds  (Figure  3-2 (b)) 
This  effort  is  as  yet  an  art  rather  than  a science.  In  a 
minimum-phase  system  the  magnitude  of  the  frequency  response 
|T(ju>) | completely  specifies  the  transfer  function  T(s),  which  in 
turn  uniquely  determines  the  system  step  response  cs(t).  Hence, 
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ij 

i 


| ! 


1.2 

0.9 

c (t) 

0.6 

0.3 

0 

0 2 4 6 Time 

Figure  3-2 (a).  Time  domain  step  response  specification. 

bounds  on  £,n|T(ju)|  suffice.  But  the  rigorous  translation  of  time- 
domain  bounds  into  bounds  on  |T(jw)|  is,  as  yet.  an  unsolved  pro- 
blem. In  practice,  however,  it  has  not  been  difficult  to  achieve 
a translation  suitable  for  any  specific  numerical  problem  en- 
countered. One  may  begin,  for  example,  by  assuming  a simple  second- 
or  third-order  system  model  for  T(s),  and  finding  the  bounds  on  the 
model  parameters  which  correspond  to  the  bounds  on  the  time  res- 
ponse. From  the  model  parameter  bounds,  one  then  determines  the 
resulting  bounds  on  £n|T(joj)|.  Suppose  this  steps  leads  to  the 
solid-line  bounds  B , B of  Figure  3-3. 

It  is  desirable,  of  course,  to  Increase  the  spread  between 
B^  and  B^,  but  it  will  be  seen  that  there  is  no  advantage  in  doing 
so  at  Isolated  points.  There  is  benefit  only  if  the  spread  in- 
creases, on  the  whole,  with  increasing  0).  One  soon  finds,  with  a 
little  experimentation,  that  indicated  modifications  B^,  B2  in 


I 1 1 1 T 

upper-7" — 
bound/  v 

L V 

I lower  bound 
I / bp 
L / / 

/ / 

/ i 

►/ ! 


I 


i La 


i 


Figure  3-3  are  generally  achievable.  It  is  very  helpful  for  such 
experimentation  to  have  a computer  programme  for  finding  time  res- 
ponse from  the  magnitude  of  the  frequency  response.  Additional 
experimentation  reveals  that  there  is  a definite  limit  to  the  per- 
missible spread  in  the  lower  frequency  range  for  a reasonably  smooth 
curve  of  | T( J co) | . Subsequent  design  details  provide  one  with  an 
appreciation  of  the  frequency  ranges  in  which  broadening  of  the 
bounds  may  or  may  not  be  important.  Such  ranges  depend  a great  deal 
on  the  nature  of  the  plant  and  its  uncertainty.  Hence,  it  is  best 
to  obtain  comparatively  quickly  estimates  of  B^,  and  proceed  with 
the  design.  The  designer  will  subsequently  understand  whether  it  is 


A 


worthwhile  to  return  for  better  determination  of  bounds  on  |t(Jid)|. 
Practitioners  of  quantitative  design  have  had  very  good  results  in 
translating  and  satisfying  time  domain  bounds.  [2,  5-16] . 


3.3.3  Derivation  of  bounds  on  1 1 + L. ^ n | from  values  of 


permitted  maximum  Cjj  power  increase 


J [ | 1+L  |*|C  | ] 2 

. ..  , * „ A o in  2s  max  du) 

Q(design  criterion)  ■ Q ■ 

/ | C | du) 
o 2s  max 


It  is  recalled  that 


A / |C  I 2 dw 
A o 22  max 

/“K  | 2 do) 
o 2s 'max 


is  the  power  increase,  which  is  apriori  given  as  part  of  the  system 

specifications.  For  the  time  being,  assume  that  from  this  given  Q, 

we  are  able  to  determine  ■ 1.22  in  this  example.  The  procedure 

for  finding  Q from  Q is  given  in  section  3.5.  The  next  step  is  to 
c 

determine  the  |l  + L ( ju>)  | specification  from  Q . The  procedure 

in  c 

for  doing  so  and  for  finding  Q from  Q,  will  be  much  better  appre- 

c 

elated  after  a few  design  examples,  so  it  is  postponed.  In  the 
meantime,  we  accept  the  curve  shown  in  Figure  3-4  for 
[ | 1 + L^n | • | ^2s 1 mex]  vs  u>,  from  which  there  is  obtained  the  bound 
on  1 1 + Lin0w)  | in  Figure  3”5. 


f 


3. 3.  A Bounds  BT  ' (a>)  on  the  inner  loop,  LJ _,  in  R , 


The  constant  |l+L^|  curves  of  Figure  2-5  are  very  useful 

here,  as  they  give  directly  the  bounds  B,  . (w)  on  L.  «*  P,  H . For 

Li  in  In  1 

example,  in  Figure  3-5,  the  bound  is  1 db  at  ii)  * 1.3,  so  the  locus 

in  Figure  3-6  of  | 1+L  | ■ 1 db  immediately  gives  BP(1.3).  The 

1 LI 

shaded  region  is  the  acceptable  range  of  L.  at  id  * 1.3.  The  set 

in 

of  bounds  on  L shown  in  Figure  3-7  are  thus  obtained  in 

V V 


3.3.5  Design  of  L P which  satisfies  B.P,  (uj)  in  R, , R„ 


L (u>)  denotes  L.  (o>)  in  R . R„.  The  bounds  B fu>)  are 
in  in  1 2 Li 

shown  in  in  Figure  3-7.  Note  how  the  bounds  move  up  with  u>  in  the 
Nichols'  chart.  The  resulting  increase  in  |L^P  | is  accompanied  by 
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Arg  Lj>0  with  magnitude  dependong  on  the  rate  of  increase  of  |L  |. 

p 

The  larger  the  rate  of  increases,  the  larger  is  Arg  L^,  and  this 

results  in  a wider  shift  of  the  template  of  P^e  - P^/d+L^)  to  the 

left  (decreasing  angle) — Figure  3-9.  The  effect  on  the  outer  loop 

bound  is  to  shift  it  to  the  right  (Figure  3~10) , so  that  the  outer 

loop  | | decreases  more  slowly — which  is  not  good  in  the  effect 

of  sensor  noise.  This  is  discussed  in  section  3.3.6,  7.  So  we  try 

to  shape  L.^  to  satisfy  B,^.  with  smallest  ^L.P  , from  w * 0 to  u 
in  Li  in  2 

(at  end  of  R ),  with  f (0)  ■ 0°  and  f (u>  ) - 0°. 

2 in  in  2 

p 

The  resulting  L is  shown  in  Figure  3-7  and  its  Bode  plot 
in 

in  Figure  3-8. 


3.3.6  Templates  of  P1  m 

p 

As  noted,  L.  does  not  contend  with  P,  uncertainty  in  R , R . 

in  12 

So  the  outer  loop  L must  do  this.  The  total  uncertainty  template 
o 

P1  A A 

is  now  that  of  P • — ■ ■ - P„*P-i_.  The  template  of  P,  * P,/(l+L  ), 
2 1+L^  2 Ae  le  A i 

can  be  found,  because  is  basically  known  in  R^,  R^.  Some  re- 
sulting templates  of  P P1p  ^or  t^'^8  exaraPle  are  shown  in  Figure  3-9. 


3.3.7  The  outer  loop  bounds  B,  (u>)  in  R, , R0 

Lo 1 2 

The  outer  loop  bounds  B,  (w)  can  be  derived  by  hand,  bv 

Lo 

manipulating  the  templates  P,  P_  to  find  the  position  of  L (ju>) 

le  2 o 

which  results  in  the  specifications  of  Figure  3-2 (b)  on  (Ln]T(juO| 


\ 
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being  satisfied.  The  details  are  the  same  as  in  the  single  loop 
system  discussed  in  section  1.2.1,  because  it  is  an  identical  pro- 
blem. This  is  pretty  fast  but  it  can,  of  course,  be  done  entirely 
on  the  computer.  The  results  in  R^,  R2>  are  shown  in  Figure  3-10. 

Note,  in  Figure  3-10,  the  left  hand  side  of  B (w)  do  not  reach  the 

Lo 

disturbance  forbidden  cylinder  in  R^,  R^.  This  is  because 

L.*5  >0  in  R , R in  Figures  3-7,  8,  so  the  templates  of  P,  P„  in 
in  12  le  2 

Figure  3-9,  all  go  upward  and  to  the  left  from  P,  (the  nominal 

len 

point  of  P,  as  defined  in  section  2.2.1).  Obviously,  this  kind 
le 

of  B,  (ui)  forces  L to  be  further  to  the  right,  i.e.,  with  larger 
Lo  on 

phase  angle  thereby  decreasing  the  rate  of  decrease  of  |L  ( jio)  j 


and  increasing  the  sensor  noise  response  in  hf  region.  This  is 
the  price  we  must  pay  for  restricting  SLVR  in  R , R^. 


3.3.8.  Design  of  L and  determination  of  B , (oj) 

on Li 

The  outer  loop  L is  designed  to  satisfy  the  B,  (u>)  ob- 
on  Lo 

tained  in  section  3.3.7  in  R^,  R^*  and  to  be  pretty  much  like  in 

an  ordinary  no  P.M.  cascaded  design,  for  R_,  R, , R . i.e.,  with 

3 4 5 

a safety  margin — chosen  as  5 db  in  this  example — see  Figure  3-10. 

It  is  Important  here  to  note  that  in  R„,  R. , R„,  L is  designed 

3 4 5 on 

(except  for  5 db  margin)  as  if  there  is  no  uncertainty  in  P,  — 

le 

as  in  the  no  P.M.  cascaded  philosophy. 

The  resulting  L design  is  final.  The  next  step  is  to  ob- 
o 

tain  the  inner  loop  bounds  B.  . (w)  in  R„,  R, , R„ , such  that  the  L 

Li  3 4 5 on 
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design  is  satisfactory,  despite  its  neglect  of  uncertainty  in 
R^>  R^,  R^.  (This  part  is  the  same  as  in  no  P.M.  cascaded  design.) 

Also,  to  be  absolutely  sure  that  the  final  design  satisfies  the  fre- 
quency specification  of  Figure  3-2(b),  this  is  also  done  in  R^,  R^ 

This  precaution  is  worthwhile  because  L as  originally  designed 
for  R^,  R^,  must  be  modified,  albeit  very  little.  The  procedure 
is  now  exactly  the  same  as  in  the  cascaded  multiple  loop  (no  P.M.) 
system.  It  can  be  done  either  graphically  as  in  [2]  or  by  a com- 
puter search  program.  Some  B corresponding  to  L of  Figure  3-10, 

Li  o 

_l 
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are  shown  in  Figure  3-11.  Note  the  patterns  of  B for  u i 40  is 

LI 

very  similar  to  those  in  the  no  P.M.  cascaded  system  (see  Figure 
1-19).  This  is  so  because  the  main  difference  between  the  P.M. 
and  the  no  P.M.  is  at  low  frequency. 


3.3.9  Completion  of  L 

Recall  (section  3.3.5,  Figures  3-7,  8)  that  L ^ (fromu*  0 

in 

to  ■ 10),  was  obtained  as  the  first  step  in  design  of  L^,  i.e. 

in  R, , R . The  final  step  is  to  complete  L so  as  to  satisfy  the 
12  in 

bounds  B , some  of  which  are  redrawn  in  Figure  3-12.  Note  that  we 
Lil. 

have  deliberately  designed  |L^(j  | to  be  near  the  maximum  value 
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needed  by  L ■ — approximately  6 db  at  W » 55.  So  from  u>  = = 10 

to  to  » 55,  L.  must  be  approximately  horizontal.  Since  LP  was 

in  in 

designed  (see  Figure  3“8)  to  be  approximately  horizontal  for  to  > » 

it  is  relatively  easy  to  achieve  a good  fitting  of  LP  with  an  L 

in  in 

which  satisfies  the  bounds  for  all  to,  i.e.,  the  low-frequency  poles 

p 

and  zeros  used  in  L are  retained  in  L.  , which,  of  course,  has 
in  in 

more  poles  and  zeros  in  the  higher  range  in  order  to  achieve  the 
shaping  needed  to  satisfy  the  B (to) . 

LI 

The  resulting  L is  shown  in  Figure  3-12.  Comparison  with 
in 

other  designs  is  done  later — Figures  3-28,31.  These  other  designs 
are  next  presented. 


3.4  More  Numerical  Examples 


3.4.1  Introduction 


Sections  3. 3. 2-9  presented  the  details  of  a design  for  the 
problem  presented  in  section  3.3.1.  We  now  present  with  much  less 
detail,  more  designs  for  the  same  problem,  differing  chiefly  in 
the  signal  level  ratio  Q and  also  in  system  structure.  This  is 
done  in  order  to  give  the  reader  a good  perspective  on  the  effect 
of  Q and  how  the  design  results  compare  with  those  of  single-loop 
and  no  P.M.  cascade  structures. 


3.4.2  Single  loop  system  (Figure  1-3 (a 


Recall  section  1.2.5,  where  the  single  loop  transmission 


Z-Li 


-360 


l8C 

(Stole  23  -180° 


Figure  3**12.  The  inner  loop  L with  B 

in  Li 


L “GPP  for  this  numerical  problem  wate  shown  in  Figure  1-9. 
sn  s In  2 

Henceforth,  L is  used  for  comparison  with  the  P.M.  2-loop  system 
sn 

with  respect  to  signal  level,  sensor  noise  effect,  etc. 


3.4.3  Cascaded,  2-loop,  no  P.M.  system  (Figure  1-1 5 ( a)_) 


In  section  1.4.3,  the  outer  loop  and  inner  loop  L^ 


^2^2  casca^e<^»  2-loop,  no  P.M.  systems  (Figures  1-18,  19) 


were  presented  for  this  numerical  problem.  These  are  used  to 
compare  the  sensor  noise  effects  at  points  X^,  X^,  with  the  2-loop 


P.M.  design  (see  Figures  l-15(a),  3-1)  later  in  section  3.9. 
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3 . A . 4 1SLF  (Ignoring-signal-level-varlatlon) , 2-loop, 
P.M.  system  (Figure  3-1) 

In  this  system,  we  Ignore  the  problem  of  signal  level  in- 


crease. The  Inner  “loop  L^  ■ P^H^  is  desi8necl  in  precisely  the  same 

manner  as  L2  (“  P2H2  t*lere^  t^'e  cascaded,  no  P.M.  design  of 

section  3.4.3,  i.e.,  to  relieve  the  burden  on  L as  much  as  pos- 

o 

sible.  The  ISLV,  2-loop,  P.M.  system  is  compared  with  the  other 


2-loop,  P.M.  systems  to  s?e  the  sacrifice  in  sensor  noise  effect 
due  to  the  problem  of  signal  level  variation  in  the  P.M.  design. 


3.4.5  More  2-loop,  P.M.  designs  (Figure  3-1) 

Four  more  2-loop,  P.M.  designs  with  different  SLVR  require- 
ments, are  presented  here.  These  designs,  denoted  by  B,  D,  A,  E, 
were  achieved  trying  different  L ^ in  R^,  R^,  which  resulted  in 
different  Q value.  Based  on  these  results,  we  derive  a relation 
between  Q and  SLVR  in  section  3.5.  This  relationship  enables  us 
to  choose  SLVR  from  a specified  Q value.  Also,  the  sensor  noise 
effects  obtained  from  these  designs  gives  us  the  trade-off  between 
Q and  sensor  noise  effect  in  a 2-loop,  P.M.  system — section  3.8. 

p 

The  assumed  L.  of  designs  B,  D,  A,  E,  are  shown  in  Figures 
in 

3-13 (a-d)  and  the  corresponding | 1 + L.  lvalues  in  Figure  3-14. 

in 

(Curve  C represents  a design  done  by  deliberate  synthesis,  unlike 

A,  B,  D,  E. ) The  templates  P.  P„  are  shown  in  Figures  3-15  (a,.  . . 

ie  2 

d),  and  the  B and  L in  Figures  3-16  (a,  . . . d).  Finally, 

Lo  on 


' 


u)=o. 2 CO- 1 C0=2  co=3.2  <0=5  <o=8  <0=10 


(a)  Design  B 


<0  = 0.2  u)-l  c0=2  *0=3.2.  <0=5  u)=8  <0=10 

(b)  Design  D 


<0=0.1  u)=l  *0=2  <0=3.2  *0*5  (,0=8  <Os  10 


(c)  Design  A 


0°  iLLe 


-180° 


~9o* 


(b)  Design  D 


_ M ’ 
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3.5  The  Relation  between  Q and  SLVK 
3.5.1  Introduction 

Recall  (section  1.6.5)  that  this  work  Is  restricted  to  the 

RMS  problem  defined  by  Q.  But  the  design  criterion  is  based  on  Sl.VK 

(u>)  (section  2.2).  So  it  is  important  to  have  a method  for  choosing 

SLVR(ui)  from  a given  Q.  In  section  3.4.5,  four  designs  (B,  D,  A,  K) 

were  made  based  on  assumed  L f*  . These  results,  together  with  the 

in 

ISLV  design,  are  analyzed  here  and  used  to  derive  a relation  between 
SLVR(to)  and  Q.  Finally,  this  derived  relation  is  used  to  do  a 
synthesis  problem — which  was  presented  as  Design  C of  section  3.3. 

It  is  re-examined  here,  in  order  to  clarify  the  procedure  for 
choosing  SLVR(uj)  from  Q. 

In  this  chapter,  we  only  consider  Q and  SLVK  at  the  output  of 
(see  Figure  3-1).  The  signal  level  at  the  input  of  is  not  con- 
sidered here.  It  will  be  considered  in  chapter  4 in  considerable 

generality,  because  plants  with  more  than  two  stages  are  treated 
there. 


3.5.2  The  nature  of  C,„  and  C„  | 

1 22 'max 1 2s ‘max 

In  the  single  loop  system  of  Figure  2-1,  |c^(.1u>)|  ” 

|R(Joi)  • T^ ( Jto) /P ^ ( Ju))  | , where  T^Qw)  - (Jui) /R(Jia) , is  the  system 

transfer  function.  Similarly,  in  2-loop,  P.M.  system  of  Figure  3-1, 

|C22(Jw)|  - |R(Ju>)  T(ju))/Pu  (jw)|,  where  P}  ^ « P^/U  + l^) , - 

P,H  . |C  I 2 
1 1 ' 22  max 


of  Design  A is  compared  here  for  various  conditions 


of  the  same  problem.  Curves  of  (|l+L  | * | C | ) 


are  also  shown.  Figure  3-19  is  for  R - unit  step,  P 


Figure  3-20  shows  the  results  for  R - unit  step  but  P * type  0 


Figure  3-21  shows  the  case  of  R ■ unit  Impulse  and  P 


Finally,  Figure  3-22  shows  R ■ unit  impulse  and  P 


Note 


that  in  Figures  3-19  -22,  the  region  of  strength  of  the  signals 


changes  with  the  input  function  and  type  number  of  P . For  example 


in  Figure  3-19,  the  peak  control  signal  Is  near  1.2  rad/sec,  and 


in  Figure  3-20.  it  is  at  ui  - 0.  In  section  1.6.5,  the  power 


22 1 max 


Figure  3-19.  Comparison  of  signal  levels  R - unit  step,  P 
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variation  Q was  defined  to  be  the  ratio  of  area  under  the  curves 

1C  I 2 , |C„  I 2 • The  designer  knows  |C,  <»|  before 
*22  max  1 2s 'max  ma* 

connencing  the  P.M.  design.  He  can  use  these  figures  to  predict 
the  resulting  1C220“)Lx- 


Figure  3-20.  Comparison  of  signal  levels  R - unit  step,  " type  0 


00.  ,2 

3.5.3  Comparison  of  / C„„  dw  and  Q 

o 1 22 'max 


In  order  to  understand  the  relation  between  |l  + L,  I and 

1 in' 

Q for  various  conditions,  the  various  2-loop  P.M.  design  are  com- 
pared in  Tables  3-1,  2.  In  these  tables,  'criterion'  means  the 

OO  2 

data  are  calculated  from  / { 1 1 + L.  |*|C  I } dw,  and  the 

o'  in'  28  max 

00 . 2 

result*  column  contain  the  actual  / C,l  dio  for  the  single 

o ' 2s  max 

CO  - 2 

loop  and  f l^22^max^U)  ^or  t*,e  P-M.  system.  In  the  Q column, 

'criterion'  Indicates  Q - / { 1 1 + L I • I C I }2d w//  |c  I 2 dw, 

o'  in'  ' 2s  max  o ' 25 'max 

and  'result'  gives  Q - / |C  I 2 dw IS  Ic  I 2 dw. 

o'  22 'max  o ' 2s 'max 

The  ISLV  design  (ignoring-slgnal-level-variation)  may  be 
impractical,  because  the  signal  level  problem  is  ignored  and  the 
only  objective  is  the  reduction  of  sensor  noise  effect.  So  in 
Tables  3-1,  2,  its  Q etc.  values  are  the  largest.  But,  it  is  use- 
ful to  know  that  even  in  this  extreme  case  Q - 2.  Design  E for  R 
a step  input  (Table  3-1)  is  excellent  in  terms  of  signal  level 
variation  because  Q (result)  ■ .98<1,  although  Q (criterion)  = 

1.17.  This  is  because  the  design  criteria  are  very  conservative, 
as  noted  in  section  2.4.2. 

Next,  the  data  of  Q in  Tables  3-1,  2 are  plotted:  Q (result) 

= Q vs  Q (criterion)  - Q in  Figure  3-23.  The  very  interesting 
r c 

point  is  that  for  each  condition  (e.g.,  P^  ■ type  1,  R ■ unit  step), 
the  data  lie  on  a straight  line.  In  Figure  3-23,  a,  b,  there  has 


unit  impulse  function 


four  straignt  lines  representing  the  four  different  cases — 

(1)  ■ type  0,  R “ unit  step;  (2)  = type  0,  R=unit  impulse; 

(3)  P^  = type  1,  R * unit  step;  (4)  P = type  1,  R = unit 

impulse.  Also,  the  data  for  each  of  Designs  A,  B,  C,  D,  E,  give 

straight  lines,  as  shown  in  Figures  3-23(c)(d).  (The  lines  were 

obtained  by  linear  regression).  By  means  of  these  lines,  the 

designer  can  find  Q from  his  known  required  Q . With  regards 
c r 

to  C,  in  section  3.3.1,  the  specification  on  was  1.05,  for 

R = unit  step,  P = type  1.  In  Figures  3-23(b,d) , the  corresponding 

Q * 1.22,  which  was  used  in  the  design  in  section  3.3.3. 
c 


Figure  3-23(a).  Relation  between  Q and  Q — P =»  type  0. 

cl 
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3.5.4.  The  nature  of  L 


in 


For  pedagogic  reason,  it  is  better  to  understand  the  nature 

of  L in  R , R before  considering  how  to  derive  |l  + L (jai)  | 

In  12  in 

from  Q . The  nature  of  L is  decisively  influenced  by  the  nature 
c in 

of  non  P.M.  cascaded  loop  design  [5],  note  previously  in  chapter  1. 

Recall  that  in  the  If  range  (R^,  R^)  there  is  little  demand  on 

L^;  the  bounds  are  upper  ones,  i.e.,  L^  must  be  below  some  boundary. 

Consider  the  loci  of  constant  |l  + L^|  in  Figure  3-24,  which  is 

replotted  from  Figure  2-5.  We  want  it  to  be  small  in  R^,  R2 

because  it  is  being  used  to  determine  SLVR,  and  in  these  ranges 

we  want  small  SLVR.  We  also  know  from  our  experience  in  cascaded 

(no  P.M.)  design  that  at  higher  frequencies  (R^,  R^)  the  bounds 

change  to  lower  bounds  (L^  must  be  above  the  boundary)  and  the 

level  of  these  lower  bounds  gradually  increases  to  a maximum  at 

(x)  in  Figure  3-24,  ar.d  then  decreases  along  MNQ  in  Figure  3-24 
M 

(see  Figure  3-12).  In  our  P.M.  problem,  we  want  tc  retain  and 
improve  the  nice  property  of  the  non  P.M.  design  in  its  drastic 
reduction  of  sensor  noise  effect.  The  above  properties  of  the 
bounds  in  the  non  P.M.  design  continue  to  exist  in  the  P.M.  design. 
The  natural  tendency  for  is  to  reach  ui^,  at  ~ -90°  in 
Figure  3-24.  And  this  is  precisely  the  region  in  which  | 1 + 
is  large  ~ 5db  in  Figure  3-24,  giving  a large  SLVR  value.  However, 
it  is  not  this  large  SLVR  at  which  bothers  us  because  is 


in' 


relatively  large — in  R^,  in  which  the  signal  level  is  small.  It 


is  the  fact  that  we  must  reach  this  region  at  w starting  from  a 

M 

small  value  of  I L I in  R.  (Recall  that  small  I 1 + L.  I is  re- 
in' 1 in' 

quired  in  R^,  R?  because  the  signal  level  is  largest  in  R^ , R.,). 

So  we  have  to  shape  L in  R,,  R„  to  join  MNQ,  but  hopefully  so 

in  12 

as  to  achieve  small  jl  + L |.  Also  Arg  L (0)  must  be  a multiple 

in  in 

7T 

of  2.  Curve  ABC  in  Figure  3-24  seems  attractive  but  it  gives  a 

negative  half  encirclement  of  the  -1  point,  so  (1  + L ) has  a 

in 

rhp  zero,  unless  L.  has  a rhp  pole.  This  is  seen  from  the  Nyquist 
in 

criterion:  no.  of  encirclements  = N = N - N , where  N = no.  of 

z p 7. 

1 + L.  in  the  rhp  and  N = no.  of  poles  of  1 + L.  in  rhp.  Con- 
in  p in 
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1.  L.  * P H has  a rhp  pole.  This  pole  must  be  due  to 
in  In  1 


H,  and  then  it  appears  as  a rhp  zero  of  L 
1 on 


GP  P / 
2n  in 


(1+P  H ),  which  is  then  nonminimum  phase  and  is  very 
In  1 

limited  in  its  feedback  capabilities.  So  it  is  ex- 
cluded. 


2.  L does  not  have  a rhp  pole.  If  so,  it  is  impossible 
In 

for  such  a stable  L.  to  achieve  the  encirclement 
in 

shown  in  Figure  3-24,  because  in  such  a case  N = -1, 

N * 0,  so  N * -1  which  is  impossible 
P z 

A second  attractive  possibility  for  L in  R , R appears 

in  1 2 

to  be  A'B'C'  or  A"B"C"  in  Figure  3-24,  which  was  seriously  attempted. 
The  problem  is  achievability  of  such  curves.  Achievability  of  such 
curves  has  been  discussed  in  detail  in  section  7.7  of  [45] 

(Figure  7.7-4,  etc.)  and  it  was  shown  to  be  impossible  to  achieve 
the  kinds  of  ratio  Im^/M^I  generally  needed  here,  where  is  the 
maximum  |L^|  needed  and  is  * -20  db  due  to  the  hf  uncertainty 
in  P,. 


This  leaves  us  only  with  the  third  possibility  a"'B"' C"'  . 

Here  Arg  L.  (0)  * 0 and  I L (0) I should  be  small  in  order  that 
in  in 

|l  + L | is  small  in  this  range  where  the  signal  level  is  large. 

But  at  higher  to,  L.  must  merge  into  MNQ  in  Figure  3-24.  Hence, 
in 

Arg  L must  be  positive  over  much  of  R, , R„,  in  order  that  |l.  I 
in  12  'in' 

may  increase  from  its  low  value  at  <i>  = 0 to  the  value  needed  at 

M , denoted  by  w 
o o 


Sal 
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We  want  to  estimate  the  u>  value  at  which  L,  across  the 

in 

zero  phase  line  on  its  way  to  MNQ.  We  know  from  our  experience 

in  cascaded  design  that  | L I (i.e.,  point  M in  Figure  3-24) 

in  max 

occurs  near  J of  L in  Figure  3-25.  Our  experience  is  that  a 
on 

reasonable  estimate  of  w,/u)  ~ 3.  This  can  also  be  approximately 

J o 

theoretically  justified  as  follows.  Arg  L is  by  definition 

in 

positive  from  0 to  u and  |L  I is  at  ui  (pt  J) . Now  M (magni- 

o 1 in  max  M M 

tude  of  L.  at  point  M)  cannot  be  much  bigger  than  M , because 
in  o 

Arg  in  this  range  is  < 0.  So  must  of  the  increase  in  | L . | 
i in 

from  its  small  |l  (0)1  to  |l  (M) | must  be  achieved  in  (0,  w ). 

in  in  o 

From  0)  to  w , |L.  | will  be  fairly  horizontal  and  cover  a range 
o M in 

o 

of  - 90  decrease  in  phase.  This  is  somewhat  similar  to  the 

situation  from  I to  J in  Figure  3-25,  where  ~ 90°  is  also  covered 

and  where  the  range  U)  /<*)  ~ 3.  This  gives  cj  near  point  1 in 

J I o 

Figure  3-25. 


3.5.5  Choice  of  "criterion"  |l+L  | from  Q 
in c 

Recall  Q » A|l  + L |-|C,  | }2du)//°°|C0  | 2 du>  (■■=  1.22 

c o in  2s  max  o 2s  max 

2 

in  design  C).  I C _ | is  known  by  executing  a single  loop  design 
2s  max 

and  is  shown  in  Figure  3-26  (which  is  replotted  from  Figure  3-4). 

2 

So  the  area  under  the  curve  ( 1 1 + L.  |*|c  | } - M* (to)  should 

in  28  max 

2 

be  1.22  times  that  under  |c„  I . This  curve  should  be  such  that 

2s  max 

the  resulting  L (jw)  looks  more  or  less  like  A'"  B'"  C'"  (Figure 
In 


3-24)  in  the  region  R^,  R2  which  covers  the  significant  range  in 
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1.  L.  - P H has  a rhp  pole.  This  pole  must  be  due  to 
in  In  1 


H,  and  then  it  appears  as  a rhp  zero  of  L = GP„  P / 
1 r on  2n  In 


(1+P  H, ) , which  is  then  nonminimum  phase  and  is  very 
In  1 


limited  in  its  feedback  capabilities.  So  it  is  ex- 


cluded. 


2.  L does  not  have  a rhp  pole.  If  so,  it  is  impossible 
In 


for  such  a stable  L to  achieve  the  encirclement 
in 


shown  in  Figure  3-24,  because  in  such  a case  N = -1, 


N “ 0,  so  N ■ -1  which  is  impossible 
P z 


A second  attractive  possibility  for  L in  R , R appears 

in  1 2. 


to  be  A'B'C'  or  A"B"C"  in  Figure  3-24,  which  was  seriously  attempted. 


The  problem  is  achievability  of  such  curves.  Achievability  of  such 


curves  has  been  discussed  in  detail  in  section  7.7  of  [45] 


(Figure  7.7-4,  etc.)  and  it  was  shown  to  be  impossible  to  achieve 


the  kinds  of  ratio  jM^/M^I  generally  needed  here,  where  is  the 


maximum  L needed  and  M„  is  ■ -20  db  due  to  the  hf  uncertainty 
1 in  2 


in  P . 


This  leaves  us  only  with  the  third  possibility  A^B'"  C"'  . 


Here  Arg  L.  (0)  ■ 0 and  I L (0)1  should  be  small  in  order  that 
in  in  1 


|l  + L I is  small  in  this  range  where  the  signal  level  is  large, 
in 


But  at  higher  w,  L.  must  merge  into  MNQ  in  Figure  3-24.  Hence, 
in 


Arg  L.  must  be  positive  over  much  of  R , R , in  order  that  |L  | 
In  1 4 in 


may  increase  from  its  low  value  at  w ■ 0 to  the  value  needed  at 


M f denoted  by  co 
o o 


Figure  3-25. 


L in  Nichols 
on 


chart . 


Figure  3-26.  This  is  necessary  in  order  that  a higher  to,  (to 

achieve  significant  sensor  noise  reduction),  L.  satisfy  the 

in 

bound  at  and  then  follow  a curve  like  MNQ  in  Figure  3-24. 

So  the  primary  concern  at  this  stage  is  to  start  L at  to  = 0 

in 

at  a point  like  A'"  , and  arrive  at  a point  like  in  Figure  3-24 
(where  to  « 15).  Obviously,  there  is  no  unique  curve  of  Y (to). 


What  are  the  trade-offs  and  constraints?  It  is  easier  to  see  the 
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effect  of  L (ju>)  shaping  on  4*  (w) , rather  than  vice  versa  and 
In 

this  is  next  done. 

First,  we  can  use  the  Bode  integral  [4,  section  7.5] 


(3.5-1) 


A -A  - - — / Arg  L.  (j(a)d  £n  u) 
o « tt  O in  J 


A - en|L  <0)|-M  A « »n|L  (-)! 


The  interval  in  (1)  is  from  0 to  but  since  the  poles  and  zeros 
introduced  beyond  the  R^,  range  are  relatively  far-off  (e.g., 
in  our  example  C,  * 26.5  ± j45.9,  -115.5  t j 309 . ) , we  can  take 

-U) 

/ o , and  A = M . So  far  a given  value  of  M , 

o OO  o A 

/“o  Arg  Lln(jw)de,nu)  is  pretty  well  fixed.  This  is  one  semi- 
constraint (because  M ,,,  is  free,  and  so  is  w ), 

a o 

Suppose  we  have  found  a curve  which  satisfies  the  1.22 

condition  and  (1),  e.g.,  that  shown  in  Figure  3-26,  for  which 

the  resulting  in  R^,  R^  is  shown  in  Figure  3-27.  Consider 

keeping  the  same  shape  of  LJ  as  in  Figure  3-27,  but  shifting  the 

in 

pattern  as  a whole  to  the  left,  i.e.,  achieving  the  transfer  from 

A'"  to  C'"  in  Figure  3-24  in  a lower  frequency  range.  Clearly, 

this  would  cause  H' (u0  to  be  the  same  for  <a>w  , but  considerably 

o 

larger  in  between  0 and  w , thereby  violating  the  1.22  condition. 

0 

So,  if  we  want  to  "start  lb.  (1u))I  off"  sooner,  we  must  flatten 

1 in  ' 

somewhat  the  curve  of  Arg  L.  (jw)  in  Figure  3-27 — lower  its  peak 

in 

and  spread  it  out  over  a larger  range,  so  that  its  area  remains  the 
same,  as  shown  by  dashed  line  in  Figure  3-27.  This  will  increase 


mmmmmmm m ■ 


• _ ...  - 
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Y at  lower  od  and  decrease  It  at  higher  w.  Note  that  lifting  A"’ 

higher,  i.e.,  increasing  |L^(0)|,  means  /(i'o  can  be  less,  and  so 

peak  Arg  L.  (jiu)  can  be  smaller.  However,  in  Figure  3-24,  it  is 
in 

seen  that  this  tends  to  increase  | 1 + L I and  therefore  't'.  How- 

in' 

ever,  reduction  of  Arg  L (jco)  is  good,  especially  in  latter  part 

in 

of  R„ — because  the  bound  on  the  outer  loop  L is  shifted  less  to 
2 on 

the  right  (Recall  section  3.3.7  and  Figure  3-10).  So  this  kind 

of  effort  (starting  L (jto)  up  sooner  and  raising  A'"  ) is  good 

in 

and  one  should  do  this  as  much  as  possible.  But  the  1.22  limit 

OO 

on  / tdw  limits  how  much  can  be  done  and  that  is  basically  the 
optimum,  which  is  seen  by  considering  the  opposite  kind  of 
variation. 

Suppose  the  pattern  in  Figure  3-27  is  pushed  to  the  right. 

T (u>)  is  thereby  decreased  but  the  peak  in  Arg  L (jw)  occurs  at 

in 

a higher  u,  which  is  more  harmful  to  L (1w)  in  the  sense  of 

on 

sensor  noise  effect  (see  section  3.3.7,  Figure  3-10).  Thus,  one 
wants  to  move  in  the  opposite  direction — to  the  left — as  pre- 
viously seen. 

The  above  presents  the  trade-offs.  The  optimum  't'(ui)  curve 

would  require  some  cut  and  try.  Our  experience  has  been  that  a 

good  way  to  start  is  to  take  the  peak  of  [H'(u))  ] * M„  in 

max  X 

Figure  3-26  at  the  same  u>  value  as  the  peak  of  |C„„|  , take 

K 1 25 'max 

4*(0)  - |C„  I at  zero,  allow  a value  * .05  larger  than  |C„  I 

1 2s 1 max  B 1 2s 1 max 

in  the  higher  to  range  (see  Figure  3-26)  and  draw  a fairly  smooth 
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curve  satisfying  the  above.  Try  different  M values  until  the  1.22 

Jv 

condition  is  satisfied.  So  now  a trial  ^((o)  is  available.  Use  it 


to  determine  the  bounds  on  1 1 + L.  I and  find  the  resulting  L 

in'  in 

(i.e.,  L in  R , R ) wilch  will  look  like  that  in  Figure  3-24. 

In  1 Z 

Use  this  as  a starter.  Now  work  from  L.^  , rather  than  from  ¥(10) . 

in 

Push  the  |L  | increase  toward  lower  frequency  as  previously 

indicated  and  modify  Arg  L f and  possibly  the  A"'  point  in  this 

in 

quest.  We  have  found  this  to  be  a good  approach. 


3.6  Comparison  of  the  Loop  Transmissions 


Bode  plots  of  the  inner  loop  L.  and  the  outer  loop  L 

in  on 

for  the  five  different  P.M.  2-loop  designs  are  given  in  Figure 


3-28,  as  well  as  the  single  loop  transmission  L . Recall 

sn 

section  2.3,  R = (0.,  to  1,|L  (jto  ) ]=  25  db;  R = (to  , to  ] , 

1 1 sn  1 2 12 

|Lsn(j<o2)|«  -M2,  M2  » hf  uncertainty  of  P?;  R^  * (to7,  0)^] , | ( Jto^) 


2’  3 ' 2’  3 1 sn  3 


+ M2),  - hf  uncertainty  of  P^  R^  * (to  , to^] , (0^  - lOuy 

R,  » (to.,  °°) . 

5 4 

|L  | in  R , R , in  Figure  3-28  is  closely  related  to  SLVR 
in  1 2 

and  Q.  The  larger  the  large  are  SLVR  and  Q.  Note,  when 

L in  R is  lower  than  some  limit,  say  Design  E,  then,  although 
in  2 

the  criterion  on  SLVR,  1 1 + L . |,  >1  the  actual  p/ (|C  I / 1 C I ) 

1 in  1 max  1 s 1 max 

<1  in  Figure  2-29(e)  in  most  of  R2.  So  Q<  1,  because  of  the 
conservative  character  of  the  design  philosophy.  It  is  reempha- 


sized that  the  signal  level  in  R^,  R^,  R,_  is  not  important,  so 


! 


t 


' 

* 

s 


i 

li 

i 
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even  large  SLVR  (p) — there  poses  no  problem.  The  curves  of  I 1+L  I 

in 

and  P/(|C  | /|C  | ) for  Design  B,  D,  A,  C,  E are  shown  in 

1 UlciX  S IuaX 

Figure  3-29  (a),  ....  (e).  Note  the  values  of  p/{|C  I / 

1 max 

|C  | > in  Figure  3-29  in  R , R , are  directly  related  to  the 

S ulaX  X iJ 

value  of  IL^J  in  Figure  3-28.  For  example.  Design  B in  Figure 

3-28  is  the  largest  in  R , R so  the  value  of  p/{ | C I / I C I } 

12  1 1 1 max  s max 

of  Figure  3-29(a)  is  the  largest.  The  comparison  of  p/{|c  I / 

1 max 

|c  | } are  shown  in  Figure  3—29 ( f > . 


5 


n*i: 


[•I 


Q value 

- 1.41 

- 1.24 

- 1.10 

- 1.047 


Curves 


(f)  Comparison  of  p/(  C I / C I ) 

1 max  s max 

Figure  3-29.  Bode  plots  of  |l+L.  | and  p/(|C. | / | C | ) 

— B in  1 max  ' s max 


The  levels  of  L and  L in  R_,  R. , R.  in  Figure  3-28, 
on  in  3 4 5 

determine  the  sensor  noise  effect.  The  overall  tendency  of  the 

effect  of  sensor  noise  depends  on  the  levels  of  I L I,  I L,  I in 

Rj.  R^»  R^  relative  to  the  corresponding  value  of  1 . / j P2 1 and 

1 . / 1 | • The  equations  are  derived  in  (3.7-3,  4)  in  section  3.7. 

The  layer  Q,  the  smaller  the  price  that  must  be  paid  in  |L^|  in 

R, , R, . Obviously,  the  smaller  L , L,  in  R_,  R. , Rr,  the  smaller 
4 5 on  in  3 4 5 

is  the  sensor  noise  effect. 


Note,  in  Figure  3-28,  all  L have  5 db  hf  overdesign. 

on 


It  is  known  that  in  the  cascaded  (no  P.M.)  design,  some  over- 


design  of  L in  this  range,  makes  the  design  of  L much  easier, 
on  in 


This  is  basically  because  L.  can  be  thereby  reduced.  In 

1 in 1 max 

a P.M.  system,  this  is  equivalent  to  shift  L from  AB'C'D'E  in 


Figure  3-30  (with  I L . | at  M„'  = 0 db)  to  ABODE.  In  our 
' in 1 max  X 

example,  the  design  is  5 db,  so  |L.  I is  at  M »~5db.  Note 

• 1 in 1 max  X 

the  corresponding  |l+L^|  value  in  Figure  3-30,  is  much  smaller 
along  ABC  than  AB'C' . Recalling  section  3.5.5,  one  wants  to  start 


L.  (jto)  up  sooner  and  raise  A'"  (see  Figure  3-24),  to  achieve 
in 

lL,  I in  a shorter  to  range.  Obviously,  curve  ABC  in  Figure 
in  max 


3-30  is  better  for  this  purpose  than  AB'C'. 


3.7  Sensor  Noise  Effects 

Recall  the  sensor  noise  effect  is  evident  in  hf  range  where 
the  signal  level  and  disturbance  response  are  usually  small.  The 
noise  transfer  functions  are  as  follows: 

A.  The  single  loop  system — Figure  2-1,  sensor  N^ 


2 A X21 


-G 


1+G  P P 
s 1 2 


L /PP, 

s 12 
1+L 


L - G P P 
s s 1 2 


~ - L /P,P„,  in  the  hf  range  where  |L  |«| 
~ s 1 2 's'1 


(3.7-1) 


1 AX11 


-G  Po 
s 2 


N,  1+G  P,P„ 
1 s 1 2 


-L  / P 
s ] 

1+L 


s -L  /P  , in  the  hf  range  where  |L  |«| 
S 1 s 


(3.7-2) 
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2 1 

The  curve  of  T is  shown  in  Figure  3-32.  And  T < 0 db  for  all 

N1  N1 

(B)  The  2-loop,  P.M.  system — Figure  3-1  sensor  N 


_ 2 A X22 

r - — 

Ni  Nx 


i+pihi+pip2g 


(1+L  ) (1+L  ) ’ Li  = P1H1’  Lo  GP2  1+L 
I O 1 

-L  /P  P 

o 1 2 

(1+L  ) 
o 

~ Lo/P^P2’  *n  hf  range  where  | | << | 


1 AX12  -(GVV 


(3.7-3) 


Ni  Nx 


1+P1H1+GP1P2 


-IL  (1+L  )/P+L/P] 
o i 111 

(1+L.) (1+L  ) 
i o 


*>  (L  /P ,+L  /P  ) in  the  hf  range  where 
o 1 i 1 

|L1I«|,  |Lj«| 

* L^/P^,  1°  the  hf  range  where  |Lq|«|l  | 


(3.7-4) 


1 2 

The  curves  of|T  | and  |T  I for  the  five  different  2-loop, 
N1  n2 

P.M.  design  are  shown  in  Figure  3-31  (a)  in  DB  vs  log  u>  and  in 

Figure  3-31 (b)  in  arithematical  scale.  In  these  figures  Design  E 

is  the  worst  and  B the  best  case,  which  is  as  it  should  be  as  the 

SLVR  constraint  in  Design  E is  the  toughest.  This  also  can  be  seen 

directly  from  (3),  (4)  and  Figure  3-28.  In  (3),  (4),  | T 1 | and 

N1 

| Tn2  | depends  on  |lJ,  |lJ  and  Ip^I,  IpJ-  Since  ! 1>1P2 1 and  IpJ 


-180*  "120* 


6o*  I2cf 


/•Li 


Figure  3-30.  The  result  of  hf  overdesign  on  L 


are  the  same  in  all  designs,  the  values  of  L , L determine 

' o i 

'Vj  1^1- 

Compare  Design  E with  the  single  loop  design  in  Figure  3-32 

(a),  (b).  T of  signal  loop  <0  db  for  all  <a,  so  it  is  not  shown 
N1 

in  Figures  3-32(a,  b),  but  T of  Design  E >0  db  as  shown  in 

N1 

Figure  3-32(a,  b) . The  single  loop  design  is  better  than  the 
p.  P.M.  design  with  respect  to  sensor  noise  effect  at  the 
H"i'  •'  P • However,  with  respect  to  the  overall  effect  sensor 
system  is  better. 

• » « 1 1 » ■ ison  is  that  it  is  the  actual  sensor  noise 

» -r-,.  ..  Pj  which  are  important,  not  the 


comparison  between  noise  effect  at  in  the  single  two-loop 
designs.  To  illustrate  this  point,  suppose  that  in  the  single- 
loop design,  the  effect  at  (Figure  3-1)  is  given  by  a number 
10,  and  at  X^  it  is  .1,  whereas  the  signal  level  capacity  of  X^ 
is  of  the  same  order  of  magnitude  or  even  larger  than  that  of  X^ 
Hence,  it  is  the  noise  at  X^  which  may  render  the  design  imprac- 
tical. The  noise  effect  at  X^  is  insignificant.  Now,  in  the 
two-loop  design,  suppose  the  effect  of  X^  is  2 and  at  X it  is 
.5  (five  times  that  in  the  single-loop  design).  The  two-loop  de- 
sign is  obviously  considerably  superior.  Thus,  the  two  important 


ratios  are 


T 2 

Nj  (2-loop) 

T 2 (1-loop) 
N1 


and 


TN1  (2-loop) 
T 2 (1-loop) 


T ^ 

Nl_C2t.io°p)  ls 

T 1 (1-loop) 

N1 


not  important. 

The  above  is  precisely  the  situation  here.  C 


22 


X„  P 
22  2 


and  in  the  relatively  low  control  signal  frequency  range,  | | >l 


generally  in  this  range,  or  at  least  ~1,  so  lc22  ^*22 


Hence, 

the  signal  level  capacity  (i.e.,  which  the  output  stage  there  can 


comfortably  handle)  at  C ^ must  be  "i  that  at  X ^ • In  the  single- 
loop design  |T  2|  is  enormous  (see  Figure  3-32)  while  |T  ^<1.  In 
N1  N1 
2 1 
the  two-loop  design  |t  | is  considerably  reduced,  while  |T  |>1 


Ni1 


N, 


but  still  small  relative  to  IT^  | ^ and  certainly  very  small 

2 

compared  to  |T  |.  Hence,  the  2-loop  design  is  decidedly  better 


N 


1 


in  terms  of  overall  system  performance.  It  will  be  seen  in 


, , , • - 


chapter  4 that  this  superiority  is  considerably  increased  if 
another  loop  is  added  around  P . 


3.8  Trade-off  between  Q and  Sensor  Noise  Effect 


The  results  of  six  different  2-loop  P.M.  designs,  ISLV,  B, 

D,  A,  C,  E are  shown  in  Figure  3-33.  The  total  noise  effects  at 

and  (in  Figure  3-1)  are  calculated  by  means  of 

oo  i22  o).  2i  2 1 

/ T ($M  ) ]dw  and  / T ) ]doo  respectively,  where  T , 

° ± 1 o N^ 

T are  defined  in  (3-7-3,  4)  and  <J>  is  the  noise  spectrum  of 
N1  N1 
sensor  N , taken  as  a uniform  distribution  (4>  = 0.001).  In 

1 N1 
Figure  3-33,  the  two  curves  give  the  relations  between  Q and  sensor 

noise  power  response  at  X^>  X^.  T^ese  are  useful  for  giving  the 

designer  a feeling  for  the  trade-off  between  signal  level  and 
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noise  effect.  Of  course,  Design  ISLV  is  the  best  one  in  sensor 
noise  effect,  but  Q is  rather  high  (=1.96).  In  Design  B,  the  noise 
is  a little  larger  than  in  ISLV  but  the  improvement  in  Q is  fairly 
large.  In  Design  E,  Q is  excellent  (=0.98)  but  the  noise  is  signi- 
ficantly more  than  in  Design  C with  its  Q = 1.047.  It  is  again 
emphasized  that  even  in  Design  E with  Q<1,  the  noise  response  (1.59) 


w 


II 
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is  still  much  smaller  than  a single  loop  system  (5.63)  in 


Figure  3-33.  Trade-off  between  Q and  sensor  noise  effect. 


3 . 9 Comparison  of  Noise  Effect  between  the  2-loop,  P.M. 

System  and  the  Cascaded,  2-loop,  no  P.M.  System. 

For  a system  with  2 plant  sections,  if  the  design  is  re- 
stricted to  2 loops,  there  are  two  options:  One  is  the  cascaded, 
2-loop,  no  P.M.  system  (Figure  3-34(b>),  the  other  one  is  the 
2-loop,  P.M.  system  (Figure  3-34 (a)).  Suppose  the  choice  of 
structure  is  dependent  on  the  savings  of  sensor  noise  effect  which 


1 


V 


r 
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means  ignoring  the  fact  that  the  2-loop  no  P.M.  design  requires 
two  sensors  instead  of  one.  A fairer  comparison  is  therefore  the 
3-loop  P.M.  design  in  chapter  4 which  uses  two  sensors  of  the 
?-loop  no  P.M.  design. 


(a)  2-loop,  P.M.  system 


The  noise  power  ratio  (N.P.R.)  is  defined  as  the 
ratio  of  total  noise  power  in  the  cascaded,  2-loop,  no  P.M. 
system  to  that  of  the  2-loop,  P.M.  system. 


W 


+ IVN2|2-|N2|2 


(3.9-1) 


I 


Since  the  noise  response  of  a cascaded,  2-loop,  no  P.M.  system  is 


(see  Figure  1-15 (a)  and  (1.4-6)) 


2,  A , 2c, 


V Ni 


in  the  hf  range 


(3.9-2) 


X L 

lTN2l  * hr^l  * Ijrl  ln  the  hf  range 
2 W2  2 


(3.9-3) 


I Tm1  I * * l;rl  in  the  hf  range 

1 1 1 


(3.9-4) 


1 A i 

Itn  I * = I L2 1 ln  the  hf  range 

2 2 


(3.9-5) 


together  with  (3.7-3,  4).  The  NPR  in  (1)  becomes 


iy’1l2-i»1i2  + iy2-iy2 


lyyMy2 


(3.9-6) 


141*  * IV.IV  , 4 n2 

; , where  R - — 

IV  "i 


(3.9-7) 


lyoyylVy2  * iyy2-iy2 


(3.9-8) 


ILo/(PiP2)|2*INi|2 


I li 1 2 + I l2pi I 2*r2 


(3.9-9) 


Figures  3-35,  36  give  graphs  W?(u>),  W^oo)  for  different 


values  of  R,  |PjJ»  using  the  design  results  of  the  cascaded. 


2-loop  system  in  section  3.4.3  and  Design  A for  the  2-loop,  P.M. 


I 


■ 


3-35,  W^,  which  is  the  important  factor  (recall  section  3.7), 
is  <1  for  most  of  u). 
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CHAPTER  4 

MULTIPLE-LOOP  P.M.  SYSTEM  SYNTHESIS 

4.1  Introduction 

Chapter  3 dealt  with  the  simple  P.M.  structure  for  the  two- 
section  (P^,  P2)  cascade  plant.  A more  general  structure  for  this 
2-section  case  admits  three  feedback  loops — one  from  C23  to  X^  in 
Figure  4-1  may  be  added.  To  determine  the  most  general  feedback 
structure,  observe  that  no  feedforward  loops  are  allowed  (i  to  i+1 
in  Figure  4-2),  other  than  in  the  original  plant  in  Figure  4-2. 

In  this  research  we  consider  only  structures  in  which  there  are  no 
crossings  of  lines  when  the  feedback  loops  are  all  drawn  below  (or 
all  above)  the  plant,  e.g.,  in  Figure  4-3(a)  feedback  from  output 
of  P^  to  input  of  P2  is  not  allowed  because  such  a line  would  cross 
the  present  one  from  P2  output  to  P^  input.  How  many  different 
structures  are  possible?  We  shall  assume  that  in  each  structure 
there  is  a local  feedback  loop  around  each  plant  section  and  one 
from  n to  0 in  Figure  4-2,  as  these  do  not  prevent  any  other  loops. 

Let  S be  the  number  of  structure  for  the  n-section  of 
n 

Figure  4-2.  Either  there  is  a loop  from  node  n to  node  1,  or 
there  isn't.  If  there  is  one,  the  remainder  on  the  right  forms 


Figure  4-1.  2-plant,  3-loop,  P.M.  system 


Figure  4-2.  n stage  cascaded-plant 

an  n-1  section  with  number  S . The  remainder  on  the  left  has 

n-1 

on'.y  one  section  giving  S possibilities  with  a total  of  S x S . 

1 n— 1 1 

If  there  isn't  one  from  n to  1,  there  are  two  possibilities: 

Either  there  is  one  from  n to  2,  or  there  isn't.  If  there  is  one, 

then  the  remainder  on  the  right  is  a n-2  section  with  S and  on 

n— 2 

the  left  a 2-section  giving  S„,  with  a total  of  S x S . If 

2 n-2  2 

there  isn't  one  from  n to  2,  there  are  two  possibilities:  Either 
there  is  one  from  n to  3,  or  not.  If  there  is  one,  the  remainders 
on  the  right  and  left  are  respectively  n-3  and  3,  giving 
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^ x S^,  etc.  Thus,  in  an  n-section  with  no  "feedback  crossings," 
for  n>2 


S - S x S + S XS  +...+S  xS  +S  x S 
n n-1  1 n-2  2 n-2  2 n-1  1 

A A 

with  “ 1,  ■ 1.  Hence,  = 2,  ■ 5,  etc.  The  five  cases 


of  n*4  are  shown  in  Figure  4-3. 


Maximum  number  of  loops 

Each  of  the  possible  structures  has  a maximum  of  2n-l 
feedback  loops,  proven  as  follows.  In  every  structure,  local  feed- 
back loop  (i  to  i-1)  is  possible  for  each  section  without  inter- 
fering with  any  other  feedback  loop,  so  each  structure  has  these 
n local  loops.  Consider  the  structure  in  which  every  section 
except  the  first  in  Figure  4-2  (for  which  is  extraneous),  also 
feeds  back  to  the  plant  input  0,  giving  an  additional  n-1  loops, 
making  a total  of  2n-l  loops.  Obviously,  no  more  loops  are  pos- 
sible without  crossing. 

Let  us  hypothesize  that  the  maximum  number  of  loops  for  a 
plant  with  n-1  cascaded  sub-plants  is  independent  of  the  specific 
structure  used,  so  is  a constant  ■ 2(n-l)-l  and  that  this  is  true 
for  each  of  n*l,  2,  ...  up  to  n-1.  We  will  prove  that  it  is  true 
for  a plant  with  n cascaded  sub-sections.  Take  any  sub-section 
structure  whatsoever  and  give  it  its  initial  n local  loops. 

Starting  from  the  plant  input  0,  proceed  to  the  right  in  Figure  4-4, 
until  the  first  poirt,  say  n^,  is  reached  which  has  at  least  one 


loop  (other  than  from  n+1)  feeding  Into  point  n^.  There  will  lu* 
a loop  from  iij  to  0 because  there  Is  nothing  preventing  it.  Pro- 
ceed forward  (to  the  right)  from  n until  the  last  point  is  reached 
which  has  a feedback  loop  terminating  on  n^.  By  hypothesis  there 
must  be  at  least  one  such  point.  Let  this  be  point  n7>n^  (of 
course).  The  part  from  n^  to  n.;  can  be  considered  as  "semi-isolated" 
(n  -n  ) structure  with  a maximum  total  of  /l  loops,  because 

2 1 Vni 

there  can  be  no  loop  from  any  point  inside  reaching  outside  or  vice 
versa  without  crossing  the  loop  from  n,  to  n . 


n,  Ti, ♦!  nx  ns 


Figure  4-4.  No  crossing  multiple  loop  system. 


Now  work  backwards  from  the  last  point  n until  the  first 

point  n^  (Figure  4-4)  is  reached  which  has  at  least  one  loop  feeding 

back  to  some  n>0(n  <n  -1)  in  addition  to  0.  If  n,  has  more  than 
i 3 4 — ~ ■ 4 

one,  take  that  one  furthest  away  (to  the  left)  from  n^.  Note  that 
n jNn , by  hypothesis  because  no  crossing  is  allowed.  The  part  from 


n,  to  n,  also  forms  a "seml-lsolated"  multiple-loop  system  with 
5 4 
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L 


t 


Tln.-n.  maximum  loops.  Let  us  now  count  the  total  loops.  By  hypo- 
A 3 


thesis  ??  ” 2(n  -n_)-l,?7  ■ 2(n  -n, )-l.  The  number  of 

Vn3  4 3 n2'ni  2 1 


loops  from  the  first  n^  sub-plants  (from  0 to  n^)  is:  n^  local 


(i  to  i-1)  and  Oj-1  (2  to  0,  3 to  0,  ...»  n^  to  0),  giving  a total 


of  " 2n  -1.  The  number  of  loops  from  the  part  consisting  of  plant 


sections  from  n„  to  n_  is:  n -n  local  and  n -n„  + 1 to  0 (i.e., 
2 3 3 2 3 2 


from  n^,  n^+l,  .....  n^  each  to  0),  totalling  2(n^-;i2)+l . Finally 


the  number  from  the  sub-plants  from  n,  to  n is:  n-n,  local  and 

A A 


n-n.+l  (i.e.,  from  n.,  n -1,  .....  n to  0),  totalling  2(n-n,)+l. 
AAA  A 


Adding  all  of  these  gives  2n-l.  It  is  easy  to  prove  the  hypothesis 
is  true  for  n * 1,  2,  3.  Hence  it  is  true\/n. 

We  shall  concentrate  first  on  the  two  structures  shown  in 
Figures  A-1.2A,  with  considerable  detail  and  numerical  examples. 
From  this  there  will  emerge  a perspective  and  understanding  which 
will  enable  us  to  tackle  any  non-crossing  structure.  In  fact  it 
will  be  possible  to  know  in  advance,  without  a detailed  design, 
a great  deal  about  the  design  problem — just  as  in  the  cascade  no 
P.M.  design  approach  (recall  section  1.5). 


4.2  The  2-plant,  3-loop,  P.M.  system 
A. 2.1  Design  perspective  (Figure  A-l). 

The  important  design  factors  are  the  increase  in  signal 


level  of  and  the  effects  of  sensor  noise  sources  N^,  at 


tuJ 
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: 


: 


' 


and  X^.  As  In  chapter  3,  the  former  is  controlled  by  careful  de- 
sign of  (i.e.,  in  the  If  range),  although  the  actual  re- 

lation between  Q and  Q (of  Figure  3-23)  is  different  in  the  3-loop 

c 

case  (see  section  4.3.9).  Just  as  in  chapter  3,  can  handle  the 

uncertainty  problem  in  the  hf  range,  with  the  help  of  some  safety 

margin  in  the  IVJ  region  of  Figure  1-8  (5  db  was  used  in  chapter  3; 

recall  Figures  3-10,  16).  But  now  L * P2^2  ava^^a^^e  to  re“ 

move  the  burden  of  P uncertainty  from  L . Thus,  except  for  the 

2 o 

(5  db)  safety  margin  in  the  IVJ  region  of  Figure  1-8,  L need  handle 

o 

no  hf  uncertainty  at  all  and  only  P^e  If  uncertainty.  One  more 
point,  however,  this  safety  margin  can  be  divided  between  L and 
This  relative  strengths  of  N^,  and  the  dynamic  ranges  of 
the  stages  at  X^,  influence  the  way  the  division  should  be 
made.  Our  objective  is  to  give  the  designer  the  perspective  and 
understanding  of  the  trade-offs  to  enable  him  to  make  an  intelli- 
gent choice. 


4.2.2  Specifications  on  numerical  examples 
Plant:  P - P^P^  “ k^/s 

Plant  uncertainty:  /2£k^l0/2,  /Jjfk^lOv'T 

independently 

Bounds  on  | T( ja>)  | : shown  in  Figure  3-2(b)  were  origin- 

ally derived  from  time  domain  bounds 
of  Figure  3-2(a) 


"T T™-  ' 2 ''■""I 
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Disturbance  Response:  Y£2.0  db 

(All  above  are  the  same  as  In  sections  1.2.5,  3.3.1) 

System  structure:  as  shown  In  Figure  4-1, 

Restructlon  on  signal  level:  Q s 1.28  for  unit  step  Input 

(The  rational  function  values  of  the  compensation  functions 
in  Figure  4-1  are  given  in  Appendix  n) 


4.2.3  Single  and  cascaded,  2-loop,  no  P.M.  design 
Since  the  first  half  of  the  specifications  in  section  4.4.1 
are  the  same  as  in  sections  1.2.5,  3.3.1,  so  the  single  loop  and 
cascaded  2-loop,  no  P.M.  designs  of  sections  3.4.2,  3,  are  valid 
here. 


4.2.4  The  P.M.  inner  loop  L - P^H^ 


The  choice  of  |l+L^n|  (for  Q ■ 1.28)  is  made  in  accordance 


with  the  philosophy  of  chapter  3,  section  3.5,  resulting  in  Figure 


4-5.  The  inner  loop  bounds  B (oj)  in  the  If  range  and  the  re- 

Lll 


suiting  L,  are  shown  in  Figure  4-6. 
lln 


4.2.5  Design  of  L 


on 


Following  section  3.3.8,  templates  of  P = P^/(l+L^)  in 


R1 , R?  are  obtained  and  shown  in  Figure  4.7 # As  the  uncertainty 
> 

2 


r *2 

in  P will  be  handled  by  L 


^2  * **2^2*  on^  Pie  uncertalnty  nee<*  be 


rrencrr.  p 
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00- -1  w:.2  U)=  4 u)  = .63  to=  I (a)=1 

Figure  4-7.  Templates  of  P 


le 


handled  by  L . This  is  the  big  difference  between  the  3-loop  de- 
o 

sign  of  Figure  4-1  and  the  2-loop  design  of  Figure  3-1.  Bounds 

B,  (w)  due  to  P,  In  R, , R_,  are  shown  In  Figure  4-8.  With  5 db 
L le  1 l 

o 

hf  overdesign,  L is  designed  to  satisfy  R,  (to)  and  is  shown  in 
on  L 

o 

Figure  4-8.  All  above  procedures  are  the  same  as  In  2-loop,  P.M. 
design  in  chapter  3. 


0 u) 


Figure  4-8.  B.  (w)  and  L 
— L i 
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4.2.6  Determination  of  Bi  , , (u)  in  R , R,  , Rr  and 
**11 3 A 5 

completion  of  L 

in 

If  we  follow  section  3.3.8  precisely,  we  would  find  Bj  (u>) 

the  bounds  on  L such  that  L remains  satisfactory  despite  the 
il  on 

hf  uncertainty  in  P,  . This  is  done  in  Figures  4-9(a,  b).  The 
le 

next  step  would  be  design  of  followed  by  finding  R^  ^ ^ (to)  the 

bounds  on  L.„,  such  that  the  design  of  L remains  satisfactory 

12  o 

despite  the  uncertainty  in  P,.  Hut  in  such  a procedure,  we  do  not 

consider  L^.t  when  we  design  . This  increases  the  burden  on  L , 

especially  for  some  frequency  ranges. 

Figure  4-l0(a)  shows  L and  which  are  templates  of 

on  ” e 

L (ho)  at  several  co  values  due  to  P.  uncertainty  only,  with  P, 
o le  2 

fixed  at  P„  , i.e.,  templates  of  P,  (1  P,  . Tn  the  "hf"  range 
2n  2en  le 

hifiA),  ■>  is  the  forbidden  region  for  L ■ G P ,P,  . Suppose  the 

o le  2e 

template  of  P^  is  that  shown  in  Figure  4-10(b))  (it  is  so  because 
/ e • 


at  this  hf  range  L.„<0,  specifically  at  w • 18,  L,  „ 

i2  i2n 


-lOdb  -ISO  > 


At  u'  - 5,  the  template  of  L » G P.  P,  remains  satisfactory.  Bui 

o le  2e 

at  to  * 12,  it  is  definitely  not  satisfactory  because  it  violates 

the  hound  > « 2 db.  Thus,  for  toc[6,  20]  in  Figure  4-l0(a),  there 

is  not  much  margin  of  uncertainty  left  for  P . To  remedy  this, 

2e 

we  restrict  the  range  of  the  templates  of  G IV  P„  as  shown  in 

2en  2 

Figure  4-11,  instead  of  that  in  Figure  4-10(a).  This  makes  it 
harder,  of  course,  on  and  easier  on  1.  in  this  range.  Thus, 
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replacing  by  P , as  shown  in  Figure  4-14,  i.e.,  bounds  B^^^w) 

are  found  in  L.„  such  that  the  design  L remains  satisfactory 
i2n  on  7 

despite  the  uncertainty  in  P^  (and  P^).  The  results  are  shown  in 

Figures  4-15 (a),  (b)  and  L in  Figure  4-16.  Note  the  general 

12n 

nature  of  bounds  B^.,(u))and  are  exactly  the  same  as  In  the  cas- 
caded, no  P.M.,  inner  loop  B , L figure  1~19.  The  Bode  plots 

of  L , L , L.  of  3-loop,  P.M.  system  together  with  the  single 
on  11  12 

loop  system,  L , are  shown  in  Figure  4-17. 

9 

In  Figure  4-15(b)  the  level  of  Bl^  in  [5,  8]  range  is  so 
much  higher  than  for  BL^  (in  Figure  4-9(b))  because  the  5 db  margin 
in  this  range  (Figures  4-10,  11)  was  used  entirely  for  L . Hence, 
iLj^l  is  correspondingly  >|L  | (see  Figure  4-17)  in  this  range 

and  its  effect  persists  for  higher  frequencies.  This  was  done  here 
to  emphasize  the  trade-off  possibilities  between  and  L^* 

4.2.8  Time  response  of  the  3-loop,  P.M.  system 
CSMP  digital  simulation  of  the  3-loop,  P.M.  system  at  four 
extreme  parameter  cases  are  shown  in  Figure  4-18.  The  time  domain 
bounds  of  Figure  3-2(a)  are  included  as  dashed  lines.  It  is  seen 
the  results  satisfy  the  time  specifications. 


I 


2 


Time  <*) 


Figure  4-18.  Time  response  of  3-loop,  P.M.  system 
at  extreme  parameter  cases 

increasing  "bandwidth"  of  each  inner  loop.  [5,  13]  Compare  the 

bandwldths  of  and  Figure  4-17,  and  it  is  seen  this 

does  not  occur  here.  The  slight  increase  in  bandwidth  in  |L  | 

IZn 

is  primarily  due  to  the  5 db  margin  being  used  entirely  for  L . 

lln 

This  is  because  is  not  "inside"  as  is  the  case  in  the  non 

t 

P.M.  cascaded  system.  If  there  is  another  inner  loop  inside 
L , then  the  bandwidth  propagation  effect  will  occur  (see  section 


r 


• . 


h 


■ 
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4.2.10  The  signal  level  variation  ratio  (SLVR),  p 
Recall  In  section  2.2,  the  SLVR  p,  is  defined 


l°23 

max 

L 

o 

* 

L 

s 

F.  (1+L  ) 
le  o 

max 

1+L 

s 

max 

'C  3 

Is, 

max 

L 

s 

max 

L 

o 

iCs 

P. (1+L  ) 

1 s 

1+L 

o 

max 

3 max 


s m axf 


(4. 2-1. 2) 


The  frequency  ranges  are  redefined  here  (refer  to  section  2.3): 

R - [0.,  a)"] , |L  (juO  | - 25  db,  w'  - 0.08;  R'  - (ur,  or], 
x x sn  X X Z 1 Z 


■ 0 db  because  is  perfect  in  Lq,  cj'  • 1.6;  R^ 


w*]»  l^gn(j<i>j)|  " ~^2  " ~^0  db  with  is  the  hf  uncertainty 


°f  P2’  W3  " 8;  R4  " ^w3»  \ “ 10w'  - 80;  R'  ■*  (u)',  °») . 


We  wish  to  simplify  (1)  (see  section  2.4) 

|c„' 


p - a*  1+L 


3 max 


iln' 


s max 


Compare  (L  „ with  (L  )„  , „ 

on  3-loop,  P.M.  on  2- loop,  P.M. 


(4.2-3) 

in  Figures  4-19 


(a,  ....  c).  They  are  defined  in  the  same  way  in  R^  , R^  and  R'.R' 


giving  oi  £ 1 in  R^,  R^  as  in  sections  2.4.1,  2.  Also,  in  the  very 


hf  region  of  R;,  |lJ,  |l | , |L  |,  | L<?n|  <<  1,  so  the  result 


sn 


on 


",i2n 1 


a - 1 of  section  2.4.5  applies.  There  remain  R',  R*. 


4.2.10.1  R',  to  - (u^.  to'] , Figures  4-19(a),  (c) 


This  region  now  corresponds  to  the  latter  part  of  R^  in 


Figure  4-19(b).  Recall  in  2.4.2,  | L^/ ( 1-l-L^)  | is  likely  to  reach  its 


maximum  at  ^or  l^snl  * w^ere  M2  *S  t*ie  ^ uncertainty 


.I 


Figure  4-19 (a). 


L of  3-loop  P.M.  system  and  1. 
on  sn 


P2*  3-loop,  P.M.  case,  |L  /(1+L  )|  definitely  cannot 

reach  its  maximum  at  |P  I , so 

le  min 


L 

o 

c i 

L 

o 

P,  (1+L  ) 
le  o 

lp  | 

max  1 le'min 

l+L 

o 

(4.2-4) 


max 


And  the  first  two  factors  in  (2)  give 
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4o 
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ILIob 

0 

-20 

-40 

“60 

-360#  -270e  -(60*  -V  0° 

(b)  L of  2-loop,  P.M.  system  (Design  A)  and  L 

on  sn 

Figure  4-19.  Outer  loops  with  defined  regions 


Figure  4-19 (c).  Bode  plots  of  loop  transmissions 
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P.  (1+L  ) 

1 le  o 'max 


\ V1+Ls> 

* 1 

max  J 

L 

o 

\l 

\ 

^le|min 

1+L 

o 

max  • 

1 

L 

s 

JPl|min 

1+L 

s 

max 

1 min 


| le | min 
i 1 1+L.  I 


I in 

So  (2)  becomes 


(4.2-5) 


p - a |l+Llln| 


3 max 


(4.2-6) 


a £ 1,  for  a)  in  0 > lLgn0w)|  >■  -M  . 

If  ^^2’  c^en  covers  the  whole  range  of  R^.  If  M^M^,  then 

|L  /(1+L  )|  may  not  be  achieved  by  |P, I . in  the  range 
I s s ] max  | 1 |min 

(a)  |-  M_>Il  (ju))l*  -M  ).  But  now  |l  I is  fairly  small  (Figure 
2 an  1 on 

4-19(a)),  so  that  at  I P_  I . * |L  /(1+L  )|  tends  to  be  further  from 

le  min  o o 

its  maximum  than  its  counterpart  in  the  single  loop  system.  So  (6) 
still  hold 


p - a | 1+L 


■'31  max 


iln  |C 


(4.2-7) 


a £ 1 in  R, 
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The  result  of  the  numerical  example  of  a 3-loop,  P.M. 
design  Is  In  Figure  4-20(a),  and  agrees  with  (7). 


4.2.10.2 


. |Lgn(jwJ  I— M^b,  |Lsn(jo)')— (M^Jdb. 
(Figures  4-19 (a),  (c)) 


Compare  (L  ) , (L  ) In  Figures  4-19 (a),  (b). 

on  2-loop  on  3-loop 

in  the  3-loop,  P.M.  system  corresponds  to  the  latter  part  of 

R,  of  the  2-loop,  P.M.  system.  The  nature  of  I L I In  both  cases 

Is  quite  similar,  because  |L  I is  fairly  small.  So  the  results  for 

1 sn' 

the  latter  part  of  R,  in  the  2-loop  system  of  section  2.4.4, 

4 

applies,  i.e.. 


P 


a*  1+L 


ISL, 


iln1  |C 


s max 


(4.2-8) 


a £ 1 in  R. 

4 

which  agrees  with  the  result  of  the  numerical  example  it  Figure 
4-20(a) . 

In  conclusion,  as  ot£l  for  all  U),  design  based  in  ll+L^I 

is  definitely  conservative  even  more  so  than  in  the  2-loop  design. 

Comparison  of  I 1+L . , I with  p/ [ I C_ I / I C I ] is  shown  in 
1 iln1  1 3 'max  1 s' max 

Figure  4-20(b),  with  Q - 1.26,  which  satisfies  the  specification 
Q<1.28. 

Comparison  of  a in  the  2-loop  and  3-loop  P.M.  systems  is 
shown  in  Figure  4-21.  In  the  2-loop  P.M.  system  (D) , Q ” 1.24  and 


in  the  3-loop  P.M.  system,  Q ■ 1.26. 


4.2.11  Sensor  noise  effect 


We  repeat  the  noise  response  function  of  section  1.6.3. 
L 


T ^ — 

N1  N1 


__c^  T1 

PlV  Nl“Nl 


(L,,+L  ) 

il  o 


2 A X2  Li2  1 A X1 

T ■ — * - , T - — » -L 

N2  N2  P2  N2  N2  i2 


(4.2.9-12) 


at  hf  where  L 


I ^ii I * |l12|«i. 


In  the  hf  range,  all  parameter  uncertainties  in  P^,  P0  are 


handled  by  the  inner  loops  of  L..,  L.^,  so  the  bandwidth  of  L is 

il  12  on 


fairly  small  (see  Figure  4-17).  Hence,  the  noise  effect  of  T in 


! 
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(9)  Is  considerably  better  than  in  the  2-loop  P.M.  design.  The 

12  1 

other  noise  effects  T , T , T in  (10,..., 12)  are  also  rela- 

N,  N„  N„ 

1^2 

tively  much  smaller  than  T . The  reason  is  that  the  "hf  asymptotic 

N1 

region,  (where  |L  | begins  its  final  downward  slope — point  X in 
on 

Figure  4-17)  occurs  at  a much  lower  frequency  range  (uT*8)  in  the 
3-loop  design,  than  in  the  2-loop  design  (u)'=55)  in  Figure  4-19(c). 
The  peaks  of  |L  [ and|L^2nl  occur  approximately  at  this  point  ( 
see  Figure  4-19(c)).  This  is  inherent  in  cascaded  design.  The 
results  are  compared  in  Figure  4-22  with  the  single  loop  and  cas- 
caded, 2-loop,  no  P.M.  systems;  in  Figure  4-23  they  are  compared 
with  the  single  loop  and  2-loop,  P.M.  (Design  D)  systems.  Note, 
in  Figure  2-23,  Design  D has  Q » 1.24  and  the  3-loop,  P.M.  system 
has  Q ■ 1.26,  so  the  signal  level  increases  are  quite  similar, 
but  the  noise  improvement  is  much  bigger  in  the  3-loop  P.M.  than 
in  the  2-loop  P.M.  system  for  the  reason  given  above. 

4.3  The  3-plant,  5-loop,  P.M.  system 
4.3-1  Design  philosophy 

Design  in  Figure  4-24  commences  exactly  as  in  the  3-loop 

P.M.  system  (section  4.2.4).  The  bounds  on  (i.e.,  on  in 

the  If  range)  are  obtained  to  control  the  signal  level  increase  at 

. Then  L f*.  is  obtained  in  the  If  range,  giving  P,  (jio)  in  the 
2d  11  le 

If  region.  The  bounds  on  the  outer  loop  L are  obtained  to  handle 

on 

the  uncertainty  in  P in  If  only  plus  the  safety  margin  (5  db  in 

le 


200 


(b)  In  Ar i thematic a 1 scale 


Figure  4-22. 


Comparison  of  noise  effects 


mmwL  k L 


Pi 


D* 


p, 
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p2e  “ P3e  = V^V*  Lo2  “ G2P2eP3e* 

T2  " P2eP3e/(1+Lo2>*  Lo  “ GlVle 
Figure  4-24.  3-plant,  5-loop,  P.M.  system 


Figures  4-8,  11)  in  the  A region.  In  the  3-loop  system  the  next 

step  was  to  find  the  bounds  on  L ^ to  handle  the  uncertainty  In  P,,. 

But  here  we  have  the  section  in  Figure  4-25(a)  to  replace  P^H^. 

Our  next  step  is  to  assign  some  safety  margin  (say  5 db)  V w to  the 

outer  loop  L . ■ G„P„  P_  of  T„,  i.e.,  as  if  the  uncertainty  in 

P„  P„  is  5 db.  This  gives  L „ . The  next  step  is  to  find  the 
2e  3e  6 o2n 

bounds  on  L.  „ to  handle  the  actual  uncertainty  in  P„,  i.e.,  to 
i2n  2 

ensure  that  L is  satisfactory  despite  the  uncertainty  in  P , (at 
oZn  z 

P^  ■ P fixed).  (We  should  also  worry  at  this  point  about  the 

effect  of  on  the  signal  level  at  C^,  but  it  will  be  shown  in 

section  4.3.5,  that  there  is  no  need  to  worry.)  The  final  step  is 

to  design  L to  handle  the  uncertainty  in  P . 

13n  3 


j 
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4.3.2  Specifications  on  numerical  examples 


(Figures 4-24.  ....  27J 


Plant : 


Plant  uncertainty: 


Bounds  on  |T(ju))| 


p - Wj.  \ - V 

1.71*k,S17.1,  1. 7Uk  <13.58, 

1.71$k^30.41  Independently 

shown  in  Figure  3-2 (b)  were  origin- 


ally derived  from  time  domain  bounds 
of  Figure  3-2 (a). 


Disturbance  Response:  y£2.0  db 


System  structure: 


Restriction  on  signal 


Figure  4-24. 

<4  1 2(103 

level:  V/TETo^*1-28’ 

o ' 2s  1 

a ric33(o,)* 2<io> 

Q * — — £ 1.30  for  unit 

/°°|C35(a))|  du 


input.  Note  the  need  for  > Q^. 
(All  above  specifications  are  the  same  as  in  section  4.2.2 
except  for  plant  and  system  structure.) 

(The  rational  function  values  of  the  compensation  functions 
in  Figure  4-24  are  given  in  Appendix  II.) 


4.3.3  Single  loop  and  cascaded.  3-loop,  no  P.M.  systems 


The  single  loop  system  (Figure  4-26)  and  cascaded,  3-loop, 
no  P.M.  system  (Figure  4-27)  are  also  designed  here  in  order  to 


2 


compare  with  the  5-loop,  P.M.  system  in  SLVR,  Q and  sensor  noise 
effects.  The  Bode  plots  are  shown  in  Figure  4-28 


4.3.4  First  P.M.  inner  loop  L._  = PH  and  outer  loop  L 


Ljh  = 5s  fjn 


A,LA 

IL..I  \ 


Figure  4-28.  Bode  plots  of  | L | , | L | , | L | , | L 


1 ™ 
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4.3.5  Second  outer  loop  L - G.P.  P_  , second  P.M. 
oz 2 ze  3e 

A 

Inner  loop  L ^ ■ P2H2  an<*  C^e  t*1*r<*  inner  loop 

L,_  - P_H 
13  3 3 

The  dash  box  (T^)  In  Figure  4-25(a)  Is  considered  next. 

At  this  stage  the  problem  is  to  control  the  variation  such  that 

the  designed  L and  L.,  are  satisfactory  despite  the  fact  that  P , 
o iln  2 

P^  uncertainties  were  neglected  in  their  design.  Three  loops  are 
available  for  this  purpose:  Lo2’  Li2*  Li3‘ 


4. 3.5.1  Design  of  (Figure  4-29) 


In  the  3-loop  system  of  sections  4. 2.4-7  and  Figure  4-1,  the 

outer  loop  L had  to  handle  the  P,  If  uncertainty  (analogous  to 
o le 

P here)  and  P,  If  design  was  dominated  by  consideration  of  signal 
2e  le 

level  at  (^35  Figure  4-24).  So,  presumably  similar  pre- 

cautions must  be  taken  here.  Fortunately  this  is  not  so.  It  is  re- 
called from  chapter  3,  section  3.3.7  (Figures  3-9,  10)  that 
design  shifted  the  template  which  Lq  had  to  handle,  to  the  left 

where  it  could  cross  B.  , facing  L to  be  more  conservative.  But 

h2 

here  ABC  in  Figure  4-29 (a)  replace  B.  , so  excursions  of  the  tem- 

h2 

plate  of  the  new  P,  (i.e.,  P„  of  Figure  4-24)  to  the  left  are 
le  Ze 

acceptable. 

The  second  point  is  that  in  the  3-loop  design  could  not 
handle  entirely  the  If  uncertainty  of  P^  (in  Figure  4-1)  because  of 
the  signal  level  problem.  But  here  because  of  the  'bandwidth 


T"r' 
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propagation  effect'  discussed  in  section  4.3.7,  we  can  neglect  the 

signal  level  problem  of  the  new  P,  (i.e.,  P_  of  Figure  4-24). 

le  2e 

Thus,  the  frequency  corresponding  to  of  Figure  3-24  is  much 

larger  (50  in  Figure  4-30(c)  cf  6 in  Figure  4-13).  Hence,  there 

is  very  little  obligation  on  the  new  outer  loop  (i.e.,  L „ of 

oZ 

Figure  4-24).  However,  as  a safety  factor,  we  design  L as  if  it 

o2 

has  to  handle  a few  db  uncertainty  over  all  u)  (somewhat  analagous 

to  the  5 db  safety  margin  used  in  Figure  4-8  for  L there).  In 

o 

this  design  example  5 db  was  used. 

A third  factor  is  the  analog  of  the  trade-off  (Figures  4-10, 
11)  between  and  in  the  3-loop  P.M.  design.  Here  the  trade- 
off is  between  L and  L , and  can  be  made  in  accordance  with  the 
11  02 

levels  of  N^,  N^.  In  this  specific  example,  the  5 db  region  in  A 

was  given  to  and  the  region  in  B was  split  roughly  equally 

between  L.,  and  L .. 

11  o2 

4. 3. 5. 2 Design  of  (Figure  4-30) 

Special  design  of  in  the  If  range  in  order  to  control 
SLVR  at  C^»  is  not  necessary  because  of  the  bandwidth  propagation 
effect  discussed  in  section  4. 3. 5.1.  is  designed  simply  to 

handle  the  uncertainty  of  P , i.e.,  bounds  B in  Figures  4-30(a, 

a Lla 

b)  are  obtained  such  that  L „ is  satisfactory  for  P uncertainty 

o2n  2 

(18  db  here)  despite  the  fact  that  L „ was  designed  on  the  basis 

ozn 

of  5 db  uncertainty.  L is  then  obtained  to  satisfy  B , as  shown 

i 2n  112 


in  Figure  4-30(c) 


Lilh/.oS 


~ ^ 
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4. 3. 5. 3 Design  of  L (Figure  4-31) 

The  design  of  L is  exactly  the  same  as  the  second  inner 
loop  of  a 3-loop,  P.M.  system  (L  of  Figure  4-1).  Bounds  B are 

X Z LI  J 

obtained  such  that  L , L are  satisfactory  despite  the  actual 

o2n  12n 

uncertainty  of  P , from  which  L, , is  as  then  found.  In  Figure  4-31 
3 IJn 

(b),  thr  level  of  B g(w)  i°r  U)£l20.,  75.)  is  much  higher  than  the 
level  of  B ^(w)  *n  F*8ure  4- 30(b).  This  is  because  we  used  all  the 
safety  margin  of  L „ for  I.  (recall  Figure  4-10,  11  in  section  4.2.6). 
Hence  lLi:jnl  is  correspondingly  >>lLi2n^  ^see  Figure  ^-32)  in  this 
range  and  the  effect  presents  for  higher  frequencies.  We  thereby 
emphasize  the  possibilities  for  trade-off  between  and  which 

is  useful  in  the  case  of  significant  differences  in  the  levels  of 
N^,  in  Figure  4-24. 


4.3.6  Time  Responses 

Eight  extreme  parameter  cases  were  simulated  on  the  digital 
computer  (CSMP)  and  shown  in  Figure  4-33,  including  the  time  domain 
bounds  of  Figure  3-2 (a) 

4.3.7  The  bandwidth  propagation  effect 

In  Figure  1-20,  in  the  single  loop  design,  the  loop  trans- 
mission L can  take  its  final  hf  slope  only  ~ at  point  E.  The 
relatively  horizontal  segment  XE  is  a preparation  for  the  final 


breakpoint.  In  the  2-loop  design,  E is  replaced  by  m 

b 


CO  *0 


'•I*!' ..-'I'  -t-V 
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for  the  outer  loop  which  Is  significantly  to  the  left  of  E.  How- 
ever, the  Inner  loop  L cannot  assume  Its  final  hf  slope  until 

2n 

04^2  which  is  slightly  to  the  right  of  E.  This  is  the  'bandwidth 
propagation'  effect.  The  critical  region  for  the  inner  loop 
(Wbi*  ui^)  *8  •hifted  to  a range  significantly  higher  than  that 
for  the  outer  loop  (~  X^,  for  the  latter).  This  property  is 
also  noted  in  Figure  3-28.  In  the  general  n-cascade  no  P.M.  design 
where  there  are  n nested  loops,  the  interval  between  points  and 
X in  Figure  1-20  would  be  split  up  between  the  inner  loops.  A 
3-loop  no  P.M.  case  is  shown  in  Figure  4-28. 

In  the  P.M.  case  it  is  possible  for  several  inner  loops  to 
have  the  same  "critical  bandwidth"  region  in  the  above  sense,  be- 
cause they  are  not  nested  but  working,  so-to-speak,  in  parallel. 
This  is  so  for  (Figure  4-17)  in  the  3-loop  design.  (It 

is  recalled  that  the  difference  there  is  due  to  the  unfair  division 
of  the  5 db  safety  margin).  Similarly,  in  Figure  4-32,  and  L 
are  in  "parallel"  and  L^*  are  in  parallel  with  the  reason  for 
the  differences  explained  in  4.3.8. 


4.3.8  Trade-offs  between  loops 

In  the  3-loop  design  of  section  4.2.6,  let  us  distinguish 
between  regions  A and  B in  Figure  4-lo(a).  The  safety  margin  in  A 
was  given  all  to  while  that  in  B was  split  between  and 


The  result  was  significantly  larger  level  of  L.„  than  L , 

i2n  iln 

(Figure  4-17),  in  the  <a  range  corresponding  to  A (Figure  4-10(a)) 
which  persisted  in  its  effect  to  higher  u>  (because  mu8t 

then  decrease  over  a larger  interval) . This  was  repeated  pre- 
cisely in  the  5-loop  design  of  section  4.3.4,  between  L , and 

iln 

L . (which  effectively  replaces  L. „ in  this  matter),  leading 
o2n  iln 

to  the  significantly  larger  | L J (Figure  4-32)  than  |L  |. 

o2n  Iln 

The  above  problem  was  repeated  in  the  design  of  L. , L,- 

12n  13n.  a2n 

is  shown  Figure  4-29 (c).  Because  of  the  'bandwidth  propagation' 

effect  A'  is  now  at  a higher  u>  range  than  A.  Again  some  extra 

margin  is  assigned  to  L , in  A'  (just  as  such  was  assigned  to  L 

o2n  o 

in  A).  And  again  we  let  L^n  have  all  of  it  in  A',  while  that  in 

B'  was  split  between  L,_  and  . The  result  is  L,_  level 

12n  13n  i3n 

significantly  larger  than  L.*  level  in  Figure  4-32  and  which  per- 

IZn 

slsts  as  before.  The  difference  between  L. _ and  LJO  in  Figure 

lJn  12n 

4-32  is  > between  L_,  and  l.„  In  Figure  4-17,  because  in  the 
iln  i2n 

former  the  uncertainty  of  > uncertainty  of  (25  db  cf  18  db), 
whereas  in  Figure  4-17,  the  uncertainty  of  P^,  P^  are  equal 
(20  db  each). 

As  previously  noted,  the  above  can  be  used  as  trade-off 
between  the  loops,  to  compensate  for  significant  differences 


between  sensor  noise  levels 


'"  I-  IW'.'W’ 
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A. 3.9  The  signal  level  variation  ratio  (SLVR),  p 
The  definition  of  SLVR  is  extended  to  account  for  the  addi- 
tional plant  modification  In  the  5-loop  case,  over  and  above  that 
in  the  3-loop  case.  (Figures  4-24,  26). 


|C 


25  max 


P3  ' 1C. 


'33 'ms 


'2s'max  ” 1 3s 'max 

Recall  section  4.2.10,  it  was  proven 


(4.3-1) 


P„  " a *|l+L 


5 'max 


2 2 ' iln'  C 


s max 


(4.3-2) 


The  frequency  ranges  for  our  numerical  example  are  (section  4.2.10): 

Rj  - [0,  up,  lLgn(Jwpl  * 25  db,  ux  - 0.08;  R*  - (u>',  u)'J , 

- 0 db,  it'  - 1.6;  R'  - (o>',  u'] , ^(jup  - -(M2+M3)- 

-43  db,  with  H^,  the  hf  uncertainties  of  P2>  P3  respectively, 

it)'  » 55;  RT  ■ (w«,  it)'] , it)'  » 10  it)'  • 550;  RjT  ■ ( id °°).  The  results 
3 434**  3 54 

of  the  5- loop  numerical  example  Is  shown  In  Figure  4-34 (a),  and 
agrees  with  (2). 

The  SLVR  for  C is 


P 


3 


lC35^max  J C35  ^ max/ ^ max 
^C3s^max  ^ C3s ^max/ ^ C2s^max 


(4.3-3) 


Q32*  ^ 1+L12n^ 


(4.3-4) 


v ;^-  m u mm 
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with  a^2  £ It  because  the  dashed  part  T 2 in  Figure  2-25 (a)  acts  as 


another  3-loop,  P.M.  system.  Combining  (2,  4) 


P3  ' a3*l1+Lilnl*l1+Li2nl 


(|CJ  / | C | ) 

5 max  1 s max 


(4.3-5) 


“3  " “23-a2  S I 


The  results  of  the  example  are  shown  in  Figure  4-34 (b)  and  agree 
with  (5).  The  nature  of  U+I^J.  VOc^/lcJ^). 

,1+Lilnl'|l+Li2nl  a“d  P3/(|C5l»ax/,Cslmax)  *re  shom  10  F1*“res  4'35 
(a) , (b) , respectively. 


4.3.10  Sensor  noise  effects 


The  sensor  noise  functions  of  the  single  loop  and  the  cas- 
caded, no  P.M.  system  for  a 3-section  plant  (Figures  4-26,  27)  are: 
(1)  single-loop  system 

™ 3 A 


N. 


X3 

L 

s 

, T2  £ 

X2 

L 

s 

N1 

P P„P, 
12  3 

V 

N1 

P1P2 

T1  - 
*1 


x. 

L 

_1 

8 

N, 

P, 

1 

1 

, all  at  hf  where  L «1. 

1 s1 


(4.3-6) 


(2)  cascaded,  no  P.M.  system 
_ 3 A 


N. 


X3 

L1 

T2  £ 

*2 

L1 

N1 

P1P2P3 

»1 

P1P2 

T1  £ 
N1 


X. 

L„ 

1 

_1 

N, 

P, 

1 

1 

all  at  hf  where  all  L «1. 


(4.3-7) 


N1  N1 


3 A ^3 
n2‘  n2 


1 A \ 

n2"  n2 


3 A X3 
r ■ — 
n3  n3 


1 A X1 
T ■ — 

n3  n3 


all  at  hf  where  all  |l|«1.  (4.3-10) 


X 

L 

T2  £ 

_2 

12 

’ Y 

N„ 

P„ 

2 

2 

2 

I ^*12 ^ hf  where  all  |l|«1.  (4.3-11) 


L 

X 

13 

T2  £ 

_2 

P„ 

’ Y 

N 

3 

3 

3 

YYl  a11  at  hf  where  a11  |L|«1.  (4.3-12) 


Results  for  the  design  example  together  with  the  single  loop  system 
are  in  Figures  4-37 (a,  b).  Comparison  of  5-loop  P.M.  system  and 
3- loop  (no  P.M.)  system  are  in  Figures  4-38  which  shows  the  big 
saving  due  to  the  5-loop  system. 


N2  5-loop,P.M.8ystem 


loop  sys 


(a)  In  Log  scale 


0 


1000 


2000 


3000  oJ 


(b)  In  Arithematical  scale 

Figure  4-37.  Noise  response  of  the  single  loop  system 
and  the  5-loop,  P.M.  system 


Figure  4~38.  Comparison  of  noise  effect  between 


4.3.11  Signal  level  at  plant  input 


In  considering  signal  level  increases,  we  have  always 


gone  backwards  up  to  and  including  the  output  of  the  first  plant 


section  (X  . in  Figure  39) , but  never  bothered  with  input  to 


The  reason  is  that  all  feedback  to  X can  proceed 


Figure  4-39  Cascaded,  no  P.M.  multiple  loop  system 


via  a compensation  network  in  cascade  with  the  plant,  e.g.  G 


in  Figure  4-39,  rather  than  via  P . The  dimensions  at  Q 


can 


be  electrical,  not  the  physical  ones  of  the  plant  (which  may  be 


chemical,  aeronautical  etc 


Thus  X 


X ,/P  and  in  this  way  need  not  be  affected  by  the  feedback 


paths  to  X , which  in  this  work  we  generally  showed  in  the  re 


turn  path 


4.3.12  Design  philosophy  and  perspective  for  second  type 


The  second  type  5-loop  PM  structure  is  shown  in  Figure 


4-40(b)(the  first  was  Figure  4-24  and  replotted  in  Figure  4-40(a)) 


The  outer  loop  L - G L P is  designed  first  with  L , - G„P0  P, 
o ± ui  je  Oi  2.  2e  ie 

<ple  ■ + Ln)»  P?p  ■ P?/(l  + Ll2)»  Ltl  “ P1H1*  Li2  “ P2H2^ 


and  ■ P3^3*  *s  assigned  to  handle  a few  db  (say  8db) 


>TT« 


3eP2eG2*  T2  “ P3eP2e/(1  + Lo2) ’ L 
(a)  Structure  I 


" P2ePleG2*  T1  “ Vie7  <1+Lol>  ’ Lo  " P3eTl 
(b)  Structure  II 

' Vl'  L12  ‘ W Li3  ’ P3H3‘  Plo  ' V<1+Lll> 

" p2/(1+l12>*  P3e  “ V<1+Li3> 


Figure  4-40  5-loop.  P.M.  systems 
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uncertainty  as  a safety  margin  for  the  next  step.  Suppose  we 
take  the  specific  problem  in  Figure  4-40(b)  where  the  hf  un- 
certainty is  63  db,  and  a single- loop  design  gives  Lgn  in  Figure 

4-41.  Measure  63-8  ■ 55  db  up  from  the  plateau  of  L at  A 

sn 

in  Figure  4-41  to  give  LQn  with  new  plateau  at  A'  the  nominal 

outer  loop  of  our  5- loop  design. 

The  loop  L . = P.  P,  C,  and  L._  ■ P_H,  now  constitute 
ol  le  2e  2 13  33 

two  "parallel"  inner  loops  (the  first  inner  loops)  next  designed. 
Normally  we  would  design  in  the  same  manner  as  was 
designed  in  Figure  4-13,  i.e.  SLVR  considerations  dominate. 
However,  here  Lq1  need  not  handle  the  entire  P.^  uncertainty 
(whereas  in  Figure  4-13  did  have  to  handle  the  P^  uncertainty 
entirely),  because  here  loops  P^H^,  P^  are  available. 


Figure  4-41  The  perspective  design  of  Structure  II 


tii 
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SLVR  considerations  due  to  the  other  loops  is  considered  later. 

We  could  assign  some  few  db  safety  margin  (say  8 db)  for  L , 

ol 

to  handle.  We  now  recall  Figure  3-24  where  we  considered  use 
of  "positive  feedback",  locus  A’B'C'MNQ  for  L^,  which  is  much 
better  than  A' ' 'B' ' 'C* ' 'MNQ  in  terms  of  SLVR.  There  it  was  not 
useable  because  L^  had  to  handle  bf  uncertainty  of  'v  20  db  and 
a positive  feedback  path  could  not  achieve  it  (recall  the  dis- 
cussion in  Section  3.5.4).  However,  here  need  handle  only 
a token  small  uncertainty  factor,  so  the  "positive  feedback" 
path  is  feasible.  Hence  we  can  pretty  well  ignore  the  SLVR 
problem  in  design  of  Lq^. 

We  can  make  a good  estimate  of  L , without  a detailed 

oln 

design.  It  is  an  inner  loop  following  Lo>  so  its  peak  will  be 

in  the  u region  (A')  of  the  plateau  of  L , as  shown  in  Figure 

4-41.  Then  |LQ^n|  will  decrease  "slowly"  (<\<  average  of 

25-30  db/decade)  until  the  8 db  margin  has  been  achieved,  followed 

by  a plateau  (at  X in  Figure  41)  and  then  rapid  decrease. 

L, , is  next  designed  to  handle  the  uncertainty  in  P„. 
un  j 

Since  it  is  in  the  first  inner  loop  group,  "paralleling"  L 
its  peak  is  also  ^ A'  location,  but  after  this  peak  it  must 
decrease  "slowly"  until  the  25  db  uncertainty  of  P^  has  been 
handled.  Thus,  its  plateau  at  CD  is  at  higher  u>  than  that  of  L^^ . 
Of  course,  just  as  in  Sec.  4.2.6,  there  is  the  choice  of  how  to 
split  up  between  and  Lt3n  t*ie  reserve  provided  by  L^. 


fl'H 
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We  next  turn  to  the  two  "parallel"  second  Inner  loops 

L11  “ P1G1’  Li2  " P2G2'  Llln  is  designed  to  handle  the 

uncertainty  of  Pj,  and  L3n  that  of  P^;  and  the  reserve  due  to 

LQln  design  can  be  spilt  up  as  desired.  In  Fig.  4-41,  note  that 

the  peaks  of  L^  , L12r  occur  in  the  id  range  (-v  X)  of  the  plateau 

of  inasmuch  as  this  parallel  inner-loop  pair  branch  off 

from  L0^n»  ^02n  were  drawn  as  if  there  was  even  division 

of  the  L , reserve, 
oln 

Note  how  the  inner  loops  can  be  drawn  approximately, 

without  the  need  for  a detailed  design,  - at  most  only  of  L , 

sn 

the  single  - loop  transmission.  We  should  next  consider  the 
relative  merits  of  Structures  I (Figure  4-40(a))  and  II  (Fig- 
ure 4-40(b)),  but  will  first  consider  more  fully  the  SLVR 
situation  in  Figure  4-40(b). 


4.3.13  SLVR  in  structures  I,  II  (Figures  4-40(a,b) 

It  is  readily  found  that 

C35/C5  - [(1+Lo1)(1+L11)(1+L12)]/(P1P2)  in  II 

C35/C5  - [(1+L11)(1+L12)]/(P1P2)  In  I 

The  effect  of  (1+L^)  has  been  previously  mentioned  as  not  being 
dominant,  so  it  would  appear  that  the  SLVR  problem  with  respect 
to  C33  is  basically  the  same  in  I,  II.  There  is,  however,  a 
difference.  In  I,  L^  is  in  the  first  inner  loop  pair,  so  the 


|l  + L^|>1  factor  is  effective  in  a lower  frequency  range  than 
in  II  where  is  in  the  secpnd  inner  loop  pair  (bandwidth 
propagation  effect).  Hence  SLVR  would  be  bigger  in  Case  I, 
because  is  in  the  second  inner  pair  for  both  I and  II. 
Consider  “ (1+1^)  ^f  ^or  k°th  As  L is  in  second 

inner  loop  pair  in  Case  II  (in  first  for  Case  I),  SLVR  tends  to 
be  bigger  for  Case  I. 

4 . 3 . 14  Noise  response  comparisons  between  structures  1,11 

The  noise  responses  of  Structures  I, II  are  compared 

f\j 

in  Table  4-1.  For  example,  at  hf  X2/N^  'v  f°r  both.  However, 

in  II,  is  in  the  first  inner  loop,  so  its  final  asymptotic 

descent  begins  earlier  than  in  I.  Hence,  in  the  hf  range 

^Li3^I  > lLi3^II*  Tlle  si8nificance  of  a,b,c,d,e,g  depend,  of 
course,  on  the  respective  L/P  ratios.  If  ] P J is  large,  then  the 
effects  are  unimportant  but  in  large  uncertainties  it  is  more 
likely  that  the  L^/Pj  are  large.  If  the  level  at  is  the 
major  factor,  then  item  (b)  (as  far  as  sensor  noise  is  concerned) 

is  the  one  to  examine,  but  of  course,  the  total  RMS  level  due 

to  all  sources  must  be  considered,  so  even  if  the  difference 
is  very  large  in  (b) , it  does  not  matter  much  if  (b)  is  still 
small  compared  to  (a).  If  the  X^  level  is  important,  then 
(d)  should  be  contrasted  with  (f). 


Table  4-1  Noise  responses  (in  hf  range  where  | l-*-L  [ ^1) 


high  frequence 
low  frequence 
Inner  loop  pair 
linear  time  invariant 
linear  time  varying 
left  hand  side 
right  hand  side 

Ignoring  signal  level  variation 
noise  power  ratio 
plant  modification 
signal  level 

signal  level  variation  ratio 


root  mean  square 


[j 

I 

•• 
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APPENDIX  I 

Patterns  of  Templates 

Seven  typical  patterns  of  templates  P^  and  P2  are  con- 
sidered. They  can  be  divided  into  three  categories.  The  first 
category  extends  upward  and  to  the  right  from  the  nominal  point 


in 


(see  Section  2.2.1),  as  shown  in  Figure  I-l(a).  And  the 


i Pi 


i 


i 


, p*  i 


i 1 

i l 

i p.  i 

i i i 


in 


in 


(a) 


(b) 


Pm 

(c) 


Figure  1-1.  Three  categories  of  template. 

second  category  extends  upward  and  to  the  left  from  P^  , as 
shown  in  Figure  I-l(b).  The  third  one  extends  both  the  right 
and  to  the  left  from  P^n  , as  shown  in  Figure  I-l(c).  The 
following  three  plant  functions  give  templates  of  the  shape  of 
Figure  1-1 (a) . 
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(1)  P (s) 


Pa(s)  - 


s7s7p+I)  ' k " [kl'k2]  ' p * [pl'p2] 


n sts/pj+l) 


(1-1) 

(1-2) 


(2)  Pu(s) 


k(s/z  •*•  1) 


s (s/p  +1) (s/p  +1) 

a d 


k - [k^] 


Figure  1-2.  Typical  templates  of  Pa(s) , Case  (1). 

in  Case  (3)  the  pole-pair  ranges  inside  MNQU  and  M’N'Q'U'  in 

S-plane  of  Figure  1-4 (a).  Three  templates.  Figures 

1-4  (b) , . . . , (d) , are  plotted  with  u * u . If  u>  < u> 

o.  « 


then 


the  template  extends  to  the  left  which  puts  it  in  the  second 
category. 

The  next  three  typical  structures  have  the  shape  of 
Figure  I- 1(b). 


(4)  Pu(s) 


s(s+a)  ' k " [kl'k2]  ' a “ ^I'V 


(1-7) 


Pn^  s(s+a_) 


(5)  Pe(s) 


2 

*lg±z.> k = Tic  k 1 

s(s+pa)  (s+p^  ' * lVk2J 


(1-8) 


Z - t«r*2]  * pa  * IPai'Pa2]  ' pb  “ lPbl'Pb2] 


(1-9) 


. Jc.(s+z.) 

e . . l i 

nts  = s(s+p  ) (s+p.  ) 


(1-10) 


p;  (<“c  > 


Q 

e 1-4  (a) . P parameter  uncertainty  in  s-plane 

(b) , . . . , (d) . Typical  templates  of  PC(s)  , Case  3 
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(6)  P (s)  - s(s+0+jWo) (s+a-jwo)  ' k " tki'k2]  ' 


o  - [o  ,o_]  , 0)  - [to  ,o>  ] 

1 2 o °i  °2 


(1-11) 


f i 

P (s)  « - 

n s(s+o  +ju>  ) (s+o_-ju>  ) 

2  2 


(1-12) 


Similarly,  the  respective  shapes  of  templates  at  three  different 
frequencies  are  shown  in  Figures  1-5,..., 7. 


IpW 


paK> 


p>a>  P>e> 


Figure  1-5.  Typical  templates  of  P (s)  , Case  (4). 


Note,  incase  (6),  the  complex  pole-pair  ranges  within  MNQU  and  M'N'Q’U' 

in  the  S-plane  of  Figure  1-4 (a).  Three  frequencies,  u>A  < < u»c 

in  Figure  1-7  are  shown  in  Figure  1-4 (a).  Again,  if  w < w , 

°1 

the  template  extends  to  the  right  which  puts  it  in  the  first 
category. 


The  third  category  is  the  template  extending  both  to  the 


Figure  1-6.  Typical  templates  of  P (s)  , Case  (5) 


Figure  I- 


. Typical  templates  of  Pf(s)  , Case  (6)  with 
parameter  uncertainty  in  Figure  2-14 (a). 


TS™-: S' "WSTT-  . 


mmmm — - • : 
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right  and  left  of  Pn(jui)  . Such  templates  definitely  occur  when 
the  plant  set  includes  both  minimum  and  nonminimum-phase  elements,  or 
both  stable  and  unstable  elements,  or  both,  as  illustrated  by 


(7)  Pg(s)  = , k - [k  ,k  ] 

s(s  +as+b)  1 


(1-13) 


and  the  complex  pole-pair  of  Pg(s)  ranging  in  MNQU  , M'N’Q’U'  of 


Figure  1-8 (a) . Three  typical  templates  at  different  frequencies 


wA<u)B<u>c  in  Figure  1-8  (a),  are  shown  in  Figures  1-8  (b) , . . . , (d) . 


Figure  1-8  (a)  Range  of  complex  pole-pair  of  Pg(s) 

(b) , (c) , (d)  Typical  templates  of  Pg(s)  , Case  (7). 


• * r~  ' Z '■ 
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APPENDIX  II 

Rational  Functions  of  Numerical  Examples 


The  following  expressions 

n(s)  - a s11  + a_sn  1 +....+  a ,s  + a 
1 z n-1  n 

d(8)  ■ b18m  + b-s10-^  +....+  b , s + b 
1 Z m-1  m 

represent  the  numerator  and  denominator  of  a rational  function 

R(s)  - £(s) 

d (s ) 

II.  1 Single  loop  system 

(1)  Figures  l-3(a),  2-1,  | AP | - lAPjP^  « 40  db 

Lg(s)  - 0.8533  E+09  x R(s)/s 

(ax  -*■  a?)  - 0. 10000000E+01,  0.  31934769E+02 , 0. 22690730E+03, 
0. 83455640E+03,  0. 15804 25 3E+04,  0 . 15747661E+04 , 
0. 65232178E+03. 

<bl  ^ bll>  " 0.10000000E+01,  0.37145410E+03,  0. 130521 8 7E+06, 
0. 18631760E+08,  0. 24542067E+10,  0. 38433726E+11 , 
0.202 72860E+12,  0. 65712659E+12 , 0. 11138572E+13, 


0. 10663364E+13 


0. 45332398E+12 


•fsrsz’rrrrrr  *— : — — *r*-*r;  - w 
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(2)  Figure  4-26,  lAP^P^  - 63  db 


Lg(s)  - 0. 2469E+13  x R(s)/s 


Uj.  ■*  a8)  - 0.10000000E+0X,  0.56891708E+02,  0. 53505396E+03, 
0.25187393E+04,  0.67933555E+04,  0. 10429051E+05, 
0. 86928672E+04,  0. 30754 795E+04. 


<bl  * b12> 


0. 10000000E+01 , 0.28619253E+04,  0. 80363180E+07, 
0.86478725E+10,  0.90781117E+13,  0. 23453811E+15, 
0. 15893765E+16,  0.63699863E+16,  0. 15216227E+17, 
0. 21976691E+17,  0. 17419455E+17,  0. 63546876E+16. 


II. 2 Cascaded  no  P.M.  system  (Figure  1-15 (a)) 


L1(s)  - 0.4513E+07  x R(s)/s 


(ax  -*•  a5)  - 0. 10000000E+01,  0. 12530328E+02 , 0.45847092E+02, 
0. 70532761E+02 , 0. 41426315E+02. 


0>1  -*■  bQ)  - 0. 10000000E+01 , 0.92955902E+02,  0.  78720469E+04, 
0.28583481E-H36,  0.91296180E+07,  0. 59384608E+08, 
0 . 18004706E+09 , 0. 22596626E+09,  0. 125876 78E+09 . 


L2(8)  - 2 


0.4  x 4Q2  x 2502 


(s‘c  + 0.7x40s  + 402)  (s2  + 0.6x250  s + 2502) 


II. 3 2- loop  P.M.  system  (Figures  2-2,  3-1) 

(1)  Design  A 


L (s)  - 0.1500E+08  x R(s)/s 
o 


1 


■ . 1 1 1 


r 


<*1  - a7) 
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0. 10000000 E+01 , 0. 2040963 7E+02,  0. 15927747E+03, 
0.67827539E+03,  0.16129690E+04,  0. 19764927E+04, 
0. 109582 35E+04. 


(bx  -*>  bn)  - 0. 10000000E+01,  0. 13237906E+03,  0. 13798289E405, 
0. 69865088E+06,  0.27376240E+08,  0. 36934195E+09, 
0. 24605775E+10,  0. 89491988E+10,  0. 19679154E+11, 
0. 20897321E+11,  0. 11361743E+11 


Lll<8)  ’ P!<8)hi(8) 


0.6(8  +2) (40s  + 1)  x 452  x 3002 


(s  + 7. 6)  (40s  + 16. 67)  (s2  + 45s  + 452)  (s2  + 0.  7 x 300s  + 3002) 


(2)  Design  B 


Lq(s)  - 0.6371E+07  x R(s)/s 


(ax  -*•  a5)  - 0. 10000000E+01 , 0.11242662E+02,  0.42154068E+02, 
0.67147141E+02,  0.40464294E+02. 

(bl  -*•  b9)  - 0 . 10000000E+01 , 0.99096939E+02,  0.86820664E+04, 
0. 33085312E+O6,  0. 10836 760E+08,  0. 71286944E+08, 
0.23882952E+09,  0. 32160640E+09,  0. 20816880E+09. 


Ln(s)  - 


0.68(1.25  + 3)  (100s  + l)(48s  + 1)  x 452  x 2502 

(1. 2s+7) (lOOs+3) (488+16. 67) ((1. 2s)2+l. 2x45s+452) (s2+0. 7x250s+2502) 
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(3)  Design  C 


Lq(s)  - 0.2174E+08  x R(s)/s 


(ax  -*■  a5)  - 0. 10000000E+01 , 0. 10000282E+02 , 0. 35892151E+02, 


0 . 56468002E+02 , 0. 33795563E+02. 


(bL  -*■  bg)  - 0. 10000000E+01,  0. 13268655E+03,  0. 14712281E+05, 
0. 72320231E+06,  0. 29226304E+08,  0. 20516770E+09, 
0.64369638E+09,  0. 89944986E+09,  0. 58078003E+09. 


Lil(8)  " 


0.5(s  + 3) (30s  + 1)  x 532  x 33Q2 

(s  + 9) (30s  + 25) (s2  + 53s  + 532) (s2  + 0. 7 x 330s  + 3302) 


(4)  Design  D 


Lq(s)  - 0.1062E+08  x R(s)/.s 


(ax  -»■  a5)  - 0 . 10000000E+01 , 0. 11204820 E+02 , 0.41860794E+02, 


0. 66633072E+02,  0.40388687E+02. 


(bx  -*■  bg)  - 0.10000000E+01,  0. 11282556E+03,  0. 11241965E+05, 
0. 48381150E+06,  0. 17928832E+08,  0. 11999072E+09, 
0. 39335680E+09 , 0. 54214912E+09,  0. 35039437E+09. 


L (s)  0. 5(s  + 2)  (65s  + 1)  x 452  x 3002 

11  (s  + 5. 5)  (65s  + 16.67)  (s2+45s  + 452)  (s2  + 0 . 7 x 300s  + 3002 ) 
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(5)  Design  E 


Lq(s)  - 0.4527E+08  x R(s)/s 


(ax  -*•  a4)  - 0. 10000000E+01,  0. 72563848E+01,  0.17384796E+02, 


0. 13756199E+02 . 


0>1  bg)  - 0 . 10000000E+01 , 0. 154594 35E+03,  0. 19788469E+05, 
0.11070540E+07,  0. 50771584E+08,  0.27372518E+09, 
0.59020390E+09,  0. 46127872E+09. 


L 0.5(s  + 4.5)  (16s  + 1)  x 602  x 4502 

11  (s  + 10) (16s  + 30) (s2  + 60s  + 602) (s2  + 0. 7 x 450s  + 4502) 


(6)  Design  ISLV 


LQ(s)  - 0. 4513E+07  x R(s)/s 


(ax  -*■  a5)  - 0.10000000E+01,  0.12530328E+02,  0.45847092E+02, 
0. 70532761E+02,  0. 41426315E+02. 


(bl  -►  bg)  - 0. 10000000E+01 , 0. 92955902E+02 , 0. 78720469E+04, 
0.28583481E+06,  0.91296180E+07,  0. 59384608E+08, 
0. 18004 706E+09,  0. 22596626E+09,  0. 12587678E+09. 


_ 0.4  x 402  x 2502 

11  (s2  + 0. 7 x 40s  + 402)  (s2  + 0. 6 x 250s  + 2502) 


4 3-loop  P.M.  system  (Figures  1-24,  4-1) 


Lq(s)  - 0.8331E+04  x R(s)/s 

(ax  -►  a4)  - 0. 10000000E+01 , 0. 66346006E+01,  0. 14799875E+02 , 
0. 11100864E+02 . 


<bl  * V 


0.10000000E+01 
0. 31348845E+04 
0. 10434300E+06 


0. 22659241E+02 
0. 19470348E+05 
0. 64274535E+05 


0. 35483521E+03 
0.66803438E+05 


Lll(8)  " 


4.125E+04  x (s  + 0.0153) 

(s2  + 7s  + 49)(s  + 0.92)(s2  + 0.4x  42s  + 422) 


L12(s)  " p2(s)H2(s) 


100  x 120^ 


s2  + 0.4  x 120s  + 1202 


x R(s) 


(a  -*•  a19)  - 0. 10000000E+01, 
0. 31943319E+06, 
0. 27825651E+10, 
0.30564963E+13, 
0. 48994571E+15, 
0. 90310750E+16, 


0. 14082841E+03,  0.83792813E+04, 
0. 83721440E+07,  0. 17572984E+09, 
0. 35667821E+11,  0. 36489101E+12 , 
0.20548717E+14,  0. 11304564E+15, 
0. 17120134E+16,  0. 44708156E+16, 
0. 11800826E+17,  0. 99326917E+16, 


0.15471337E+16 


<bl  - b21> 


0. 10000000E+01 , 0. 14528056E+03,  0. 946552 J4E+04 , 
0. 39388550E+06,  0. 115404 78E+08,  0. 26542619E+09, 
0.48147415E+10,  0. 70719570E+11,  0. 85274945E+12, 
0.85201749E+13,  0. 70939019E+14,  0.49207494E+15, 
0. 28397371E+16,  0.13527785E+17,  0. 52607418E+17, 
0.16437946E+18,  0. 39°02693E+18,  0. 73407265E+18, 
0.91156008E+18,  0. 70823165E+18,  0. 12667604E+18. 


II. 5.  5-loop  P.M.  system  (type  I)  (Figure  4-24) 


L0(s)  a®d  Lji(s)  are  the  same  In  3-loop  P.M.  system 


. /nx  . 2 x 82  x 602 

^02 222  o 

(s*  + 8s  + 8*)  (sz  + 0. 4 X 60s  + 6<T) 


(8) klIE^7(s^0  01)(s^  0.08)  x R(s) 

(0.  Is  + 1)  (10s  + 1)  (s^  + 0. 4 x 160s  + 160  ) 


(*i  - •„> 


- 0. lOOOOOOOE+Ol , 0. 21468179E+03,  0. 21766156E+05, 
0. 16546280E+07 , 0. 80722400E+08,  0. 29849879E+10, 
0.82077090E+11,  0. 16920483E+13,  0. 27214205E+14, 
0.31784369E+15,  0. 27451995E+16,  0. 15689971E+17 , 


0.55914070E+17,  0. 89149159E+17 


(b. 


b16) 


0 . 10000000E+01 , 
0.27528480E+07, 
0.28532185E+12, 
0.26835680E+16, 
0. 22259569E+19, 
0. 41061024E+20 


0.25456444E+03, 
0.17333214E+09, 
0. 77741834E+13, 
0.33811868E+17, 
0. 10466948E+20 , 


0 . 31895910E+05 , 
0 . 80436388E+10 , 
0. 16361775E+15, 
0. 32018864E+18, 
0. 30991538E+20, 


5.4E-H?9(s  180) 

(s  + 250)  (s2  + 2 x 0. 38  x 36s  + 362)  (s2  + 400s  + 106) 
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